First we need to build the HLM models in R (done) Then do the AIC testing (not possible. We’ll do a simple variance explained argument) Then do assumption checking Then build descriptives for gain aggregated and disagregated and add the means to the descriptives table

Load stuff (There are some packages that we don’t need here, but I didn’t clean them up. Sorry.)

#load("~/Box Sync/work/Research/R stuff/hmiout40") #Ben's computer
load("~/Documents/LASSO Data/Old/hmiout40")
library(tidyr)
library("ggplot2")
library(gvlma)
library("HLMdiag")
library("DHARMa")
library("car") #for the Levene test which we will not discuss here
library("Matrix")
library(mitools)
library(stargazer)
library(lme4)
library(nlme)
library(mice)
library(mitml)
library(multcomp)
library(foreach)
library(ggplot2)
library(stringr)
library(dplyr)  #I load dplyr last because some of its functions (select) will be masked by plyr and it is a PITA to debug

Creating mitml and extra variables

MIdata<-mids2mitml.list(hmi.out40) #converts file type

thing <- list()
for (i in 1:10){
  temp <- MIdata[[i]]
  class_means <- temp %>% group_by(crse_id) %>% summarise(pre_mean_class = mean(pre_scor))
  class_means$class_pre_cent <- mean(class_means$pre_mean_class) - class_means$pre_mean_class
  temp <- left_join(temp,class_means, by="crse_id")
  temp$stud_pre_cent <- temp$pre_mean_class - temp$pre_scor
  temp$stud_pre_grand <- temp$pre_scor-mean(temp$pre_scor)
  temp$gain <- temp$pst_scor - temp$pre_scor
  temp$collabnla <- ifelse(temp$colablrn==1,ifelse(temp$used_las==0,1,0),0)
  temp$coll <- temp$collabnla + temp$used_las
  #assign(df.names[i], temp)
  thing[[i]] <- temp
  }
MIdata <- as.mitml.list(thing)

Calculating the descriptive statistics

student_means <- foreach(i=1:10, .combine=rbind) %do% {
  temp <- MIdata[[i]]
  temp$instruction <- ifelse(temp$used_las==1,"Used_LAs",ifelse(temp$collabnla==1,"Collab_No_LA","Lecture"))
  student_means <- temp %>% group_by(instruction) %>% summarise(mean_gain = mean(gain))
  return <- student_means
  }

"Used_LAs"
## [1] "Used_LAs"
mean(student_means$mean_gain[student_means$instruction=="Used_LAs"])
## [1] 19.27524
"Collab_No_LA"
## [1] "Collab_No_LA"
mean(student_means$mean_gain[student_means$instruction=="Collab_No_LA"])
## [1] 24.999
"Lecture"
## [1] "Lecture"
mean(student_means$mean_gain[student_means$instruction=="Lecture"])
## [1] 12.00881
"All"
## [1] "All"
mean(student_means$mean_gain)
## [1] 18.76102
class_means <- foreach(i=1:10, .combine=rbind) %do% {
  temp <- MIdata[[i]]
  temp$instruction <- ifelse(temp$used_las==1,"Used_LAs",ifelse(temp$collabnla==1,"Collab_No_LA","Lecture"))

  class_means <- temp %>% group_by(crse_id) %>% summarise(gain = mean(gain))
  class_means <- left_join(class_means,unique(temp[c(13,50)]), by = "crse_id")
  class_means <- class_means %>% group_by(instruction) %>% summarise(gain = mean(gain))
  return <- class_means
}
 
"Used_LAs"
## [1] "Used_LAs"
mean(class_means$gain[class_means$instruction=="Used_LAs"])
## [1] 19.45648
"Collab_No_LA"
## [1] "Collab_No_LA"
mean(class_means$gain[class_means$instruction=="Collab_No_LA"])
## [1] 15.36021
"Lecture"
## [1] "Lecture"
mean(class_means$gain[class_means$instruction=="Lecture"])
## [1] 14.17819
 "All"
## [1] "All"
 mean(class_means$gain)
## [1] 16.33163
inst_means <- foreach(i=1:10, .combine=rbind) %do% {
  temp <- MIdata[[i]]
  class_means <- temp %>% group_by(crse_id) %>% summarise(gain = mean(gain))
  class_means <- left_join(class_means,unique(temp[c(13,14)]), by = "crse_id")
  class_means <- class_means %>% group_by(FMCE) %>% summarise(gain = mean(gain))
  return <- class_means
}
 
"FCI"
## [1] "FCI"
mean(inst_means$gain[inst_means$FMCE==0])
## [1] 17.48376
"FMCE"
## [1] "FMCE"
mean(inst_means$gain[inst_means$FMCE==1])
## [1] 19.12894
 inst_means_stu <- foreach(i=1:10, .combine=rbind) %do% {
  temp <- MIdata[[i]]
  class_means <- temp %>% group_by(FMCE) %>% summarise(gain = mean(gain))
  return <- class_means
}
 
"FCI"
## [1] "FCI"
mean(inst_means_stu$gain[inst_means_stu$FMCE==0])
## [1] 19.40524
"FMCE"
## [1] "FMCE"
mean(inst_means_stu$gain[inst_means_stu$FMCE==1])
## [1] 19.18769

Models (Note, we can’t actually use these since they don’t work with the with function)

mod1 <- (gain~1 + (1|crse_id))
mod2 <- (gain~1 + used_las + collabnla + (1|crse_id))
mod3 <- (gain~1 + stud_pre_cent+  used_las + collabnla + (1|crse_id))
mod4 <- (gain~1 + stud_pre_cent+  used_las + collabnla + (1+ stud_pre_cent|crse_id))
mod5 <- (gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent + (1|crse_id))
mod6 <- (gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent + (1+ stud_pre_cent|crse_id))
mod7 <- (gain~1 + stud_pre_cent+  used_las + collabnla + FMCE + (1|crse_id))
mod8 <- (gain~1 + stud_pre_cent+  used_las + collabnla + fst_time + (1|crse_id))

Run HLM models

HLM1<-with(MIdata,{lmer(gain~1 + (1|crse_id))})
HLM1.1<-with(MIdata,{lmer(gain~1 +  stud_pre_cent + (1|crse_id))})
HLM1.2<-with(MIdata,{lmer(gain~1 +  stud_pre_cent + class_pre_cent + (1|crse_id))})
HLM2<-with(MIdata,{lmer(gain~1 + used_las + collabnla + (1|crse_id))})
HLM3<-with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + (1|crse_id))})
HLM4<-with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + (1+ stud_pre_cent|crse_id))})
HLM5<-with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent + (1|crse_id))})
HLM6<-with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent + (1+ stud_pre_cent|crse_id))})
HLM7 <- with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + FMCE + (1|crse_id))})
HLM8 <- with(MIdata,{lmer(gain~1 + stud_pre_cent+  used_las + collabnla + fst_time + (1|crse_id))})

Run MLR models

MLR1<-with(MIdata,{lm(gain~1)})
MLR2<-with(MIdata,{lm(gain~1 + used_las + collabnla)})
MLR3<-with(MIdata,{lm(gain~1 + stud_pre_cent+  used_las + collabnla)})
MLR4<-with(MIdata,{lm(gain~1 + stud_pre_cent+  used_las + collabnla)})
MLR5<-with(MIdata,{lm(gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent)})
MLR6<-with(MIdata,{lm(gain~1 + stud_pre_cent+  used_las + collabnla + class_pre_cent)})

HLM Model outputs

testEstimates(HLM1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    18.284     0.872    20.977   669.237     0.000     0.131     0.119 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   60.620 
## Residual~~Residual            409.799 
## ICC|crse_id                     0.129 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    13.696     2.131     6.428   891.880     0.000     0.112     0.102 
## used_las        6.014     2.352     2.557  1215.587     0.011     0.094     0.088 
## collabnla       3.512     3.008     1.168   269.731     0.244     0.223     0.189 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   56.570 
## Residual~~Residual            409.830 
## ICC|crse_id                     0.121 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM3, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM3, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      13.760     2.144     6.419   812.980     0.000     0.118     0.107 
## stud_pre_cent     0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## used_las          5.936     2.363     2.512  1160.173     0.012     0.097     0.090 
## collabnla         3.256     3.023     1.077   247.880     0.283     0.235     0.197 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   58.827 
## Residual~~Residual            343.497 
## ICC|crse_id                     0.146 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM4, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM4, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      13.347     2.052     6.505  1402.786     0.000     0.087     0.081 
## stud_pre_cent     0.426     0.025    16.833    33.212     0.000     1.086     0.547 
## used_las          6.662     2.289     2.911   951.492     0.004     0.108     0.099 
## collabnla         2.952     2.951     1.001   201.219     0.318     0.268     0.219 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           58.825 
## Intercept~~stud_pre_cent|crse_id        0.437 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.331 
## ICC|crse_id                             0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM5, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM5, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       13.659     2.209     6.184  1545.906     0.000     0.083     0.077 
## stud_pre_cent      0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## used_las           6.036     2.445     2.469  1426.050     0.014     0.086     0.081 
## collabnla          3.388     3.135     1.081   306.195     0.281     0.207     0.177 
## class_pre_cent     0.015     0.098     0.158   337.055     0.874     0.195     0.168 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   59.394 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM6, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM6, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       12.932     2.122     6.093  2448.080     0.000     0.065     0.061 
## stud_pre_cent      0.426     0.025    16.876    33.600     0.000     1.073     0.544 
## used_las           7.091     2.405     2.948   636.817     0.003     0.135     0.122 
## collabnla          3.493     3.032     1.152   279.923     0.250     0.218     0.185 
## class_pre_cent     0.060     0.096     0.623   218.282     0.534     0.255     0.210 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           59.644 
## Intercept~~stud_pre_cent|crse_id        0.455 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.323 
## ICC|crse_id                             0.149 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM7, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM7, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      13.715     2.160     6.349   769.754     0.000     0.121     0.110 
## stud_pre_cent     0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## used_las          5.789     2.398     2.414  1269.656     0.016     0.092     0.086 
## collabnla         3.194     3.031     1.054   259.469     0.293     0.229     0.192 
## FMCE              0.821     2.153     0.381  4056.408     0.703     0.049     0.048 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   59.400 
## Residual~~Residual            343.493 
## ICC|crse_id                     0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM8, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM8, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      10.995     2.169     5.069   943.235     0.000     0.108     0.100 
## stud_pre_cent     0.436     0.023    18.717    19.150     0.000     2.180     0.714 
## used_las          5.172     2.351     2.200   944.587     0.028     0.108     0.100 
## collabnla         2.467     3.007     0.820   227.714     0.413     0.248     0.206 
## fst_time          4.244     0.892     4.756    88.991     0.000     0.466     0.333 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   57.251 
## Residual~~Residual            341.744 
## ICC|crse_id                     0.143 
## 
## Unadjusted hypothesis test as appropriate in larger samples.

MLR Model outputs

testEstimates(MLR1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    19.337     0.343    56.366    81.196     0.000     0.499     0.349 
## 
##                    Estimate 
## Residual~~Residual  467.805 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(MLR2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    12.009     0.993    12.089    51.731     0.000     0.716     0.438 
## used_las        7.266     1.063     6.836    58.355     0.000     0.647     0.413 
## collabnla      12.990     1.340     9.694    45.976     0.000     0.794     0.465 
## 
##                    Estimate 
## Residual~~Residual  454.978 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(MLR3, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR3, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      12.009     0.951    12.627    43.459     0.000     0.835     0.479 
## stud_pre_cent     0.440     0.024    18.456    20.828     0.000     1.918     0.686 
## used_las          7.266     1.016     7.154    48.651     0.000     0.755     0.452 
## collabnla        12.990     1.285    10.106    38.930     0.000     0.926     0.506 
## 
##                    Estimate 
## Residual~~Residual  389.853 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(MLR4, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR4, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      12.009     0.951    12.627    43.459     0.000     0.835     0.479 
## stud_pre_cent     0.440     0.024    18.456    20.828     0.000     1.918     0.686 
## used_las          7.266     1.016     7.154    48.651     0.000     0.755     0.452 
## collabnla        12.990     1.285    10.106    38.930     0.000     0.926     0.506 
## 
##                    Estimate 
## Residual~~Residual  389.853 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(MLR5, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR5, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       12.392     1.014    12.226    41.175     0.000     0.878     0.492 
## stud_pre_cent      0.440     0.024    18.457    20.823     0.000     1.919     0.686 
## used_las           6.915     1.077     6.421    43.681     0.000     0.831     0.477 
## collabnla         12.427     1.408     8.827    34.829     0.000     1.034     0.534 
## class_pre_cent    -0.050     0.035    -1.435   160.813     0.153     0.310     0.246 
## 
##                    Estimate 
## Residual~~Residual  389.726 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(MLR6, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = MLR6, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       12.392     1.014    12.226    41.175     0.000     0.878     0.492 
## stud_pre_cent      0.440     0.024    18.457    20.823     0.000     1.919     0.686 
## used_las           6.915     1.077     6.421    43.681     0.000     0.831     0.477 
## collabnla         12.427     1.408     8.827    34.829     0.000     1.034     0.534 
## class_pre_cent    -0.050     0.035    -1.435   160.813     0.153     0.310     0.246 
## 
##                    Estimate 
## Residual~~Residual  389.726 
## 
## Unadjusted hypothesis test as appropriate in larger samples.

Assumption checking as recommended here: https://ademos.people.uic.edu/Chapter18.html Note: I think we may just need to do this for each of the 10 datasets independently.

for (i in 1:10) {
  D1 <- lmer(mod3,data=MIdata[[i]])
  print(plot(D1, xlab="Fitted Value", ylab="Residual Variance"))
}

for (i in 1:10) {
  D1 <- lmer(mod3,data=MIdata[[i]])
#visual homogeneity of variance
MIdata[[i]]$Model.F.Res <- residuals(D1) #extracts the residuals and places them in a new column in our original data table
MIdata[[i]]$Abs.Model.F.Res <-abs(MIdata[[1]]$Model.F.Res) #creates a new column with the absolute value of the residuals
print(boxplot(MIdata[[i]]$Model.F.Res ~ MIdata[[i]]$crse_id), xlab = "Course", ylab = "Residuals" )
}

## $stats
##            [,1]       [,2]         [,3]       [,4]       [,5]       [,6]
## [1,] -46.626554 -43.396099 -33.36069090 -21.969755 -28.883951 -21.717970
## [2,] -13.520209 -16.237893  -8.52206222 -14.101199  -3.042342 -11.276046
## [3,]  -1.573926   1.948848   0.08743497  -6.020259   6.624863  -7.308162
## [4,]  14.627506  17.648130   8.13088313   6.961936  14.904044  14.096205
## [5,]  55.179459  36.285938  21.93974961  36.666446  26.234585  40.110791
##            [,7]      [,8]       [,9]       [,10]       [,11]     [,12]
## [1,] -31.857444 -8.892331 -42.181440 -40.0885979 -26.6327687 -42.26774
## [2,] -13.655341 -6.389056 -11.048939 -15.9526662 -15.5697580 -11.86983
## [3,]  -2.073521 -4.542794   2.661469  -0.6394002   0.2183416   2.57054
## [4,]  13.421639 -2.572436  11.309141  13.8285954  12.2217939  13.14444
## [5,]  27.506707 -2.572436  42.573088  43.9801944  18.2664497  40.92921
##           [,13]      [,14]      [,15]      [,16]      [,17]         [,18]
## [1,] -42.559275 -39.293224 -32.165401 -22.512412 -34.247775 -35.844352138
## [2,] -10.024324  -8.474646  -8.716382  -7.901208 -14.165349 -14.692055349
## [3,]   1.422061   3.161244   4.060702   2.649683  -5.595163   0.009651716
## [4,]  12.031453  12.436845  15.241085   7.781019   2.648998  13.608987251
## [5,]  37.336445  40.380257  30.519554  17.928337  27.767054  31.397849738
##           [,19]      [,20]      [,21]     [,22]       [,23]      [,24]
## [1,] -26.424563 -30.100529 -24.827803 -8.748752 -42.7566919 -40.909356
## [2,] -16.128861 -15.034273 -17.570945 -3.147701 -14.6240408 -12.547129
## [3,]  -7.026092  -6.511989  -3.647421  1.134874  -0.8115566   2.116302
## [4,]  -2.674471  12.196615  10.200372  5.324401  14.0868299  13.726306
## [5,]   9.610885  43.476429  21.727449 17.139676  32.5217768  49.859444
##           [,25]      [,26]      [,27]      [,28]      [,29]      [,30]
## [1,] -43.883248 -26.872969 -52.238275 -28.406176 -19.902530 -26.672924
## [2,] -12.727889  -6.955277 -15.279503  -8.256906  -8.238086  -7.911494
## [3,]   4.319961  -2.390544   3.869433   2.194991   3.810957  -1.465109
## [4,]  17.192103   9.080970  12.587000  11.107540  12.572387   6.018462
## [5,]  36.856481  22.604159  41.413869  36.431729  18.856678  21.822504
##           [,31]      [,32]       [,33]      [,34]       [,35]     [,36]
## [1,] -37.270857 -28.624540 -35.9660908 -33.410471 -24.9617863 -8.403255
## [2,] -22.075291 -12.528352  -8.7058514 -14.767647 -13.0443015 -8.403255
## [3,]  -3.682146  -3.177811   0.6281689  -1.226499  -6.8971571 -8.208694
## [4,]  13.227708   3.011335  10.3763866   9.421880  -0.9207281 -8.014134
## [5,]  57.858280  23.327566  34.3580984  37.917308  -0.2996140 -8.014134
##           [,37]       [,38]      [,39]        [,40]      [,41]      [,42]
## [1,] -31.788802 -35.8805053 -26.017221 -39.42191453 -49.138826 -29.157211
## [2,] -31.788802  -7.3647122 -10.250248  -9.68431176 -14.359562 -10.154358
## [3,] -16.444698  -0.4807583  -2.955112   0.04608107   1.068337  -5.380110
## [4,]  -1.100594  14.0626148   6.720699  10.47417788  11.982916   7.339124
## [5,]  -1.100594  33.8708632  20.577891  34.55979724  40.744184  18.424869
##           [,43]      [,44]       [,45]      [,46]      [,47]       [,48]
## [1,] -39.684356 -30.919936 -46.7859023 -30.150428 -29.142735 -24.8496443
## [2,] -15.823320 -13.513450 -16.5020183 -12.083097 -13.924395  -8.6424592
## [3,]  -2.792950  -6.637725  -0.7600222   1.678333  -3.597413  -0.4712205
## [4,]   7.858985  13.224712  10.4524598   9.201194  11.402587   8.1500511
## [5,]  25.151988  42.315088  36.6281539  30.409801  36.704288  32.1433677
##           [,49]       [,50]      [,51]      [,52]      [,53]     [,54]
## [1,] -26.868973 -27.3884471 -34.694260 -25.203893 -35.172903 -28.32820
## [2,] -11.837412  -7.0642579 -10.755183 -10.365987  -7.306399 -11.83751
## [3,]  -2.103414  -0.6178729  -1.739826  -5.528082   3.027089   3.28280
## [4,]  10.315539  11.2345725  14.810467  11.406547  13.519198  11.62511
## [5,]  29.792948  31.1069808  28.908552  28.777819  29.824179  35.75918
##           [,55]      [,56]      [,57]     [,58]       [,59]       [,60]
## [1,] -37.008193 -30.268685 -21.342626 -57.37577 -42.4161435 -37.2675414
## [2,] -17.685391  -7.357643  -3.984895 -17.25939 -12.5563153 -13.0151294
## [3,]   1.420402   1.993979  13.930423  -3.27768  -0.5197389   0.1998605
## [4,]  13.083711  12.318168  21.642951  15.21775   8.6766653  10.8404138
## [5,]  42.019702  31.022026  36.345354  55.21775  36.4711170  36.8967705
##          [,61]       [,62]      [,63]      [,64]       [,65]      [,66]
## [1,]  1.419309 -17.9261291 -20.269998 -41.825387 -46.1236186 -29.273815
## [2,]  8.490438 -10.6036640 -11.566755 -17.941359 -11.0190234  -7.511820
## [3,] 15.561567  -0.5071764   0.719256  -4.549503   0.3513302  -0.305641
## [4,] 17.091123   9.9980510   7.252863  13.623153  13.8442773  12.447076
## [5,] 18.620680  32.1535702  31.558764  46.514326  38.9273157  42.111608
##            [,67]       [,68]      [,69]      [,70]      [,71]      [,72]
## [1,] -41.6935416 -26.3520979 -38.311558 -25.925445 -37.653727 -20.153846
## [2,] -12.5102652  -9.5873401 -11.215698 -10.399791  -9.871532 -11.961448
## [3,]  -0.8301939  -0.7677229  -1.769326  -1.255984  -1.152704  -1.010039
## [4,]  14.5538379   6.2231330  11.127418   6.915254   9.071519   9.139589
## [5,]  36.8064607  29.4383008  36.138227  29.018491  29.258860  26.714253
##           [,73]     [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -27.952161 -29.72295 -29.706559 -26.893885 -37.997025 -34.755194
## [2,]  -9.332629 -11.95539 -19.174555 -12.736553 -14.356589 -15.270747
## [3,]  -1.335305  -3.84214   2.213643  -3.561891  -3.752743  -7.988033
## [4,]   6.438526  11.46939  13.167808   6.438109   9.790961  20.129432
## [5,]  25.388903  35.67280  39.401815  33.960969  45.391064  38.057338
##            [,79]      [,80]       [,81]      [,82]     [,83]       [,84]
## [1,] -38.9365605 -39.772429 -33.9692047 -36.282706 -51.75386 -34.4734284
## [2,] -10.7020574 -11.718850  -9.2352069 -12.551399 -12.69512  -9.6355230
## [3,]   0.3140579  -3.060842   0.3366963  -3.199777  -1.09275   0.8964814
## [4,]  12.5986738  14.428602  12.8595565  10.989750  15.93255  17.2389589
## [5,]  40.1541363  40.258898  23.0216512  46.268218  46.44868  41.2955901
##            [,85]      [,86]      [,87]      [,88]      [,89]      [,90]
## [1,] -33.7366598 -49.762474 -39.036757 -34.826853 -38.306507 -17.685949
## [2,] -13.5745652 -11.310026 -11.206999 -18.555131 -12.329683  -8.600330
## [3,]   0.4528671   3.151824  -1.189708   1.292147   2.225497  -6.181377
## [4,]  13.1378220  14.694114  10.272886  13.262224  12.283684   2.314051
## [5,]  39.9208627  38.379580  37.269484  44.661961  35.060574  17.893954
##           [,91]      [,92]       [,93]      [,94]         [,95]      [,96]
## [1,] -47.189471 -34.337573 -35.1600525 -18.073295 -35.339907641 -33.187922
## [2,] -13.155681 -16.539882 -14.6695529 -13.686154  -9.837521665 -11.513578
## [3,]   2.131844  -2.558049   0.6538721  -6.524558   0.004955097  -2.183469
## [4,]  14.159364  11.307168   9.5157582   3.755466   7.552862433  17.284527
## [5,]  34.691369  49.802596  37.8490916  18.874537  31.858763429  47.608716
##            [,97]      [,98]       [,99]     [,100]      [,101]     [,102]
## [1,] -38.6194230 -21.691743 -32.2820470 -31.546337 -38.0651415 -24.985920
## [2,] -13.4573284 -12.975173  -8.4070609 -16.058832 -12.2998687 -12.575071
## [3,]   0.4080093  -4.994653   0.4029079  -5.081117  -0.3780761  -5.584215
## [4,]   6.2575722   7.214648   9.2225251  10.718356  12.3703124   9.577879
## [5,]  29.8210431  37.358757  33.9187294  50.394167  38.0825015  30.708081
##         [,103]     [,104]    [,105]     [,106]     [,107]     [,108]
## [1,] -10.65423 -37.317623 -40.27863 -34.805996 -39.651942 -59.628197
## [2,] -10.65423 -11.145389 -10.39500 -14.208555 -13.334119 -17.107621
## [3,] -10.65423  -2.974150   2.52270  -3.503933  -2.396116  -0.574484
## [4,] -10.65423   9.635768  12.93833  10.909502  11.783937  17.258595
## [5,] -10.65423  36.053283  44.35146  36.987923  41.549667  53.651652
##           [,109]     [,110]     [,111]      [,112]
## [1,] -61.4367683 -32.963608 -29.020863 -44.5982572
## [2,] -15.6931455 -13.788579 -12.620374 -16.2935775
## [3,]  -0.2574137  -1.459036  -0.533571  -0.6104715
## [4,]  17.3098087  15.648194  10.104727  14.6529937
## [5,]  50.9938856  32.963239  32.888194  54.2295824
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]      [,4]     [,5]       [,6]
## [1,] -4.506415 -3.892836 -6.490478 -13.46184 -1.23948 -17.330201
## [2,]  1.358562  7.790532  6.665348   1.42132 14.48921   2.713878
##            [,7]      [,8]       [,9]     [,10]     [,11]     [,12]
## [1,] -13.939010 -7.239609 -0.2325302 -7.228333 -10.75932 -2.222273
## [2,]   9.791968 -1.845980  5.5554688  5.949533  11.19600  7.363353
##           [,13]    [,14]     [,15]     [,16]       [,17]     [,18]
## [1,] -0.7775337 1.063192 -2.968386 -1.800563 -10.4455513 -9.117892
## [2,]  3.6216549 5.259296 11.089790  7.099928  -0.7447754  9.137195
##           [,19]     [,20]     [,21]     [,22]     [,23]     [,24]
## [1,] -12.514868 -15.68491 -20.23200 -2.321356 -8.708272 -4.141871
## [2,]  -1.537317   2.66093  12.93716  4.591103  7.085159  8.374476
##          [,25]     [,26]     [,27]     [,28]     [,29]     [,30]
## [1,] -3.348827 -7.792469 -3.006754 -2.417503 -3.045111 -5.133330
## [2,] 11.988750  3.011380 10.745619  6.807485 10.667025  2.203113
##           [,31]     [,32]     [,33]     [,34]      [,35]     [,36]
## [1,] -13.700304 -7.902981 -4.876431 -7.509738 -12.672681 -8.643431
## [2,]   6.336011  1.547360  6.132769  5.056741  -1.121633 -7.773957
##          [,37]     [,38]     [,39]     [,40]     [,41]      [,42]
## [1,] -50.73045 -7.391417 -8.806428 -3.349184 -2.699862 -11.721092
## [2,]  17.84105  6.429901  2.896203  3.441347  4.836536   0.960872
##          [,43]       [,44]     [,45]     [,46]     [,47]     [,48]
## [1,] -9.624518 -13.9918617 -8.288601 -5.333835 -9.146721 -3.419239
## [2,]  4.038619   0.7164107  6.768557  8.690502  1.951895  2.476798
##          [,49]     [,50]     [,51]        [,52]     [,53]     [,54]
## [1,] -7.637663 -5.576266 -8.207999 -11.10859702 -2.241839 -2.506708
## [2,]  3.430835  4.340520  4.728347   0.05243291  8.296016  9.072309
##          [,55]     [,56]     [,57]      [,58]     [,59]     [,60]
## [1,] -6.917032 -1.548802 -1.374113 -13.540456 -5.466138 -5.240476
## [2,]  9.757837  5.536760 29.234959   6.985096  4.426660  5.640197
##          [,61]     [,62]     [,63]      [,64]     [,65]     [,66]
## [1,]  7.715906 -9.903757 -6.289350 -13.231092 -2.866942 -5.880304
## [2,] 23.407228  8.889404  7.727862   4.132085  3.569602  5.269022
##          [,67]     [,68]     [,69]     [,70]     [,71]     [,72]     [,73]
## [1,] -7.428403 -3.614520 -5.490497 -4.295737 -5.117032 -7.961834 -5.049924
## [2,]  5.768015  2.079075  1.951844  1.783768  2.811623  5.941756  2.379315
##          [,74]     [,75]     [,76]      [,77]      [,78]     [,79]
## [1,] -10.59942 -10.18015 -7.981012 -10.201799 -20.193458 -3.854438
## [2,]   2.91514  14.60744  0.857230   2.696313   4.217391  4.482553
##          [,80]     [,81]     [,82]     [,83]     [,84]     [,85]
## [1,] -8.070776 -4.755411 -8.744483 -8.072157 -7.595855 -5.638982
## [2,]  1.949092  5.428803  2.344928  5.886657  9.388818  6.544717
##             [,86]     [,87]     [,88]     [,89]     [,90]     [,91]
## [1,] 0.0006301753 -6.140105 -6.656462 -4.001746 -9.096267 -4.869271
## [2,] 6.3030179362  3.760690  9.240757  8.452740 -3.266487  9.132959
##          [,92]     [,93]      [,94]     [,95]     [,96]     [,97]
## [1,] -9.791331 -6.700183 -16.267698 -2.503323 -10.10418 -4.579909
## [2,]  4.675234  8.007927   3.218581  2.513233   5.73724  5.395928
##          [,98]     [,99]     [,100]    [,101]      [,102]    [,103]
## [1,] -14.20337 -2.108670 -11.319084 -7.268635 -11.3384531 -10.65423
## [2,]   4.21406  2.914486   1.156851  6.512483   0.1700224 -10.65423
##          [,104]    [,105]     [,106]     [,107]    [,108]    [,109]
## [1,] -6.9269263 -2.197592 -7.3586291 -4.6009232 -3.683614 -5.206767
## [2,]  0.9786267  7.242993  0.3507638 -0.1913089  2.534646  4.691939
##         [,110]    [,111]    [,112]
## [1,] -9.105250 -8.188683 -5.315451
## [2,]  6.187177  7.121541  4.094508
## 
## $out
##  [1]  43.35023  20.62217 -81.61720 -45.28892  48.18595 -46.95699 -43.86278
##  [8] -56.50827 -42.23121 -40.40180 -53.24530 -42.18599  44.07973 -48.46708
## [15] -46.51414 -48.82447 -44.01476  35.87725  35.65883 -52.26692  32.22930
## [22]  41.90724  46.91599  69.64057 -23.47768  61.01720  40.50768 -38.25203
## [29] -39.82312  46.13059  40.31790 -37.61457  66.84112  43.02301  41.34895
## [36]  46.08855  17.75525  21.62059  39.72984  34.89194  43.40853  48.50736
## [43] -38.48917  38.55621 -38.14054 -36.20058 -40.18978 -38.36101  57.19870
## [50] -47.98755 -40.17571  52.09177 -52.62743  53.42180  33.24104 -39.57326
## [57]  59.16560  33.55384 -38.48184 -54.32147  30.99156  38.97455 -56.99149
## [64] -54.03037  48.42117 -65.06394 -50.70330  24.67482  37.80033 -32.73530
## [71]  46.01526 -47.45532 -39.38558  40.67838 -38.14124  40.95347  44.18191
## [78] -40.70771 -34.92861  50.35219  59.21378 -69.73165  54.48010 -52.11548
## [85] -54.90872  63.69304  61.79909  68.83895 -54.96996  55.84400  55.84400
## 
## $group
##  [1]   3   8   9   9   9  13  13  14  14  14  14  14  14  14  14  14  14
## [18]  16  16  16  17  17  17  19  22  23  26  26  29  30  30  32  32  32
## [35]  33  34  35  35  39  39  40  43  48  48  48  48  50  50  52  53  56
## [52]  65  65  68  68  68  70  70  70  73  73  76  81  82  82  86  86  90
## [69]  90  90  90  94  95  95  95  95  97  99  99  99 100 101 101 105 105
## [86] 107 107 107 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##            [,1]        [,2]       [,3]       [,4]       [,5]       [,6]
## [1,] -45.691977 -41.0858813 -20.252348 -27.700243 -15.154509 -21.177193
## [2,] -13.527528 -14.5322277 -10.015755 -17.178699  -8.227087 -10.927071
## [3,]  -1.169404  -0.0777659  -3.450004  -4.505999  -1.170727  -5.926246
## [4,]  12.312663  16.3715115   3.623988  11.275797   9.452587  14.869563
## [5,]  47.094883  36.4309821   9.100288  33.142198  35.085012  40.785866
##            [,7]       [,8]       [,9]       [,10]      [,11]      [,12]
## [1,] -26.465156 -10.324073 -36.754741 -38.6876348 -34.067382 -33.711069
## [2,] -17.212440  -9.654698  -9.400722 -12.9074557 -18.570017  -8.446475
## [3,]   1.160567  -3.847626   1.526856  -0.2426152   1.634022   1.979811
## [4,]  11.152251   8.069651  10.362382  13.8875417  13.220539  14.011838
## [5,]  29.762984  15.235619  39.417551  32.1285221  38.102958  30.648602
##           [,13]      [,14]      [,15]       [,16]      [,17]      [,18]
## [1,] -33.293373 -37.341830 -30.590587 -27.8083760 -36.744135 -35.554617
## [2,]  -7.972197  -7.845016 -12.553646  -7.5409657 -16.525858 -15.168963
## [3,]   1.510902   4.176251   3.889648   0.5272927  -5.084206  -2.207828
## [4,]   9.159098  11.833741  15.852707   7.0824899  21.133166  11.845015
## [5,]  33.289255  32.191854  31.557177  20.3457591  40.282102  44.959693
##          [,19]      [,20]      [,21]      [,22]       [,23]      [,24]
## [1,] -26.16146 -28.270556 -21.309987 -10.546906 -44.4397365 -35.575324
## [2,] -17.99430 -13.991107 -18.039597  -3.933998 -14.2163640  -9.676759
## [3,] -13.63572  -6.918341  -9.939713  -1.389777  -0.8830307   2.658708
## [4,]   3.46716   9.481684   5.892376   3.333164  15.3755576  12.962048
## [5,]  27.80965  31.974663  19.467503   7.820041  51.9123202  32.276911
##           [,25]       [,26]       [,27]      [,28]      [,29]      [,30]
## [1,] -29.739220 -45.2484884 -35.8653936 -26.561069 -17.041896 -28.814805
## [2,]  -8.993484 -14.9659110 -15.6389631 -12.461420  -6.745504  -7.104216
## [3,]  -1.099006   0.1137525  -0.7453461   3.610635   5.105869   3.570824
## [4,]  15.030947   9.5264982  15.7549793  12.016988  10.883359   8.929080
## [5,]  34.933201  33.3869808  44.6922509  39.148834  22.399079  26.760877
##           [,31]      [,32]      [,33]      [,34]      [,35]      [,36]
## [1,] -31.781245 -29.915602 -34.937359 -31.152441 -28.603701 -26.494137
## [2,] -19.700745 -10.518531 -10.233751 -16.225461 -19.787929 -26.494137
## [3,]  -7.467195  -5.258331   1.432916  -2.299343  -5.800774 -18.151606
## [4,]   8.231393  12.120157  16.248210   9.255637  11.404537  -9.809074
## [5,]  49.187716  34.032082  42.914877  45.922304  33.656326  -9.809074
##           [,37]      [,38]      [,39]       [,40]      [,41]     [,42]
## [1,] -57.001543 -24.914682 -20.526702 -30.9223489 -48.340969 -20.66644
## [2,] -57.001543  -9.243757  -8.044708  -9.2750155 -12.599681 -11.48615
## [3,] -27.882780   3.461279  -2.340237  -0.7119329  -1.311697  -1.75078
## [4,]   1.235983  11.391157   1.951006   9.5329644  13.882280   4.33053
## [5,]   1.235983  39.717839  11.547214  35.3561573  43.434745  21.42734
##           [,43]      [,44]      [,45]      [,46]      [,47]      [,48]
## [1,] -47.426714 -33.389792 -44.654973 -29.132079 -33.700218 -31.847318
## [2,] -16.031842 -11.703626 -12.254318 -10.594707  -9.008298  -9.885426
## [3,]  -2.029511  -5.681429   1.579124   1.645411  -1.042152  -2.796003
## [4,]  10.025901   7.417973  13.673816   8.219724   7.615828   8.929006
## [5,]  48.470920  24.677424  33.618843  23.423764  26.438992  36.373163
##            [,49]      [,50]       [,51]      [,52]      [,53]      [,54]
## [1,] -36.6593963 -40.240341 -31.8684828 -27.991512 -32.853222 -31.570084
## [2,] -12.6522231 -11.833988  -8.9238946  -9.227135  -5.630969 -11.095545
## [3,]  -0.7815172   1.735207  -0.9062867  -3.002750   1.400445   2.787135
## [4,]  11.3662475  10.537149  10.2086817   9.402590  13.156240  14.388464
## [5,]  24.4031887  32.130796  29.3418002  29.995375  32.229582  37.458920
##           [,55]       [,56]      [,57]      [,58]       [,59]       [,60]
## [1,] -29.572213 -37.8396044 -19.029817 -28.719279 -34.8924269 -45.1010917
## [2,] -17.025729  -9.2245848 -17.955699 -15.794024 -10.7816034 -13.0395647
## [3,]  -0.480877  -0.3313283   3.282823  -2.052612   0.3034766   0.8037111
## [4,]  12.750898  12.2952322   6.324332  13.243780   6.6625535  13.3011076
## [5,]  45.530114  37.8679470  22.952711  54.323500  29.6264751  31.9159023
##           [,61]      [,62]      [,63]      [,64]      [,65]      [,66]
## [1,] -28.718717 -22.165686 -19.732764 -33.649104 -39.541993 -23.881044
## [2,] -17.770646 -13.569788 -12.103902 -13.956414 -11.405304  -9.992580
## [3,]  -6.822575   1.769090   2.769823  -2.694164   1.750305   1.333668
## [4,]  -2.325858   7.872343   8.225898  13.061979  13.614613  10.526026
## [5,]   2.170859  23.321572  32.230294  41.315252  37.943807  29.306708
##          [,67]      [,68]      [,69]      [,70]      [,71]      [,72]
## [1,] -30.55967 -28.896442 -32.524048 -27.885999 -27.151557 -21.806688
## [2,] -11.64477  -9.550632 -10.317432 -10.712469 -11.881788 -12.676531
## [3,]  -2.77050  -1.701827  -1.687706   0.473100  -3.067357   2.119430
## [4,]  12.19248   7.297310   8.905078   7.732551   9.488486   7.785233
## [5,]  43.29916  30.731660  36.655563  29.422868  39.966721  26.936352
##           [,73]      [,74]      [,75]      [,76]       [,77]      [,78]
## [1,] -24.793670 -27.137831 -24.730254 -19.195786 -32.3723293 -26.308868
## [2,]  -9.312572 -11.841439 -15.138332 -10.801643 -15.7940344 -15.149724
## [3,]  -2.053121  -2.084781  -7.696179  -4.105424  -0.5783387  -7.917934
## [4,]   8.514449  12.236988  14.577945   5.272239  12.9396245   6.862896
## [5,]  27.169389  32.547549  34.603185  25.272239  43.9405633  30.305373
##           [,79]     [,80]       [,81]      [,82]     [,83]      [,84]
## [1,] -41.415797 -40.76603 -34.6103099 -37.663955 -48.65151 -36.854474
## [2,] -11.026487 -12.68871  -9.4062707 -12.571647 -11.90229 -16.280993
## [3,]  -1.026487  -1.20387   0.4404565  -2.359137   1.55870  -2.636299
## [4,]  14.269905  14.47929  12.3527489  11.480164  15.69036  19.238771
## [5,]  37.010454  32.90359  25.0390709  46.650863  46.95080  41.692318
##            [,85]      [,86]      [,87]     [,88]      [,89]      [,90]
## [1,] -32.9152136 -49.075832 -38.687988 -31.45333 -36.956699 -16.722260
## [2,] -12.6188215 -12.018671 -10.855086 -16.14493 -12.556501  -7.018652
## [3,]  -0.9521548   2.400604   3.459776  -4.18432   2.635222  -5.351985
## [4,]  12.5126920  13.984834   9.713934  12.33682  12.654556   3.277740
## [5,]  41.0109040  40.027874  38.993183  48.70709  35.006360  16.202995
##           [,91]      [,92]      [,93]      [,94]       [,95]      [,96]
## [1,] -29.900979 -32.512355 -36.872125 -23.036804 -35.8856385 -40.070919
## [2,] -13.126213 -15.161556 -15.266822 -12.355630  -9.3653903 -11.263561
## [3,]   3.519687  -2.421007  -4.377646  -2.180559   0.5229234  -2.179875
## [4,]  14.726022  13.244797  12.158002   5.871549   9.4490410  11.820695
## [5,]  40.494125  50.912326  39.140394  25.115843  31.8008450  37.635419
##           [,97]      [,98]      [,99]     [,100]      [,101]     [,102]
## [1,] -34.051047 -26.773364 -32.255661 -37.672697 -47.5658342 -36.240996
## [2,] -12.192336 -12.033178  -9.350255 -15.866234 -14.7236859 -12.560573
## [3,]  -1.495204  -1.309566   0.261515  -9.087881   0.1426011  -4.115554
## [4,]   5.736242   8.389123   9.076891  13.171570  12.8713907  11.007183
## [5,]  21.060639  16.246277  29.892103  52.578786  50.0949047  32.926698
##         [,103]     [,104]    [,105]    [,106]     [,107]     [,108]
## [1,] -11.83354 -41.026611 -43.51939 -34.87249 -48.013605 -62.617625
## [2,] -11.83354 -10.795734 -10.97033 -13.07494 -12.923798 -16.399750
## [3,] -11.83354  -2.758793   2.77381  -3.81990  -3.302721  -1.494751
## [4,] -11.83354  10.389834  11.36214  11.70022  13.115292  17.117163
## [5,] -11.83354  35.462854  32.95579  37.35004  43.246230  66.899907
##           [,109]      [,110]     [,111]      [,112]
## [1,] -56.3308122 -36.0724871 -33.061180 -46.7782603
## [2,] -15.4844781  -7.8121734 -11.393504 -17.8628372
## [3,]  -0.3481338  -0.9608008  -1.688798  -0.3661955
## [4,]  17.6533430  14.1508854  11.018210  14.7486259
## [5,]  49.0843067  32.0031207  34.800839  53.7655766
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]       [,4]      [,5]       [,6]
## [1,] -3.861490 -5.405327 -8.837702 -14.558936 -8.918193 -16.115917
## [2,]  1.522681  5.249796  1.937694   5.546939  6.576739   4.263424
##           [,7]       [,8]      [,9]     [,10]     [,11]     [,12]
## [1,] -11.26921 -16.371606 -1.031253 -6.170854 -10.92325 -2.323272
## [2,]  13.59035   8.676355  4.084966  5.685624  14.19129  6.282894
##           [,13]    [,14]     [,15]     [,16]      [,17]      [,18]
## [1,] -0.1975799 2.201879 -4.444737 -3.622499 -15.947599 -10.920272
## [2,]  3.2193842 6.150622 12.224033  4.677084   5.779188   6.504617
##           [,19]       [,20]      [,21]     [,22]     [,23]     [,24]
## [1,] -22.391009 -14.8253180 -24.231501 -4.354447 -9.022073 -2.733719
## [2,]  -4.880423   0.9886361   4.352075  1.574892  7.256011  8.051135
##          [,25]     [,26]     [,27]     [,28]     [,29]      [,30]
## [1,] -7.256704 -8.136690 -8.491944 -2.219974 -0.702009 -0.6512773
## [2,]  5.058692  8.364195  7.001252  9.441244 10.913747  7.7929253
##            [,31]      [,32]     [,33]     [,34]     [,35]       [,36]
## [1,] -15.3936754 -12.142102 -6.206262 -8.918068 -20.66049 -36.7926378
## [2,]   0.4592864   1.625441  9.072093  4.319382   9.05894   0.4894266
##          [,37]     [,38]     [,39]     [,40]     [,41]     [,42]     [,43]
## [1,] -92.94742 -3.193814 -5.786602 -3.879734 -5.099848 -7.483961 -9.546316
## [2,]  37.18186 10.116372  1.106128  2.455868  2.476454  3.982401  5.487293
##            [,44]     [,45]     [,46]     [,47]      [,48]     [,49]
## [1,] -10.9406850 -5.662790 -4.553056 -4.684607 -6.0989805 -6.781810
## [2,]  -0.4221727  8.821038  7.843877  2.600304  0.5069754  5.218776
##          [,50]     [,51]     [,52]     [,53]    [,54]     [,55]     [,56]
## [1,] -4.326650 -5.746876 -7.777732 -3.352765 -3.50116 -8.549383 -4.206137
## [2,]  7.797064  3.934303  1.772232  6.153656  9.07543  7.587629  3.543481
##          [,57]      [,58]     [,59]     [,60]      [,61]    [,62]
## [1,] -11.21682 -11.228558 -3.760284 -5.203368 -20.911518 -8.01081
## [2,]  17.78247   7.123334  4.367238  6.810790   7.266368 11.54899
##          [,63]      [,64]     [,65]     [,66]     [,67]     [,68]
## [1,] -4.801189 -10.125376 -1.488239 -4.397326 -8.582007 -4.735429
## [2,] 10.340835   4.737047  4.988849  7.064662  3.041006  1.331775
##          [,69]     [,70]     [,71]     [,72]     [,73]     [,74]
## [1,] -4.889150 -2.765026 -7.539645 -4.621756 -6.251964 -9.030618
## [2,]  1.513738  3.711226  1.404931  8.860615  2.145722  4.861056
##           [,75]      [,76]     [,77]       [,78]     [,79]     [,80]
## [1,] -19.083643 -7.8099182 -8.252201 -15.5075389 -5.552006 -6.409345
## [2,]   3.691285 -0.4009293  7.095524  -0.3283292  3.499032  4.001604
##          [,81]     [,82]     [,83]      [,84]     [,85]      [,86]
## [1,] -4.574273 -8.024120 -5.168369 -13.860545 -6.683481 -0.7505132
## [2,]  5.455186  3.305847  8.285768   8.587947  4.779171  5.5517204
##          [,87]      [,88]     [,89]     [,90]     [,91]     [,92]
## [1,] -1.280698 -11.299628 -3.743237 -8.101830 -3.619115 -9.799569
## [2,]  8.200250   2.930987  9.013682 -2.602141 10.658489  4.957554
##           [,93]      [,94]     [,95]     [,96]     [,97]      [,98]
## [1,] -12.716743 -12.362522 -2.190750 -8.529031 -6.031179 -10.624314
## [2,]   3.961452   8.001405  3.236597  4.169281  3.040771   8.005183
##          [,99]     [,100]    [,101]    [,102]    [,103]    [,104]   [,105]
## [1,] -2.363686 -15.852478 -7.564902 -10.23729 -11.83354 -6.788493 -1.74401
## [2,]  2.886716  -2.323283  7.850104   2.00618 -11.83354  1.270906  7.29163
##           [,106]    [,107]    [,108]    [,109]    [,110]    [,111]
## [1,] -7.62197617 -5.588375 -4.527044 -5.317712 -6.665714 -9.238344
## [2,] -0.01782443 -1.017068  1.537542  4.621445  4.744112  5.860747
##         [,112]
## [1,] -5.324298
## [2,]  4.591907
## 
## $out
##  [1]  68.56582  55.03351 -56.19452  46.58311 -40.62345  34.58560  59.62393
##  [8] -41.43351 -83.77559  41.16684 -48.40559  48.68493 -41.72507  36.62259
## [15] -38.16378 -38.99232 -45.93325 -34.66257 -34.14713 -37.10013 -42.38954
## [22] -38.45764 -42.96275 -42.88271  52.92046 -38.53244 -41.69406 -39.73615
## [29] -44.92581 -60.63356 -41.81585  30.02607  35.23090  55.77892  59.51291
## [36]  17.03042  61.08003  47.57330 -33.82651 -34.58068  52.07212  52.54111
## [43]  58.19308  36.65890  31.95529  25.57077 -42.70070  33.00581  38.01076
## [50]  42.79479 -51.95851 -52.10624  39.96641  47.09931 -39.82029 -36.76943
## [57]  56.24366 -43.67188 -61.49849  33.88853  33.73730 -49.23307  51.81621
## [64]  39.05946  33.87278 -36.58964 -39.71718  38.69677 -39.27612  42.10430
## [71]  36.96515 -39.41693  40.33716  32.71640  44.04512 -49.46229  49.20671
## [78] -64.66606 -65.61321 -52.08234  24.83272  39.16692 -23.59222 -39.26652
## [85]  41.02335 -38.19916 -51.99974  40.34830  48.82099  49.69036  61.80130
## [92] -46.56122  46.28912 -55.11982  68.60874  55.39710  55.39710
## 
## $group
##  [1]   1   1   1   3   3   3   4   9   9   9   9   9  13  13  13  13  13
## [18]  13  13  13  14  14  14  14  14  14  14  14  14  14  14  16  16  19
## [35]  19  22  23  24  30  30  31  31  32  39  39  39  40  42  44  44  45
## [52]  45  48  48  51  53  56  56  58  59  59  65  68  68  68  68  70  70
## [69]  70  71  73  73  76  76  78  82  82  86  86  86  90  90  90  95  95
## [86]  95  95  95  96  97 100 105 105 105 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##            [,1]       [,2]       [,3]      [,4]       [,5]       [,6]
## [1,] -38.200747 -40.825999 -17.252445 -32.71051 -15.351290 -21.901330
## [2,] -14.809169 -14.847053  -1.649083 -16.33815  -8.925348 -10.046860
## [3,]  -1.956534  -1.723715   1.845443  -1.92851   1.609199  -3.588798
## [4,]  14.147836  17.286920  13.809389  11.74525   8.420746  13.590412
## [5,]  50.595944  36.063514  34.068108  33.14205  22.565479  40.393324
##            [,7]       [,8]       [,9]       [,10]      [,11]      [,12]
## [1,] -3.6793379 -17.047447 -37.096920 -35.5462037 -42.811173 -42.757974
## [2,] -3.6476791  -6.747462  -9.935577 -12.9900835 -19.496160 -10.133308
## [3,] -0.2711341  -3.991306   1.200953   0.8843957   5.723933   2.199528
## [4,]  2.0598182   2.065863  10.577938   9.6053479  12.415481  12.346279
## [5,]  2.0990384   6.457753  40.043504  32.6611084  28.713572  43.436207
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -33.023483 -39.687658 -34.091205 -28.026140 -29.580250 -24.039494
## [2,]  -9.523083  -8.519944  -7.864722  -9.019501 -13.450484 -13.423780
## [3,]   1.234577   3.275365   5.468611   2.275028  -7.207990   3.735621
## [4,]  10.028829  13.084275  14.148058   6.434810   8.894622  17.884188
## [5,]  37.492311  43.886829  28.801944  16.791713  39.319552  27.561779
##          [,19]      [,20]      [,21]       [,22]      [,23]      [,24]
## [1,] -45.14404 -30.616428 -21.223432 -10.2477518 -44.228704 -35.242559
## [2,] -14.36410 -23.547830 -15.533473  -2.2285816 -14.684863 -11.917213
## [3,]  -7.02917  -2.212079  -9.666812   0.4444751  -2.908150   2.642603
## [4,]  10.75813  16.161660   1.265963   6.2364620   8.276752  12.596231
## [5,]  26.59607  34.611597  20.113096  12.0315224  30.425183  46.974338
##            [,25]     [,26]      [,27]     [,28]      [,29]      [,30]
## [1,] -26.0830596 -39.11273 -50.897545 -28.49711 -19.927936 -25.179871
## [2,] -14.2903716 -14.27100 -16.539652 -10.83205  -9.178449  -8.876196
## [3,]  -0.3384946  -4.67894   5.469956   3.43054   4.419424   2.667278
## [4,]  16.2150056  11.83526  14.522139  13.03367  11.371781   7.145780
## [5,]  36.8809911  47.62920  43.315785  36.39284  18.046435  27.096734
##           [,31]      [,32]      [,33]      [,34]      [,35]     [,36]
## [1,] -34.632203 -39.815992 -35.504912 -42.926441 -26.405857 -26.50477
## [2,] -18.402023 -12.816882 -10.322867 -16.753667 -16.048595 -26.50477
## [3,]  -6.872592  -4.022802   1.456333  -3.200241  -5.620167 -18.52198
## [4,]  12.905255   9.763758  13.188737  16.579666  13.657409 -10.53919
## [5,]  53.375685  37.115348  43.849013  46.139482  19.477521 -10.53919
##            [,37]      [,38]      [,39]      [,40]      [,41]      [,42]
## [1,] -35.2947698 -22.323961 -28.815425 -40.620187 -46.393030 -20.762800
## [2,] -35.2947698  -5.963661  -8.578004 -12.539337 -14.054221 -14.820491
## [3,] -18.0948127   1.698506   1.201904   0.378215  -2.258456   1.041926
## [4,]  -0.8948557  13.430780   8.443345  10.330281  12.988341   9.120110
## [5,]  -0.8948557  29.594557  29.425191  33.153330  47.955069  23.465408
##           [,43]      [,44]      [,45]      [,46]      [,47]      [,48]
## [1,] -30.382156 -47.648323 -56.490587 -43.073604 -26.364882 -32.557494
## [2,] -13.916568 -16.056642 -16.428936 -11.159455 -10.439283 -10.384686
## [3,]  -1.249261  -5.074353   1.575308   3.958578  -1.344375  -1.554807
## [4,]   6.275462  16.014745  10.634100  10.405153   9.658410   7.618509
## [5,]  25.514402  43.058044  42.175777  24.618855  24.372347  34.011492
##            [,49]     [,50]      [,51]      [,52]      [,53]      [,54]
## [1,] -35.3437898 -39.96864 -37.959150 -25.470307 -36.088484 -25.971787
## [2,] -12.1025977 -10.14409  -8.985862 -10.360260  -8.925703 -14.443542
## [3,]  -0.3982405   1.82376   1.383990  -5.250214   1.376679   2.195808
## [4,]  10.7825636  13.14460  10.435727  14.765761  14.972298  12.456373
## [5,]  29.0614874  30.91173  29.009345  31.116008  30.415739  34.997328
##           [,55]       [,56]      [,57]      [,58]      [,59]       [,60]
## [1,] -22.903583 -31.1392313 -12.162365 -58.506273 -37.910186 -34.9417331
## [2,] -14.460204  -6.0022203  -5.746893 -14.292571 -13.114601 -10.6220984
## [3,]  -3.242117   0.6644463   1.112919  -3.943067   1.885399  -0.4462605
## [4,]  10.462103  11.1046307  10.320881  20.817475  15.564846  11.2329364
## [5,]  44.203268  32.5901194  19.117587  53.254464  34.888594  32.0163853
##          [,61]       [,62]      [,63]      [,64]      [,65]       [,66]
## [1,] -8.225215 -21.2604650 -20.996731 -39.039331 -39.349246 -33.3319282
## [2,] -5.805848  -9.5654846 -11.877100 -15.607570 -10.619308  -9.9906073
## [3,] -3.386482   0.3687193   1.855621  -6.465905   1.743119   0.2758621
## [4,] 18.357478   5.8097951  11.456233  15.436801  13.640593  12.6533197
## [5,] 40.101437  13.2175875  30.779982  45.533508  37.107285  38.9686327
##            [,67]      [,68]      [,69]        [,70]      [,71]       [,72]
## [1,] -31.1270084 -33.092309 -33.384904 -33.96257493 -38.287104 -27.0530629
## [2,] -10.6814311 -10.595196 -10.198654  -9.14279026 -11.185895 -12.5414549
## [3,]  -0.8922643  -1.758091  -1.087010  -0.03913421  -2.962607   0.4537526
## [4,]  14.8143991   6.905381   9.801346   7.96406081   7.976864   9.2432455
## [5,]  41.4770295  29.729576  29.817321  30.03934034  29.735792  16.9163022
##           [,73]      [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -26.409375 -40.351945 -39.002481 -28.725940 -39.831311 -21.220032
## [2,] -12.562882 -14.750868 -19.789970 -10.593629 -14.948530 -16.888006
## [3,]  -1.907140  -4.670668  -8.739155  -3.250711  -3.414067  -6.219212
## [4,]   6.140284  15.022646  11.372462   6.073038   6.585933  10.992143
## [5,]  27.020653  32.521028  26.205684  29.452801  37.262185  41.922087
##            [,79]      [,80]       [,81]      [,82]       [,83]       [,84]
## [1,] -27.5001861 -40.248253 -35.6822831 -32.960916 -35.7139899 -32.0035789
## [2,] -11.5407297 -14.345029 -10.0660654 -13.751617 -14.3774616 -15.1008449
## [3,]  -0.7209196  -3.527523  -0.6242835  -4.631985  -0.9260946  -0.6510754
## [4,]  12.5023676  13.050117  13.3853016  12.034681  14.9462867  21.8095956
## [5,]  42.5023676  52.119803  30.5188387  46.468476  41.8330455  39.3329495
##            [,85]      [,86]      [,87]       [,88]      [,89]      [,90]
## [1,] -31.5529718 -43.124423 -41.471878 -48.3896017 -36.849519 -16.686707
## [2,] -12.9995462 -12.192121 -11.130538 -16.7355674 -13.846324  -7.575063
## [3,]   0.1891464   3.098732   1.572911   0.8472997   2.982645  -5.130087
## [4,]  12.8888101  15.414735  10.632598  13.1736353  12.482217   2.934024
## [5,]  42.0004538  48.932626  36.913332  45.9653333  45.520079  17.542970
##             [,91]      [,92]      [,93]      [,94]         [,95]     [,96]
## [1,] -29.73106436 -28.518744 -35.630146 -35.410032 -39.060433294 -52.15202
## [2,] -12.84430549 -15.405503 -19.117387 -15.172985 -11.124878894 -12.68733
## [3,]  -0.09758755  -4.068974  -2.482302  -2.030087  -0.008442737  -1.36661
## [4,]  12.48583282  12.157508  15.714370   2.762913   8.753486894  18.42911
## [5,]  33.38217677  51.041072  37.483095   7.363508  31.320098055  48.83769
##            [,97]        [,98]      [,99]     [,100]     [,101]     [,102]
## [1,] -48.9204869 -26.50146209 -33.241607 -45.082185 -32.387670 -20.260946
## [2,] -12.6378260  -7.71209048  -8.950244 -16.186962 -15.279369 -12.037659
## [3,]  -0.6250459  -0.06288295  -0.128366  -6.107562  -1.993649  -1.377382
## [4,]  12.5597941   7.36386661   8.315013  12.696602  10.385678   6.201604
## [5,]  46.0256456  21.01333761  31.823008  52.476510  46.063946  20.798564
##         [,103]     [,104]    [,105]     [,106]     [,107]      [,108]
## [1,] -11.32921 -35.621536 -47.05328 -33.256831 -39.499912 -52.7481912
## [2,] -11.32921 -11.251357 -11.49142 -13.124504 -13.263206 -19.2859929
## [3,] -11.32921  -3.028070   1.40772  -4.182216  -2.775317   0.4248809
## [4,] -11.32921   9.865079  12.57714   9.005195  12.485790  16.0311167
## [5,] -11.32921  36.311653  46.33330  38.370976  47.394309  62.1081409
##         [,109]     [,110]     [,111]      [,112]
## [1,] -50.10609 -36.943086 -29.649814 -39.8169778
## [2,] -18.52197 -13.383874 -13.060187 -18.2724801
## [3,]  -2.38332   1.044134  -1.428302   0.6822288
## [4,]  19.10689  15.222899  10.645431  16.6452390
## [5,]  52.57331  28.150985  32.192798  50.9597650
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]       [,4]      [,5]       [,6]
## [1,] -4.973336 -7.263359 -4.260654 -11.850338 -5.992087 -12.925520
## [2,]  1.060268  3.815929  7.951539   7.993318  9.210486   5.747925
##           [,7]       [,8]      [,9]     [,10]     [,11]     [,12]
## [1,] -2.772234 -10.218780 -1.454289 -4.114714 -6.881165 -2.107632
## [2,]  2.229966   2.236167  3.856194  5.883505 18.329030  6.506687
##           [,13]    [,14]      [,15]     [,16]       [,17]     [,18]
## [1,] -0.7153098 1.107812 -0.9899087 -2.110540 -13.6538200 -6.361704
## [2,]  3.1844644 5.442919 11.9271305  6.660596  -0.7621591 13.832946
##          [,19]     [,20]       [,21]     [,22]     [,23]     [,24]
## [1,] -17.27789 -15.58850 -19.6991636 -3.008875 -9.223575 -3.196352
## [2,]   3.21955  11.16435   0.3655404  3.897826  3.407274  8.481557
##          [,25]      [,26]     [,27]     [,28]     [,29]     [,30]
## [1,] -8.157323 -13.473021 -2.194681 -2.254130 -2.350906 -1.551843
## [2,]  7.480334   4.115142 13.134594  9.115209 11.189755  6.886398
##           [,31]      [,32]     [,33]      [,34]      [,35]       [,36]
## [1,] -15.756858 -10.888923 -5.325994 -11.858588 -19.771749 -36.3591781
## [2,]   2.011673   2.843319  8.238660   5.458105   8.531415  -0.6847766
##          [,37]     [,38]     [,39]     [,40]     [,41]     [,42]     [,43]
## [1,] -56.52738 -4.556514 -4.666790 -3.473682 -6.126799 -7.635988 -7.074000
## [2,]  20.33776  7.953526  7.070597  4.230112  1.609887  9.719840  4.575478
##           [,44]     [,45]     [,46]     [,47]     [,48]     [,49]
## [1,] -13.895353 -5.983592 -3.145942 -5.747912 -4.715368 -6.115410
## [2,]   3.746648  9.134208 11.063098  3.059162  1.605753  5.318929
##          [,50]     [,51]      [,52]     [,53]     [,54]     [,55]
## [1,] -4.486724 -3.529721 -11.690261 -4.669575 -4.441869 -9.995259
## [2,]  8.134244  6.297700   1.189832  7.422933  8.833484  3.511024
##          [,56]     [,57]      [,58]     [,59]     [,60]     [,61]
## [1,] -2.415774 -8.482496 -15.037842 -4.795715 -5.430375 -25.42859
## [2,]  3.744667 10.708335   7.151708  8.566514  4.537854  18.65563
##          [,62]     [,63]      [,64]     [,65]     [,66]     [,67]
## [1,] -6.644049 -6.833935 -15.004433 -1.397049 -6.048749 -7.108130
## [2,]  7.381488 10.545178   2.072622  4.883288  6.600473  5.323601
##          [,68]     [,69]     [,70]     [,71]     [,72]     [,73]
## [1,] -4.909205 -4.417943 -3.042337 -6.972915 -6.723277 -6.312343
## [2,]  1.393023  2.243922  2.964069  1.047700  7.630783  2.498064
##           [,74]      [,75]      [,76]     [,77]      [,78]     [,79]
## [1,] -13.259351 -20.680795 -7.0918222 -9.165249 -15.831849 -5.022224
## [2,]   3.918014   3.202484  0.5904007  2.337115   3.393425  3.580385
##          [,80]     [,81]      [,82]     [,83]     [,84]     [,85]
## [1,] -8.776519 -6.029042 -10.705496 -8.075204 -12.31477 -5.714780
## [2,]  1.721473  4.780475   1.441526  6.223015  11.01262  6.093073
##           [,86]     [,87]     [,88]     [,89]     [,90]     [,91]
## [1,] -0.2466792 -3.442767 -6.624615 -3.678541 -7.936735 -6.589951
## [2,]  6.4441442  6.588589  8.319215  9.643831 -2.323438  6.394776
##           [,92]     [,93]      [,94]     [,95]     [,96]     [,97]
## [1,] -11.228477 -13.07363 -12.049337 -2.875571 -9.924958 -7.000106
## [2,]   3.090529   8.10903   7.989163  2.858686  7.191739  5.750015
##          [,98]     [,99]      [,100]    [,101]    [,102]    [,103]
## [1,] -6.939128 -2.588040 -12.8362275 -9.162081 -6.115038 -11.32921
## [2,]  6.813362  2.331308   0.6211038  5.174783  3.360274 -11.32921
##          [,104]    [,105]     [,106]     [,107]    [,108]    [,109]
## [1,] -7.0446201 -3.461308 -7.5783093 -5.0355071 -2.770277 -8.026405
## [2,]  0.9884796  6.276749 -0.7861219 -0.5151278  3.620039  3.259765
##         [,110]    [,111]    [,112]
## [1,] -6.386487 -9.413708 -4.626507
## [2,]  8.474755  6.557105  5.990964
## 
## $out
##  [1]  63.18393  64.25334  59.70159 -32.88438  46.10803 -16.69747  11.47638
##  [8] -23.51526  14.08075 -44.75445 -83.88504 -41.38391 -48.71986  44.88822
## [15]  54.96714 -46.66107 -39.76218 -39.96873  43.12632 -40.04562 -46.32131
## [22] -46.01716 -45.21200 -39.12908 -45.26640 -48.19719 -41.05542 -42.63958
## [29] -44.75688 -54.08954 -41.20475 -44.41629  37.16498  44.42797  65.18773
## [36] -20.77435  58.86856  56.42832 -38.50360  37.57308  46.72561  60.80057
## [43] -48.22192  34.31514  45.94854  47.45823  43.60241  35.62170  51.67934
## [50] -52.12424  38.00998 -40.04941 -34.93626  37.95524 -56.22971  52.15126
## [57]  29.02662 -53.81783 -60.74926  54.46839  47.11595  36.21487  43.46274
## [64] -45.72854 -42.07211  61.59544  55.62109 -38.80857  33.58558  33.85614
## [71]  54.20243  48.14255 -48.91057 -69.32060  56.18151 -70.97646 -67.00959
## [78] -63.52672  25.98635  19.42135  38.64343 -28.12053  39.98357  49.86235
## [85] -45.23715  35.32115  47.28247 -44.88759  61.13998  46.71125 -54.44717
## [92] -52.97309  52.19699  54.03378 -53.66238  56.83314  56.83314
## 
## $group
##  [1]   1   1   1   3   5   7   7   7   7   9   9   9   9   9   9  12  13
## [18]  13  13  13  13  13  13  13  14  14  14  14  14  14  14  14  16  17
## [35]  19  22  23  23  30  30  32  32  38  39  40  40  43  48  48  48  56
## [52]  56  56  57  59  60  62  65  67  68  68  68  69  70  70  71  71  76
## [69]  76  76  77  77  79  79  80  83  86  86  90  90  90  90  95  95  98
## [86]  99  99  99 100 102 105 105 107 107 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##            [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
## [1,] -44.849988 -34.798855 -47.703259 -32.554513 -19.027244 -25.641310
## [2,] -13.255096 -15.656388 -10.252491 -14.896232  -5.093029 -11.611999
## [3,]  -2.408684  -1.100659  -1.366682   1.332328   4.560226  -5.726601
## [4,]  15.836016  15.407849  24.144682   9.508757  13.315249  17.782971
## [5,]  45.591986  42.395475  43.062772  33.734755  35.017931  41.684859
##            [,7]      [,8]       [,9]       [,10]       [,11]      [,12]
## [1,] -23.214562 -4.870598 -42.922463 -43.2337145 -24.4508751 -35.541907
## [2,] -14.474947 -4.870598 -10.731939 -12.0641530 -16.4778036  -9.171177
## [3,]   3.060509 -2.185159   1.304703  -0.7319663  -0.4645649   3.546549
## [4,]   4.035466  7.242261  12.170627  13.4543463   9.4418306  11.855875
## [5,]  31.089368 14.069324  44.765789  50.3548537  19.6576241  37.489769
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -31.969708 -41.187450 -30.941534 -30.501026 -29.640261 -42.580245
## [2,]  -8.246582  -8.620063 -12.785280  -9.661116 -14.932527 -12.073471
## [3,]   1.122888   4.345022   4.117934   2.134906  -3.447606   5.207823
## [4,]   9.814678  13.218271  12.464749  11.965535   5.487483  12.254686
## [5,]  36.704605  37.826731  31.683481  33.978652  30.923393  27.932648
##           [,19]       [,20]      [,21]      [,22]       [,23]      [,24]
## [1,] -20.822963 -36.6929240 -23.037904 -13.920312 -45.3341088 -27.544828
## [2,] -14.820963 -15.2667964 -15.803523  -6.024332 -15.0528993  -6.885790
## [3,]  -8.317687  -0.8917635  -1.902476   2.528081  -0.2299789   5.039119
## [4,]  -2.597323  11.1761144   1.099873   6.356178  15.3742098  11.896792
## [5,]   7.249178  30.6374777  18.097524  16.109492  58.6345933  31.319744
##           [,25]      [,26]      [,27]      [,28]      [,29]      [,30]
## [1,] -32.337857 -18.066944 -40.801828 -25.345715 -18.792487 -27.360184
## [2,]  -6.864364  -7.712785 -17.249631  -8.501969  -6.026868  -9.373863
## [3,]   1.320820  -2.496813   3.262419   3.447462   4.540846   2.543100
## [4,]  11.489158   6.047658  15.210508  12.315774  11.811701   8.464973
## [5,]  37.196107  25.096310  47.889214  39.477205  24.999992  25.980563
##           [,31]      [,32]       [,33]      [,34]      [,35]     [,36]
## [1,] -47.128885 -41.733735 -36.0104067 -28.875110 -24.531309 -20.95185
## [2,] -18.076807 -10.160856  -8.2812621 -18.517693 -16.346297 -20.95185
## [3,]  -4.370616  -3.478285   0.4791803  -2.913504  -5.753345 -15.61902
## [4,]  12.680991  11.286692  14.6979119  16.024018  12.437967 -10.28619
## [5,]  50.496415  36.269069  41.3645786  48.221924  15.951465 -10.28619
##          [,37]     [,38]       [,39]        [,40]       [,41]      [,42]
## [1,] -9.121349 -28.54546 -27.9591966 -29.73308364 -33.7488148 -24.571976
## [2,] -9.121349 -10.96972  -9.5865927  -9.25674311 -11.2426170 -15.453718
## [3,] -1.830375   1.08244   0.4803111  -0.06394054   0.5287722  -9.238129
## [4,]  5.460598  17.11339   8.2884515  10.29339853  12.8719847   4.206869
## [5,]  5.460598  31.39545  32.6842628  38.57445867  38.2632012   9.128084
##           [,43]      [,44]        [,45]      [,46]      [,47]      [,48]
## [1,] -32.156575 -31.040261 -52.59005115 -28.490762 -26.640490 -28.316357
## [2,] -12.464080 -15.609836 -13.84521713 -12.185134 -11.452079  -8.403191
## [3,]  -5.136011  -3.689996  -0.08477469   1.067799  -2.511278  -1.736524
## [4,]   7.972908  11.583924  14.17043180   8.957173  10.817372   7.200998
## [5,]  37.172566  47.785497  36.80250316  23.707144  30.026855  29.866340
##            [,49]      [,50]       [,51]      [,52]      [,53]      [,54]
## [1,] -27.6380076 -41.659716 -31.2935401 -23.431234 -49.316577 -32.223498
## [2,] -11.5025799 -11.940925  -8.5279206  -8.608314 -10.585201 -10.592749
## [3,]  -0.9348661   2.168382  -0.5852215  -3.608314   2.818352  -1.608977
## [4,]   8.6660398  15.227648   8.4929054  10.256259  16.312729  11.814142
## [5,]  35.0514502  34.277008  28.2793508  30.610418  46.771572  44.813859
##           [,55]      [,56]      [,57]      [,58]      [,59]       [,60]
## [1,] -26.712339 -29.967435 -10.686303 -58.057320 -39.768766 -34.9657755
## [2,] -13.509890  -6.127730  -3.292852 -15.957417 -11.790305 -14.4152193
## [3,]  -1.096506   1.601414  13.513823  -3.734458  -1.009036  -0.2382903
## [4,]  11.516864  10.893096  15.775595  14.775955   8.760620  12.5423147
## [5,]  40.662646  19.954238  16.895247  54.775955  36.158986  32.1881554
##          [,61]       [,62]       [,63]      [,64]      [,65]      [,66]
## [1,] -2.514170 -20.2419191 -19.9815360 -36.040158 -40.951325 -33.387919
## [2,]  1.410202 -12.1552462 -11.3981732 -18.635068 -11.453932  -7.206881
## [3,]  5.334573  -0.6526317  -0.1084481  -5.875116  -0.924603   4.399223
## [4,] 30.591763  14.9714507   7.4663990  16.958992  13.334397  12.788838
## [5,] 55.848953  24.2731986  31.8622103  49.773623  42.802545  20.974725
##          [,67]      [,68]      [,69]      [,70]      [,71]      [,72]
## [1,] -30.33615 -37.361619 -23.217229 -29.274673 -24.039667 -22.102701
## [2,] -11.01599 -11.640259  -8.979578  -7.572826  -8.228354 -10.472509
## [3,]  -1.68320  -3.129846  -1.539248   0.510478  -2.721977   1.584792
## [4,]  14.48423   5.734726   9.138565   9.093841   7.985284   7.771797
## [5,]  44.16608  29.804857  31.778127  27.958413  26.754683  12.439010
##           [,73]      [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -26.065143 -37.354073 -30.214981 -23.779887 -29.736486 -33.739431
## [2,]  -9.227504 -17.254033 -19.787872 -12.170232  -8.680727 -24.239204
## [3,]  -2.206678  -1.735567   1.420506  -2.524391   2.280644  -5.435001
## [4,]   5.241257  11.880785  11.347556   6.217629   8.146762  16.224488
## [5,]  26.657894  49.112396  26.281480  28.145329  30.796699  39.557822
##           [,79]      [,80]       [,81]      [,82]      [,83]      [,84]
## [1,] -29.299932 -49.883516 -32.8746176 -37.849075 -57.253059 -32.815603
## [2,]  -8.657169 -13.429823  -9.7548388 -13.100348 -12.822297 -10.845908
## [3,]  -2.130881  -2.119103   0.6851528  -4.413162   1.469384   2.080267
## [4,]   6.148861  12.669130  13.1878603  11.263977  20.865195  18.746934
## [5,]  25.440542  45.949825  24.1462084  45.940998  45.969325  38.319824
##           [,85]      [,86]       [,87]       [,88]      [,89]      [,90]
## [1,] -32.562135 -49.640776 -38.9345363 -40.0540240 -35.815476 -19.934101
## [2,] -12.385055 -12.947007 -10.6488199 -14.9948337 -12.208950  -9.757021
## [3,]  -0.858993   2.939424   0.9430151  -0.2982915   1.676076  -4.757021
## [4,]  12.047231  14.533072  10.9248589  11.7083079  13.049096   8.700111
## [5,]  41.125358  43.055229  40.3758772  48.1978949  36.028271  35.597138
##           [,91]      [,92]      [,93]      [,94]       [,95]      [,96]
## [1,] -32.814726 -29.413113 -39.249385 -27.589460 -37.2000124 -28.022917
## [2,] -12.926798 -12.319337 -22.072736 -15.452473 -10.6062955 -13.550204
## [3,]  -1.564101  -3.454765  -2.365012  -2.779869   0.7035079  -1.810357
## [4,]  13.708502  10.659837  14.735781   6.061444   8.4041764  16.800003
## [5,]  36.542833  26.264026  37.989148  24.085301  32.7999876  47.154162
##           [,97]      [,98]      [,99]     [,100]      [,101]      [,102]
## [1,] -38.626500 -26.711440 -33.179811 -40.569655 -48.7949005 -20.1275938
## [2,] -13.449420  -7.227412  -9.189575 -16.769635 -11.8363213 -11.9713401
## [3,]  -3.045147  -2.083931   1.242031  -8.906239  -0.4298052  -0.3918898
## [4,]   9.102645  10.204446   8.914084   7.389010  13.0915833   6.8701882
## [5,]  21.913854  23.827375  33.123850  36.607741  47.8307764  23.5187103
##        [,103]     [,104]     [,105]     [,106]     [,107]     [,108]
## [1,] -13.2529 -33.977835 -43.887805 -34.854228 -44.650230 -63.326999
## [2,] -13.2529 -11.276479 -11.976970 -14.122877 -13.603620 -17.440170
## [3,] -13.2529  -3.120226   2.496804  -4.525128  -3.214685  -1.319338
## [4,] -13.2529  12.056854  11.688741  10.963532  12.717514  16.561259
## [5,] -13.2529  35.817296  44.340550  36.987516  48.949447  52.546537
##           [,109]      [,110]    [,111]      [,112]
## [1,] -54.0121321 -37.9659473 -29.43030 -49.6252128
## [2,] -15.8073871 -10.3097525 -11.76976 -16.5973053
## [3,]  -0.6714748   0.8256752  -1.78636  -0.5528152
## [4,]  16.3883748  14.8673271  11.63449  13.5301856
## [5,]  49.6513777  32.0670281  34.36057  52.5794866
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##            [,1]      [,2]      [,3]      [,4]      [,5]       [,6]
## [1,] -5.4394573 -6.455889 -14.95357 -7.289924 -3.506524 -17.337614
## [2,]  0.6220898  4.254571  12.22020  9.954579 12.626976   5.884412
##           [,7]       [,8]      [,9]     [,10]      [,11]      [,12]
## [1,] -5.050998 -10.744075 -1.659774 -6.377788 -10.702820 -0.4823006
## [2,] 11.172015   6.373757  4.269180  4.913855   9.773691  7.5753979
##           [,13]    [,14]    [,15]     [,16]     [,17]    [,18]      [,19]
## [1,] -0.6783383 2.153980 -3.29039 -4.002227 -9.338109 -2.63840 -13.304373
## [2,]  2.9241142 6.536064 11.52626  8.272038  2.442898 13.05405  -3.331001
##          [,20]     [,21]     [,22]     [,23]     [,24]     [,25]     [,26]
## [1,] -9.799247 -11.99691 -2.522601 -8.598733 0.5652216 -3.383368 -7.132116
## [2,]  8.015720   8.19196  7.578763  8.138775 9.5130168  6.025009  2.138490
##          [,27]     [,28]    [,29]     [,30]     [,31]     [,32]     [,33]
## [1,] -4.747268 -1.511199 -1.33612 -2.154460 -13.09895 -9.999866 -6.149558
## [2,] 11.272106  8.406122 10.41781  7.240661   4.35772  3.043295  7.107919
##           [,34]      [,35]      [,36]     [,37]     [,38]     [,39]
## [1,] -11.885727 -19.465821 -27.535004 -18.12175 -7.974817 -5.682723
## [2,]   6.058719   7.959131  -3.703034  14.46099 10.139696  6.643345
##          [,40]     [,41]      [,42]       [,43]      [,44]     [,45]
## [1,] -3.356743 -2.920737 -16.364637 -11.0314118 -11.169441 -7.909747
## [2,]  3.228862  3.978281  -2.111621   0.7593903   3.789449  7.740197
##          [,46]     [,47]     [,48]     [,49]     [,50]     [,51]     [,52]
## [1,] -5.897593 -7.390662 -4.475926 -5.973389 -5.193426 -4.891533 -8.443489
## [2,]  8.033191  2.368105  1.002878  4.103657  9.530189  3.721090  1.226862
##          [,53]     [,54]     [,55]     [,56]     [,57]      [,58]
## [1,] -3.986891 -7.137979 -7.877950 -1.463317  2.126454 -13.446204
## [2,]  9.623595  3.920025  5.684937  4.666145 24.901192   5.977287
##          [,59]     [,60]     [,61]     [,62]     [,63]     [,64]     [,65]
## [1,] -5.796545 -6.386046 -21.28524 -13.02530 -7.133795 -15.66500 -4.133171
## [2,]  3.778472  5.909466  31.95438  11.72004  6.916899   3.91477  2.283965
##          [,66]     [,67]        [,68]     [,69]     [,70]      [,71]
## [1,] -1.185725 -7.900135 -6.258346464 -4.556763 -2.415448 -6.1151038
## [2,]  9.984171  4.533734 -0.001346319  1.478268  3.436404  0.6711505
##          [,72]     [,73]      [,74]     [,75]     [,76]     [,77]
## [1,] -4.425845 -5.614542 -10.140007 -10.51079 -6.762180 -2.213453
## [2,]  7.595429  1.201185   6.668873  13.35180  1.713398  6.774740
##           [,78]      [,79]     [,80]     [,81]      [,82]     [,83]
## [1,] -19.386245 -4.7796762 -7.119745 -4.602375 -10.151752 -6.743603
## [2,]   8.516244  0.5179145  2.881538  5.972681   1.325428  9.682371
##          [,84]     [,85]      [,86]     [,87]     [,88]     [,89]
## [1,] -7.271071 -6.430858 -0.3906252 -4.028999 -6.969269 -4.714272
## [2,] 11.431605  4.712872  6.2694725  5.915029  6.372686  8.066424
##           [,90]     [,91]     [,92]      [,93]      [,94]     [,95]
## [1,] -9.6863442 -8.390991 -9.423614 -13.557420 -14.797852 -2.038441
## [2,]  0.1723021  5.262789  2.514085   8.827396   9.238114  3.445457
##           [,96]     [,97]      [,98]     [,99]     [,100]    [,101]
## [1,] -10.157960 -8.750876 -10.034718 -1.337085 -14.534195 -7.392348
## [2,]   6.537246  2.660581   5.866857  3.821147  -3.278282  6.532737
##         [,102]   [,103]    [,104]    [,105]     [,106]     [,107]
## [1,] -5.285984 -13.2529 -7.558451 -2.290729 -8.3749674 -5.5250960
## [2,]  4.502205 -13.2529  1.317999  7.284336 -0.6752877 -0.9042747
##         [,108]    [,109]    [,110]    [,111]    [,112]
## [1,] -4.395466 -5.499775 -5.714081 -9.670251 -5.133265
## [2,]  1.756789  4.156826  7.365432  6.097530  4.027635
## 
## $out
##   [1]  64.00793 -42.47062 -47.58270 -83.20830 -57.38913  46.89136 -47.18357
##   [8] -41.97081 -45.08824 -42.19188 -37.33107 -37.44843 -41.86327 -51.27219
##  [15] -45.65804 -44.75183  65.92098 -38.57039 -64.92163  52.64001 -51.85716
##  [22] -41.76786 -43.00101 -52.01892 -42.27318 -49.25894 -42.79269 -60.40236
##  [29] -64.94312 -44.11693 -52.71560  40.61089  45.59095  64.23438 -59.94599
##  [36]  46.49682 -36.05524 -45.17952  48.03461 -41.48721  30.90937 -39.12942
##  [43]  60.12544 -41.63380 -39.67013 -48.26720 -73.19724  55.52238  34.11175
##  [50]  37.58725  41.98821  41.68894 -47.78212  35.34569 -32.70491  35.00309
##  [57]  51.93558 -48.69316 -38.46182  54.57460 -45.55609  52.82850  47.27648
##  [64] -38.44264  50.43579  54.48939 -48.97239  41.78457 -46.30084 -37.79071
##  [71] -36.15054 -44.79063  28.98739  37.01902  34.92354 -52.10307  34.38552
##  [78]  45.80216  35.56411  39.90930  30.69057  48.53996  40.20478 -31.77744
##  [85] -38.23757 -54.88023  48.13890  50.04414 -63.47083 -62.33530 -52.02784
##  [92]  37.74708  49.87857  41.66456  50.62404 -43.77437  44.27972  49.68294
##  [99] -56.28527 -43.34826  40.35806  41.71304 -41.33067  54.40984  45.11815
## [106]  47.84236  47.38901  52.21193  63.27441 -62.61695 -49.04688 -55.15939
## [113]  83.37622  54.82446  66.30318  55.31396  55.31396
## 
## $group
##   [1]   1   7   8   9   9   9  12  13  13  13  13  13  13  13  13  13  13
##  [18]  13  14  14  14  14  14  14  14  14  14  14  14  14  14  17  17  19
##  [35]  22  24  24  24  25  25  26  26  32  40  40  41  41  41  42  42  42
##  [52]  43  47  48  48  48  54  56  56  57  59  68  68  68  68  68  69  70
##  [69]  70  70  70  71  73  73  76  77  77  77  79  79  79  79  79  79  79
##  [86]  82  82  82  86  86  87  90  92  95  95  95  97  97  98  99  99  99
## [103]  99 100 100 100 100 100 100 101 102 105 107 107 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##            [,1]        [,2]       [,3]       [,4]       [,5]       [,6]
## [1,] -46.605373 -40.7121146 -26.006805 -24.820949 -18.854659 -21.373253
## [2,] -14.468281 -13.3059384  -9.785766 -15.082308  -3.825141 -10.745241
## [3,]  -1.598594  -0.1949637   3.447884  -6.522804   2.845541  -2.169764
## [4,]  14.317177  12.3472325  14.584322   9.574823  11.051193  16.176519
## [5,]  56.731220  35.6891845  40.507052  43.344495  26.884772  40.587394
##             [,7]       [,8]       [,9]       [,10]      [,11]       [,12]
## [1,] -26.6001708 -21.950310 -43.014535 -38.2008179 -26.260244 -36.6641962
## [2,] -16.6133910  -6.067722 -10.974048 -10.8720075 -16.086132  -8.9484758
## [3,]   0.8531892  -2.837440   1.666316   0.6647383  -5.511023   0.4055835
## [4,]  16.3524446   5.322219  12.306099   9.6406084  11.941955  13.5173571
## [5,]  21.2800052   6.862891  47.097756  32.5815796  22.138718  46.4720029
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -33.110580 -35.762581 -32.370119 -21.867621 -32.239436 -16.694907
## [2,]  -9.090444  -6.904195  -6.879375  -5.403157 -11.760114 -10.929785
## [3,]   1.382855   2.168545   4.032640   3.153188  -5.138016  -3.689069
## [4,]   9.840803  12.682705  14.102112   6.995023   7.174975   9.895751
## [5,]  36.676836  39.993396  29.787291  20.483229  30.543376  29.636126
##           [,19]      [,20]      [,21]      [,22]       [,23]      [,24]
## [1,] -21.614851 -29.777005 -21.668228 -7.2349544 -45.7283167 -43.071856
## [2,] -13.521326 -12.978448 -19.269716 -3.4186943 -14.1611948 -10.227030
## [3,]  -7.860798  -6.490597 -10.295542 -0.9571597   0.1536256   2.940699
## [4,]  -2.184863  11.658049   9.116496  2.3750102  12.8540665  13.333061
## [5,]   2.773322  32.432271  20.398485  7.0232847  31.1327859  35.587690
##            [,25]       [,26]       [,27]      [,28]      [,29]      [,30]
## [1,] -27.8866818 -39.0461996 -35.0270117 -27.339130 -16.922171 -21.507919
## [2,] -12.7292829  -6.6943533 -14.6250732  -9.791978  -6.628190  -7.924935
## [3,]  -0.2238874  -0.4313787   0.8669976   2.759252   5.235239  -1.040026
## [4,]  13.5099857  15.3681410  14.3805739  11.927616  11.332543   8.334329
## [5,]  39.9030503  34.0519144  43.7787741  37.366889  20.288109  27.742564
##           [,31]      [,32]      [,33]      [,34]      [,35]      [,36]
## [1,] -37.226328 -27.412739 -34.855132 -30.660038 -42.588197 -13.559145
## [2,] -16.532736 -16.502589 -10.149113 -15.209897 -20.369066 -13.559145
## [3,]  -6.988769  -4.345659   1.517554  -2.937262  -5.299406 -10.825140
## [4,]   7.439378  13.935560  11.571807  10.280846   6.286874  -8.091135
## [5,]  28.857419  36.542641  42.987458  46.947513  23.854147  -8.091135
##           [,37]      [,38]      [,39]      [,40]      [,41]      [,42]
## [1,] -31.546027 -23.050754 -21.880500 -29.960482 -42.482306 -18.617302
## [2,] -31.546027  -7.112726 -10.217263  -9.960482 -13.490434 -11.427344
## [3,] -16.685431   1.209893  -2.406746   0.310731  -1.124137  -3.723266
## [4,]  -1.824835  12.358694  11.594244   8.841175  14.399199   4.793947
## [5,]  -1.824835  32.798190  36.390559  32.784889  42.063135  18.939374
##           [,43]      [,44]      [,45]      [,46]      [,47]       [,48]
## [1,] -26.903794 -32.671557 -33.450307 -28.492006 -34.852010 -34.3220352
## [2,] -13.788812 -12.563034 -13.710043 -11.629740 -12.193580  -9.4201777
## [3,]  -2.429550  -5.018076   3.137973   3.468642  -3.275613  -0.5487816
## [4,]   3.499007   6.543211  14.600010   8.861004   8.294187   8.8136109
## [5,]  13.761861  34.101306  32.993935  24.056603  34.939290  30.3830595
##            [,49]       [,50]      [,51]      [,52]      [,53]     [,54]
## [1,] -26.1055931 -39.7300359 -29.854009 -25.172463 -33.129378 -40.41665
## [2,] -10.3208688 -11.2971578  -8.579704  -9.878482  -7.386786 -12.96902
## [3,]  -0.4568611   0.1715827  -1.297901  -4.584501   1.165847   5.42784
## [4,]  10.4638555  12.0361756  10.341591   9.711299  14.430464  16.71627
## [5,]  23.5032081  32.6218103  30.537190  29.336794  31.951886  33.83652
##            [,55]       [,56]      [,57]      [,58]      [,59]       [,60]
## [1,] -26.0813801 -29.3002580 -13.694444 -27.303675 -34.392021 -29.5555780
## [2,] -17.9403636  -6.6340685  -9.756045 -14.361541 -10.667489 -12.7682180
## [3,]  -0.2987818   0.8927321   1.158626  -3.970342  -2.332922  -0.9408487
## [4,]  13.4799099   9.3256102  12.056417  14.656972   8.910081  10.9351656
## [5,]  43.0863844  32.2156395  20.806965  54.656972  34.235293  32.6018322
##           [,61]      [,62]      [,63]      [,64]      [,65]       [,66]
## [1,] -0.5170388 -11.679803 -18.601864 -34.855419 -43.612766 -27.0502554
## [2,]  0.9988586  -6.610586 -10.953711 -17.687769 -10.730591 -10.2659729
## [3,]  2.5147560   1.563870  -2.621541   2.190348   2.443037   0.7414039
## [4,] 14.8720157   5.034477   8.261565  17.014070  13.387466  14.5653629
## [5,] 27.2292753  19.002023  33.358783  42.434670  39.059052  37.1647029
##           [,67]      [,68]      [,69]       [,70]       [,71]      [,72]
## [1,] -38.652064 -30.978237 -25.911421 -34.9602484 -28.3848666 -21.849092
## [2,] -13.639984  -9.506634 -11.795690  -9.0922921  -7.2084478  -9.769018
## [3,]  -3.611335  -2.643204  -2.677633   0.7265446  -0.9668129   2.072203
## [4,]  14.526399   7.774103   7.616348   8.1759494   7.2158709   6.974985
## [5,]  33.153713  32.498646  34.513120  31.5880632  25.7652061  15.713867
##           [,73]     [,74]      [,75]      [,76]     [,77]      [,78]
## [1,] -25.113274 -34.30052 -28.950313 -25.668912 -30.48747 -38.342644
## [2,] -10.279051 -11.91133 -15.153103 -10.679492 -12.96396 -11.451129
## [3,]  -1.353371  -5.26880  -3.773168  -2.140981  -4.18779  -2.263905
## [4,]   5.789305  10.75938   3.948561   6.287276  10.33352  11.627547
## [5,]  23.447808  37.14907  19.276965  25.997939  41.10619  40.509518
##            [,79]       [,80]       [,81]      [,82]       [,83]      [,84]
## [1,] -28.0116720 -54.8531052 -33.0155456 -38.326025 -37.9511934 -29.343274
## [2,] -10.9541317 -13.6107551  -6.3477156 -11.722611 -12.5540556 -14.438165
## [3,]  -0.2830986   0.7555571   0.3177877  -2.507336  -0.4696289  -1.301594
## [4,]  12.6131307  14.4195582  11.8872364  11.219523  15.1183327  17.128957
## [5,]  44.6494031  48.1557062  22.4751979  31.610722  45.9007302  41.441451
##             [,85]      [,86]      [,87]      [,88]      [,89]      [,90]
## [1,] -39.19606225 -48.115308 -38.174337 -50.810349 -48.906621 -16.233495
## [2,] -12.33710468 -11.601731  -9.928634 -17.058932 -11.685656  -7.606181
## [3,]   0.01590493   2.292338   1.498792   1.990095   5.471754  -4.860809
## [4,]  14.91868694  13.807561  11.249965  11.373946  14.490267   3.199027
## [5,]  44.03558121  39.086805  29.619205  53.896991  37.823600  16.708639
##            [,91]      [,92]      [,93]      [,94]       [,95]       [,96]
## [1,] -49.5937826 -29.862785 -33.669614 -32.272088 -34.1951658 -36.6991006
## [2,] -13.4459038 -15.086817 -21.414986 -16.107402 -10.1754895 -10.3005212
## [3,]   0.5737725  -4.302092  -4.454338  -2.385596   0.0224367   0.6077885
## [4,]  14.3966860  12.364574  16.919631   2.743699   7.7666451   7.7388314
## [5,]  55.4981610  50.991888  39.075758  20.846487  32.0803021  30.4842031
##           [,97]      [,98]       [,99]     [,100]      [,101]     [,102]
## [1,] -37.179738 -23.538235 -31.7906731 -27.528420 -51.6709758 -38.892472
## [2,] -11.885757  -9.266532  -8.2550057 -13.734199 -14.8114839 -11.323577
## [3,]  -2.864918  -2.137869   0.7602626  -7.136057  -0.2672417   1.812994
## [4,]   9.246497   7.596836   8.4060896   8.688074  15.9446059  10.440308
## [5,]  34.353469  17.343252  33.2295364  35.804914  33.5784999  35.502088
##         [,103]     [,104]     [,105]     [,106]     [,107]     [,108]
## [1,] -11.14048 -29.286176 -39.079962 -34.831170 -40.093377 -50.005146
## [2,] -11.14048 -11.834785 -11.237372 -13.055171 -12.952869 -16.183640
## [3,] -11.14048  -1.834785   2.540935  -4.231292  -3.918795  -1.994285
## [4,] -11.14048   9.343465  12.399793  11.226752  12.543068  16.476875
## [5,] -11.14048  34.440683  46.605218  37.383663  50.774298  61.195932
##           [,109]      [,110]     [,111]    [,112]
## [1,] -53.5888880 -34.3598915 -30.426923 -51.72400
## [2,] -14.1038204 -10.0033335 -11.411327 -17.22989
## [3,]   0.4772987   0.7373265  -1.712008  -3.02332
## [4,]  15.6841296  13.9793961  10.992638  15.27388
## [5,]  54.4138631  33.6771343  34.770703  52.48199
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]      [,4]      [,5]       [,6]
## [1,] -4.597524 -4.617369 -6.178301 -15.23414 -3.673464 -12.803859
## [2,]  1.400335  4.227441 13.074068   2.18853  9.364546   8.464331
##           [,7]       [,8]      [,9]     [,10]      [,11]     [,12]
## [1,] -13.59287 -10.885543 -1.347035 -3.873560 -16.582118 -3.898941
## [2,]  15.29925   5.210664  4.679666  5.203037   5.560071  4.710108
##           [,13]     [,14]    [,15]      [,16]       [,17]      [,18]
## [1,] -0.5051338 0.2033893 -2.12330 -0.3651225 -10.6001680 -10.405640
## [2,]  3.2708442 4.1337006 10.18858  6.6714984   0.3241369   3.027503
##           [,19]      [,20]      [,21]     [,22]     [,23]     [,24]
## [1,] -12.485556 -14.789577 -27.247330 -3.320726 -7.276725 -2.671172
## [2,]  -3.236039   1.808383   6.656246  1.406407  7.583976  8.552571
##          [,25]     [,26]     [,27]     [,28]      [,29]     [,30]
## [1,] -6.949270 -7.863287 -6.290278 -2.414225 -0.6819737 -5.321633
## [2,]  6.501495  7.000530  8.024273  7.932728 11.1524527  3.241580
##            [,31]      [,32]     [,33]     [,34]      [,35]      [,36]
## [1,] -13.7914885 -13.601022 -4.748219 -9.558493 -17.997973 -16.934157
## [2,]  -0.1860494   4.909704  7.783327  3.683969   7.399162  -4.716122
##          [,37]     [,38]     [,39]     [,40]     [,41]     [,42]     [,43]
## [1,] -49.89080 -5.069954 -9.927011 -2.856005 -5.113650 -9.603108 -7.416519
## [2,]  16.51994  7.489740  5.113518  3.477467  2.865377  2.156576  2.557419
##            [,44]     [,45]    [,46]     [,47]     [,48]     [,49]
## [1,] -10.2731092 -4.769228 -3.28209 -7.764618 -3.749824 -5.649300
## [2,]   0.2369567 11.045174 10.21937  1.213391  2.652261  4.735577
##          [,50]     [,51]      [,52]     [,53]     [,54]     [,55]
## [1,] -6.150998 -6.085036 -9.6055550 -4.353972 -1.897141 -8.812663
## [2,]  6.494164  3.489233  0.4365524  6.685665 12.752822  8.215100
##          [,56]     [,57]      [,58]     [,59]     [,60]     [,61]
## [1,] -1.980931 -11.86742 -13.140192 -6.893679 -6.346485 -10.14052
## [2,]  3.766395  14.18468   5.199508  2.227835  4.464788  15.17004
##          [,62]     [,63]    [,64]      [,65]     [,66]      [,67]
## [1,] -3.747521 -9.777493 -7.35414 -0.6787719 -6.194166 -10.478279
## [2,]  6.875261  4.534412 11.73483  5.5648455  7.676974   3.255609
##           [,68]      [,69]     [,70]     [,71]     [,72]     [,73]
## [1,] -5.7547342 -5.9106421 -2.304991 -3.985478 -3.444156 -5.137991
## [2,]  0.4683256  0.5553768  3.758080  2.051852  7.588561  2.431249
##           [,74]     [,75]     [,76]      [,77]      [,78]     [,79]
## [1,] -11.808558 -11.09305 -6.051256 -10.409819 -10.221069 -4.499276
## [2,]   1.270958   3.54671  1.769293   2.034239   5.693259  3.933079
##          [,80]     [,81]     [,82]     [,83]      [,84]    [,85]
## [1,] -4.615139 -3.884761 -7.910954 -7.216138 -11.276805 -6.19987
## [2,]  6.126253  4.520337  2.896282  6.276880   8.673616  6.23168
##           [,86]     [,87]     [,88]    [,89]     [,90]     [,91]
## [1,] -0.7867715 -3.382170 -5.113004 -1.15082 -7.746543 -6.562557
## [2,]  5.3714483  6.379753  9.093195 12.09433 -1.975076  7.710102
##           [,92]      [,93]     [,94]     [,95]     [,96]     [,97]
## [1,] -11.432602 -16.110789 -12.91609 -2.565422 -4.353804 -8.211431
## [2,]   2.828417   7.202112   8.14490  2.610295  5.569381  2.481595
##          [,98]     [,99]     [,100]    [,101]    [,102]    [,103]
## [1,] -9.829364 -1.613340 -12.359512 -8.857638 -3.840184 -11.14048
## [2,]  5.553626  3.133866  -1.912603  8.323155  7.466171 -11.14048
##         [,104]    [,105]     [,106]    [,107]     [,108]    [,109]
## [1,] -5.863093 -2.240822 -7.9576731 -6.156772 -4.9490994 -3.989910
## [2,]  2.193522  7.322693 -0.5049115 -1.680819  0.9605291  4.944507
##         [,110]    [,111]    [,112]
## [1,] -5.492197 -9.258944 -7.965049
## [2,]  6.966850  5.834927  1.918408
## 
## $out
##   [1]  79.59584 -58.78642  61.85259  48.09907 -85.31506 -49.31245 -46.84232
##   [8] -52.19405  49.30989 -39.80109 -47.23072 -45.87177 -44.53999 -70.12637
##  [15] -47.42996 -56.80315 -42.52026 -44.56104 -45.14494 -41.53373 -36.54770
##  [22]  43.24252 -40.06476 -54.23187 -36.70701 -43.49438 -40.35548  36.17173
##  [29] -30.72340 -32.72340  40.30549  45.06170  68.80587 -32.96559 -14.32382
##  [36]  26.25836  59.76010  57.52391  50.88167 -33.55722  51.76079  48.88254
##  [43]  45.00251  43.42856  60.69749 -41.72437 -58.46840 -36.12291  39.31332
##  [50] -51.14726  48.27112  34.05645  39.87316  57.56501 -61.81680  56.33230
##  [57] -46.76999  38.42241 -46.66293 -41.71311 -33.07176 -60.04672 -24.30112
##  [64]  49.51421 -50.15653  62.81226  34.31416  51.26782  58.13875 -35.64072
##  [71]  32.54272 -33.22274  29.68249  29.26378  32.66358  34.55296  32.90971
##  [78]  30.89687 -40.79338  32.17384  46.75669  47.03186 -48.29451 -34.77943
##  [85] -57.36717  50.62923  48.08063 -66.29291 -57.71080 -45.32999  25.33595
##  [92]  39.64845 -26.44031  23.91685 -37.62984  41.29558 -37.91970  40.48420
##  [99]  51.66013  46.54479 -53.40967  49.80095 -33.99703  53.25630  44.43223
## [106] -53.40076  46.00168  50.70770  62.47158 -54.76514 -47.67349  51.44713
## [113]  54.90508  54.90508 -53.48593 -66.72693 -71.31379
## 
## $group
##   [1]   1   1   1   9   9   9  10  12  13  13  13  13  13  14  14  14  14
##  [18]  14  14  14  14  14  14  14  14  14  14  16  16  16  17  17  19  19
##  [35]  22  22  23  23  24  30  31  31  31  31  32  40  41  42  42  43  43
##  [52]  43  43  44  45  47  47  48  48  51  56  58  62  68  68  68  68  69
##  [69]  69  70  71  71  71  71  72  73  73  73  73  76  77  77  79  81  82
##  [86]  82  82  86  86  87  90  90  90  90  95  95  95  96  96  97  97  99
## [103]  99 100 100 100 100 100 100 105 105 107 107 107 107 108 109
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##            [,1]       [,2]       [,3]        [,4]       [,5]       [,6]
## [1,] -52.975368 -42.483623 -31.615075 -24.1785464 -38.891781 -40.268728
## [2,] -14.766890 -14.062275 -11.004036 -11.6711175  -8.642923  -9.174204
## [3,]  -1.071269   1.255501   2.078377  -0.9755958   1.273921  -5.638004
## [4,]  15.762654  15.964593  16.524297   9.8548419  11.734454  17.750130
## [5,]  58.253433  34.796875  40.066902  32.1909896  33.941314  41.746862
##            [,7]       [,8]       [,9]       [,10]      [,11]      [,12]
## [1,] -33.594717 -20.881351 -35.665085 -43.4534392 -35.780627 -48.125763
## [2,]  -8.271240 -14.695615  -8.593425 -14.1791582 -21.367179 -13.556462
## [3,]   2.527303  -9.530361   1.267939   0.3271607   1.948445   4.343739
## [4,]  13.676325   6.873214  12.264459  15.1159080  16.632979  16.478707
## [5,]  29.841829  25.215415  39.131937  44.0730204  31.074670  37.889143
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -29.762224 -40.871046 -31.349702 -26.881556 -28.769073 -40.394753
## [2,]  -8.994243  -9.243719 -11.608194  -8.598017 -13.594922 -14.345708
## [3,]   1.981624   3.590450   4.192207  -1.950826  -4.695103   2.067992
## [4,]   9.281876  12.840619  16.604475   5.476874  11.572090  14.245293
## [5,]  34.339091  43.604565  32.701006  18.569815  40.548398  30.709093
##           [,19]      [,20]      [,21]     [,22]       [,23]      [,24]
## [1,] -21.976171 -42.699093 -23.336613 -5.560738 -43.4632765 -33.677497
## [2,] -15.134039 -18.709848 -16.182013 -1.815781 -12.5970113  -9.704370
## [3,] -11.712047  -4.559494  -3.298997  2.246199  -0.7000782   1.908978
## [4,]   7.590295  13.483572   5.863983  6.029353  16.0213895  10.731628
## [5,]  26.037249  33.997356  17.639254 16.042912  39.3034066  40.206824
##           [,25]      [,26]      [,27]      [,28]      [,29]       [,30]
## [1,] -37.237884 -17.272519 -42.799667 -27.367180 -32.085278 -24.7553276
## [2,] -14.292975  -7.684720 -15.391056  -8.856320  -6.756859  -8.0576627
## [3,]   2.834395  -1.700300   5.590011   3.149691   5.425142  -0.2239945
## [4,]  15.794544   1.470838  15.900944  11.686415  12.091808   7.8633266
## [5,]  41.588898   9.171424  44.771327  38.068014  29.260744  25.1933924
##          [,31]      [,32]       [,33]      [,34]      [,35]     [,36]
## [1,] -53.94805 -28.078671 -36.0683244 -36.615112 -40.136156 -25.21884
## [2,] -19.74636 -13.790458  -8.8608544 -16.169791 -11.672879 -25.21884
## [3,]  -4.67079  -3.463607   0.8267174  -1.915391  -2.397263 -17.80452
## [4,]  10.56925  14.308342  10.1378092  14.345543  15.920087 -10.39020
## [5,]  49.10861  36.080053  34.3922086  46.703009  22.229061 -10.39020
##            [,37]        [,38]      [,39]      [,40]      [,41]      [,42]
## [1,] -36.4611028 -20.46513015 -16.250172 -30.513364 -50.646200 -32.715797
## [2,] -36.4611028  -7.85124159  -9.157051  -9.607160 -14.505953 -13.284959
## [3,] -18.6875723  -0.01308331  -3.025096   0.303661   1.700181  -6.380527
## [4,]  -0.9140417  13.12682983   6.939513   8.565570  13.248268   3.899480
## [5,]  -0.9140417  38.56101491  29.495429  32.359437  44.503411  13.674274
##           [,43]      [,44]      [,45]      [,46]      [,47]       [,48]
## [1,] -27.281505 -48.894179 -30.404198 -30.549565 -36.047808 -37.7030032
## [2,] -13.097480 -14.157698 -14.826364 -12.243632 -12.422824 -10.0688123
## [3,]  -2.977890  -5.208029  -2.102355   3.704834  -2.584391  -0.1008466
## [4,]   5.679052  10.761056  12.304989   9.017302   9.222880   9.5210547
## [5,]  22.089267  37.562786  46.884940  21.807901  29.504689  30.8750205
##           [,49]       [,50]        [,51]      [,52]      [,53]      [,54]
## [1,] -27.694523 -39.3813376 -30.59405831 -31.470689 -33.465625 -30.987994
## [2,] -11.718656 -11.7936048  -7.77028891  -9.316940  -7.452333  -9.453667
## [3,]  -3.155056   0.3609954  -0.08979693  -5.007739   1.010803   0.598394
## [4,]   9.514878  13.4914627   9.37854122  14.911877  13.195907  13.292716
## [5,]  22.415079  32.4607947  29.35114074  29.707193  41.040794  36.298076
##           [,55]       [,56]      [,57]      [,58]       [,59]       [,60]
## [1,] -41.836257 -32.6889607 -29.229231 -59.854488 -41.4680677 -38.3759475
## [2,] -14.119341  -8.4907617  -4.853725 -14.679023  -9.9013556 -17.4453538
## [3,]  -1.842946   0.1722116   4.768139  -4.733824   0.9300629  -0.1313211
## [4,]  10.233873  11.0932777  18.250167  20.551244  13.1157526  13.8068034
## [5,]  28.227338  33.5888799  39.414586  53.829776  33.5214847  35.8065866
##            [,61]       [,62]      [,63]       [,64]      [,65]       [,66]
## [1,] -11.4391730 -17.1788864 -19.985541 -35.0627581 -45.318289 -35.0275871
## [2,]  -7.3428775  -7.5313639 -11.827673 -12.5665811 -11.528328  -8.7099849
## [3,]  -3.2465821   0.2059804   2.829742   0.4773994   0.888533   0.1483081
## [4,]   0.6196307   6.3918652   7.196460  11.7807318  13.511719  14.3915860
## [5,]   4.4858435   9.3129836  31.911393  40.0282209  37.617513  43.4102521
##           [,67]       [,68]       [,69]      [,70]      [,71]        [,72]
## [1,] -29.376805 -36.3164489 -29.3709456 -30.822993 -28.729728 -19.97878997
## [2,] -16.549037 -11.0479097  -9.4911395  -9.762167 -10.082771 -11.89885658
## [3,]  -2.958037  -0.9468153  -0.4944071  -1.253368  -1.640213   0.05161738
## [4,]  16.839097   7.9762233   9.0176594   6.389166  12.632131   8.66701134
## [5,]  36.063159  30.6645696  29.5329934  26.058337  29.215383  19.94398493
##            [,73]      [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -27.9217054 -26.854735 -28.276277 -26.540917 -23.741826 -29.400482
## [2,] -10.7975486 -10.576833 -24.148846 -11.155871 -11.718016 -13.664366
## [3,]  -0.7975486  -4.015084  -8.080971  -3.637495  -1.287158  -2.026386
## [4,]   6.8449851   6.906252  11.064937   6.228996   4.764205   8.326003
## [5,]  31.1609465  32.295967  44.003958  31.561195  28.339719  40.334865
##            [,79]       [,80]      [,81]      [,82]      [,83]      [,84]
## [1,] -30.7236197 -35.2323159 -34.385920 -31.799458 -53.603500 -31.858900
## [2,] -12.2075404 -11.3778203  -8.530172 -12.853672 -13.051310 -15.488363
## [3,]  -0.3062288  -0.8230332   0.329012  -4.488977  -3.820169   2.505102
## [4,]  13.8039173  11.6021892  12.511013   9.152863  17.581321  21.209322
## [5,]  39.2572994  32.1403148  22.741279  36.821365  48.654727  39.116967
##           [,85]      [,86]      [,87]      [,88]      [,89]      [,90]
## [1,] -37.353527 -41.326862 -44.454341 -48.642560 -36.325201 -16.950722
## [2,] -14.225962 -12.925444 -12.955854 -16.997509 -13.682667  -7.399963
## [3,]  -2.040144   2.697659   0.379911   1.752994   3.214266  -5.514322
## [4,]  13.241106  14.159919   9.337040  14.938263  13.457030   3.049278
## [5,]  39.337638  47.431232  35.071071  46.089595  35.571732  16.437412
##           [,91]      [,92]      [,93]      [,94]      [,95]      [,96]
## [1,] -30.602877 -28.801316 -37.094809 -23.095846 -35.686720 -57.119752
## [2,] -11.666861 -14.130411 -18.761476 -16.990403  -9.449752 -13.276002
## [3,]   1.167473  -2.069948   4.911726  -4.049624   1.350649   1.245341
## [4,]  13.162416  11.045055  17.608693   1.987423   8.220182  19.122555
## [5,]  34.078739  26.330122  35.777991  12.010344  32.501982  47.912008
##          [,97]      [,98]       [,99]     [,100]      [,101]     [,102]
## [1,] -35.88038 -27.598920 -32.1919465 -38.123685 -41.2751123 -37.555575
## [2,] -12.34418  -8.381802  -9.0255068 -15.999337 -13.5198603 -11.682992
## [3,]  -4.70864  -2.163242   0.6808537  -5.249700  -0.7739482   1.189809
## [4,]   9.32249   7.672681   9.4451775   8.905981  15.2040070   7.694588
## [5,]  33.31514  25.837214  36.7278386  44.055464  46.8387872  26.972253
##         [,103]     [,104]       [,105]     [,106]     [,107]     [,108]
## [1,] -10.78964 -28.939746 -40.46937911 -35.393965 -46.498035 -52.105066
## [2,] -10.78964 -11.297213 -10.98471313 -13.397755 -13.038862 -15.699790
## [3,] -10.78964  -3.194146  -0.03632628  -3.858266  -2.232676  -2.013436
## [4,] -10.78964   9.623854  11.99768907  10.559436  12.951479  16.256899
## [5,] -10.78964  35.424255  45.56096029  37.960032  50.896321  55.433372
##           [,109]     [,110]    [,111]     [,112]
## [1,] -54.6137055 -36.644436 -28.57179 -41.137548
## [2,] -16.3630478  -8.030733 -12.43238 -15.985096
## [3,]  -0.3031973   1.030983  -1.60414  -1.494281
## [4,]  19.4762484  13.958585  10.41894  14.183543
## [5,]  54.3591168  31.952050  33.42430  54.521627
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]      [,4]      [,5]       [,6]
## [1,] -4.251902 -3.920895 -8.795315 -8.580690 -7.655714 -16.273116
## [2,]  2.109363  6.431897 12.952068  6.629499 10.203556   4.997108
##           [,7]      [,8]      [,9]     [,10]     [,11]     [,12]     [,13]
## [1,] -7.090409 -24.77084 -1.431877 -6.154204 -13.06162 -1.411094 0.1589705
## [2,] 12.145015   5.71012  3.967755  6.808526  16.95851 10.098572 3.8042782
##         [,14]     [,15]     [,16]      [,17]     [,18]      [,19]
## [1,] 1.374726 -4.085351 -5.944948 -11.954961 -7.153068 -20.982536
## [2,] 5.806174 12.469766  2.043295   2.564754 11.289051  -2.441559
##           [,20]      [,21]      [,22]     [,23]     [,24]     [,25]
## [1,] -15.404077 -16.464510 -0.9542567 -8.571360 -2.958753 -4.877332
## [2,]   6.285089   9.866516  5.4466549  7.171204  6.776709 10.546122
##          [,26]     [,27]     [,28]      [,29]     [,30]      [,31]
## [1,] -4.784415 -2.131433 -1.743464 -0.7846041 -4.416522 -13.273646
## [2,]  1.383815 13.311454  8.042846 11.6348874  3.968533   3.932067
##           [,32]     [,33]     [,34]     [,35]      [,36]     [,37]
## [1,] -12.007641 -4.653774 -9.841762 -15.54222 -34.371499 -58.40177
## [2,]   5.080427  6.307209  6.010979  10.74769  -1.237532  21.02663
##          [,38]     [,39]     [,40]     [,41]       [,42]     [,43]
## [1,] -6.778850 -8.574938 -2.757146 -2.269963 -12.6094880 -8.394304
## [2,]  6.752684  2.524746  3.364468  5.670325  -0.1515652  2.438524
##           [,44]     [,45]    [,46]     [,47]     [,48]     [,49]     [,50]
## [1,] -12.061751 -9.680336 -3.29964 -7.327107 -3.539957 -8.459616 -6.490443
## [2,]   1.645692  5.475627 10.70931  2.158325  3.338263  2.149504  7.212434
##          [,51]      [,52]     [,53]     [,54]     [,55]     [,56]
## [1,] -4.428494 -11.217823 -4.213253 -5.014379 -8.441882 -3.354046
## [2,]  4.248900   1.202345  6.234859  6.211167  4.755990  3.698469
##          [,57]      [,58]     [,59]     [,60]      [,61]     [,62]
## [1,] -9.029133 -15.866588 -4.431963 -7.258480 -10.510089 -6.144498
## [2,] 18.565410   6.398941  6.292089  6.995838   4.016925  6.556459
##          [,63]     [,64]     [,65]     [,66]      [,67]     [,68]
## [1,] -4.255027 -6.219151 -2.352617 -6.304127 -11.098040 -4.372257
## [2,]  9.914511  7.173950  4.129683  6.600743   5.181967  2.478627
##          [,69]     [,70]     [,71]     [,72]     [,73]     [,74]
## [1,] -3.576985 -4.088824 -6.393899 -6.723865 -4.952939 -9.058381
## [2,]  2.588171  1.582088  3.113473  6.827100  3.357841  1.028213
##           [,75]      [,76]     [,77]     [,78]     [,79]     [,80]
## [1,] -21.575114 -7.6441272 -5.689044 -9.608319 -4.959673 -5.226075
## [2,]   5.413172  0.3691381  3.114728  5.555547  4.347215  3.580009
##          [,81]      [,82]      [,83]     [,84]     [,85]      [,86]
## [1,] -4.520280 -9.6722313 -11.288383 -9.091367 -8.304101 -0.5845578
## [2,]  5.178304  0.6942769   3.648045 14.101570  4.223813  5.9798759
##          [,87]     [,88]     [,89]     [,90]     [,91]     [,92]     [,93]
## [1,] -4.757859 -6.225198 -3.652145 -8.304988 -5.196515 -8.609285 -6.147393
## [2,]  5.517681  9.731187 10.080677 -2.723657  7.531461  4.469389 15.970845
##           [,94]     [,95]     [,96]       [,97]     [,98]     [,99]
## [1,] -14.650910 -1.197949 -7.665645 -10.1903604 -9.485799 -1.950550
## [2,]   6.551662  3.899248 10.156327   0.7730805  5.159316  3.312258
##        [,100]    [,101]    [,102]    [,103]     [,104]    [,105]
## [1,] -11.0516 -8.796730 -3.843525 -10.78964 -7.1735344 -4.685626
## [2,]   0.5522  7.248834  6.223143 -10.78964  0.7852428  4.612974
##          [,106]      [,107]     [,108]    [,109]    [,110]   [,111]
## [1,] -7.5348125 -4.51405053 -4.9045752 -5.677908 -4.680751 -9.30177
## [2,] -0.1817199  0.04869813  0.8777024  5.071513  6.742717  6.09349
##         [,112]
## [1,] -6.080987
## [2,]  3.092425
## 
## $out
##   [1]  67.47810 -47.91383  47.15812 -43.55283 -84.75530 -48.13962  46.39679
##   [8] -38.24724 -51.99871  37.67242 -51.29309  42.37838  42.74237 -47.21244
##  [15] -44.68504 -43.76981 -44.09865 -43.70183 -42.87758 -64.97848 -45.89395
##  [22] -51.90434 -44.77451 -46.37190  46.04725  36.19507  37.24527  59.86596
##  [29] -18.65990 -26.40004 -62.80132  61.19686  49.75576  33.20866  24.45649
##  [36] -24.30588  21.17376  25.42677 -36.05884  59.43033  60.07222  41.51941
##  [43]  34.26516  33.30888  39.74221  45.61145 -41.89496  35.31543  46.60012
##  [50]  48.83086  44.92790  38.97809 -49.51950  51.51834 -48.01947 -49.24727
##  [57]  46.79094  53.23380 -50.92665  54.69769  41.61566 -44.07452  36.96700
##  [64]  34.09092 -39.14477 -39.42212  32.69280  41.36479  33.44017  36.59719
##  [71]  45.28231  47.60686  45.24939 -56.66817 -64.47799  25.00101  36.22109
##  [78]  38.74008 -35.54569  50.52972  41.52611 -37.49802 -56.19056  47.39816
##  [85] -46.12778 -37.73101 -38.20355 -39.38764  42.46731  40.61005  51.18266
##  [92]  48.82520  60.20680 -45.47071 -48.88523 -55.64484  53.45694  54.32884
##  [99]  54.32884  65.52114
## 
## $group
##   [1]   1   4   4   9   9   9   9  13  13  13  13  13  13  13  13  14  14
##  [18]  14  14  14  14  14  14  14  14  16  16  19  22  22  23  23  24  26
##  [35]  26  26  26  26  26  31  32  33  39  39  40  40  40  42  43  44  47
##  [52]  48  48  49  51  54  55  62  65  68  68  69  70  70  73  74  76  76
##  [69]  77  77  77  82  82  86  86  90  90  90  90  92  95  95  95  97  99
##  [86]  99  99  99  99  99 100 100 100 105 105 105 107 107 107 108
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##             [,1]       [,2]       [,3]        [,4]          [,5]
## [1,] -44.9195591 -49.508220 -12.299130 -29.9603882 -14.279046065
## [2,] -13.5622502 -15.088586  -4.347759 -10.5061891  -8.048586501
## [3,]  -0.7675704  -1.364823  -1.107263  -0.4403256   0.003209472
## [4,]  14.3612828  15.779121   6.423455   7.1836843   8.912400560
## [5,]  53.5438592  44.925984  14.754567  32.4115424  24.549731664
##             [,6]       [,7]       [,8]       [,9]      [,10]     [,11]
## [1,] -21.4629785 -34.822029 -23.112709 -34.911666 -52.757100 -33.26434
## [2,] -13.0446182 -13.428316 -17.536071  -9.544144 -16.237094 -15.56088
## [3,]  -0.7365461   1.280924  -3.604942   1.912969   2.090976  -2.94235
## [4,]  16.5277781   5.118820  -2.200021  12.278328  10.912840  11.16728
## [5,]  40.4174900  26.369438  18.685346  39.937903  33.971445  20.80717
##          [,12]      [,13]      [,14]      [,15]      [,16]      [,17]
## [1,] -43.25572 -34.120025 -40.890916 -29.782889 -31.290161 -37.341740
## [2,] -14.57905  -8.309608  -8.396345 -11.663358 -11.080158 -12.942688
## [3,]   1.45869   1.927182   2.415386   4.318156   1.223954  -5.459741
## [4,]  15.24919   9.704137  13.341372  16.796282   6.736585   9.229369
## [5,]  55.26811  36.713380  40.732052  32.695237  20.346749  39.429986
##           [,18]      [,19]      [,20]      [,21]      [,22]       [,23]
## [1,] -27.820821 -19.554785 -37.732089 -20.022271 -13.433952 -45.6557988
## [2,] -14.060483 -14.985403 -15.689455 -15.410298  -3.906793 -14.9706455
## [3,]  -3.673715 -11.084600  -8.207341  -9.067062  -0.545234  -0.1844473
## [4,]  13.760293  -4.618517   7.533236  -3.783208   3.956814  13.7465436
## [5,]  31.153762   1.159666  32.700839   1.002990   8.361135  41.3578748
##           [,24]      [,25]      [,26]       [,27]      [,28]      [,29]
## [1,] -30.321364 -52.396638 -38.366681 -30.3873239 -27.680728 -15.177936
## [2,] -10.188291 -12.133473  -8.302198 -14.8392848  -6.397007  -8.810098
## [3,]   2.902168   1.243949  -3.322933  -0.3206393   2.844510   6.444983
## [4,]  11.806223  18.534462  18.218187  15.7667647  10.998655  13.325452
## [5,]  37.200510  35.707570  31.824194  45.2810602  30.866407  31.372548
##           [,30]      [,31]      [,32]       [,33]      [,34]      [,35]
## [1,] -24.621876 -40.987917 -33.080149 -34.7611604 -46.851592 -34.775719
## [2,]  -7.923616 -17.802171 -12.320799  -9.9749622 -15.541297 -19.706984
## [3,]   2.171097  -2.024872  -6.444051   0.9432697  -1.610306   4.965546
## [4,]   7.586599  11.461775  12.489648  16.0940469  18.389694  13.305373
## [5,]  26.935339  50.620741  42.208422  42.7607136  47.534487  16.742231
##           [,36]     [,37]      [,38]       [,39]       [,40]       [,41]
## [1,] -13.409309 -8.091771 -29.160324 -21.1610142 -39.3111896 -46.3085450
## [2,] -13.409309 -8.091771  -6.847114  -9.2481087  -9.3962165 -13.0901459
## [3,] -10.859033 -2.183585   5.166119   0.1542337   0.8175853  -0.3643807
## [4,]  -8.308758  3.724601  16.129231   4.8172336  10.8542609  13.1489172
## [5,]  -8.308758  3.724601  27.390729  15.0930753  40.4507399  37.0703380
##           [,42]      [,43]      [,44]      [,45]      [,46]      [,47]
## [1,] -30.337018 -38.197969 -34.481125 -44.065226 -29.918139 -40.529813
## [2,] -11.790727 -10.592405 -13.692323 -15.371388 -11.932945 -11.249463
## [3,]  -5.669079  -1.436936  -2.409512   2.649933   1.962330  -2.144495
## [4,]   4.236086   8.122261   8.576025  13.568977   9.221477  12.017216
## [5,]  19.244567  18.969673  34.180449  55.334954  22.389933  32.154521
##           [,48]      [,49]       [,50]       [,51]      [,52]       [,53]
## [1,] -28.891450 -26.565654 -38.0442299 -36.6402261 -25.079545 -34.3105024
## [2,]  -8.900120 -10.540393 -12.4028245  -9.6147879  -9.865744  -7.1146082
## [3,]  -2.233453  -1.993258   0.4178783   0.5832641  -4.046858   0.4391695
## [4,]   7.766547  10.399841  14.6039501  11.4637068  13.433273  14.1340634
## [5,]  32.292713  23.774454  33.6661847  31.6775086  29.108995  30.7195839
##            [,54]      [,55]      [,56]      [,57]      [,58]       [,59]
## [1,] -27.3331743 -24.905762 -28.721769 -21.636749 -58.170063 -35.4928245
## [2,]  -9.6543239 -14.853963  -6.927680  -3.986993 -17.484909 -11.7350163
## [3,]  -0.1379106  -1.980275   0.594194   1.675254  -3.296369   0.2336079
## [4,]  13.1482060  13.514274  13.106750  16.811405  15.250766  10.5098374
## [5,]  35.9538877  42.616112  28.767942  24.879329  55.250766  33.0543107
##           [,60]      [,61]       [,62]       [,63]      [,64]      [,65]
## [1,] -36.780703 -16.997837 -20.7012283 -19.4258983 -33.162019 -36.614761
## [2,] -14.428539 -15.465374 -11.6568174 -11.1363129 -15.696257 -12.299802
## [3,]   1.695103 -13.932911  -0.1433574   0.8848838  -7.137346  -1.991512
## [4,]  13.510746  -4.839572   7.0486011   7.8384259  13.755036  12.439395
## [5,]  32.934166   4.253768  29.8812470  32.4545702  44.644265  46.840519
##            [,66]       [,67]      [,68]       [,69]       [,70]
## [1,] -32.1991265 -36.9457648 -26.664929 -35.7200413 -28.2899592
## [2,]  -9.9879698  -9.1068016 -11.594133 -10.6257506  -9.3890259
## [3,]   0.7619659   0.5353764  -1.166529  -0.8989444   0.4409202
## [4,]  11.7384015  14.1390359   7.808210   7.9213846   7.7052444
## [5,]  27.6758668  30.8931984  30.360129  29.6450780  31.7325225
##            [,71]      [,72]      [,73]       [,74]       [,75]      [,76]
## [1,] -22.5581294 -21.907342 -26.252983 -32.6875321 -42.9813393 -27.640085
## [2,]  -6.1548636 -10.111900  -8.979956 -17.0418676 -17.8875007 -11.263776
## [3,]  -0.3474055   1.853595  -3.250758  -0.6444172   0.2097993  -3.019923
## [4,]   5.2980012   7.322478   6.603836  13.1281863  14.3330756   5.421669
## [5,]  21.8527855  21.867580  29.339512  38.1885459  27.3460645  27.538860
##           [,77]      [,78]      [,79]      [,80]       [,81]      [,82]
## [1,] -42.839234 -35.961930 -34.613500 -38.202454 -33.6509347 -34.205310
## [2,] -14.383173 -12.904033  -9.443623 -12.645651 -10.4026283 -11.592529
## [3,]  -2.545647  -2.451289  -1.622252  -2.476021   0.5376059  -3.045394
## [4,]  13.461998   4.840570  12.162976  13.806789  13.2295338   8.835075
## [5,]  47.147997  26.854330  38.675908  41.090237  28.0589547  31.313201
##           [,83]       [,84]       [,85]      [,86]      [,87]      [,88]
## [1,] -59.721751 -29.1845277 -34.8792039 -47.860522 -38.661932 -42.743063
## [2,] -13.603892 -10.4701754 -12.7655359 -11.712603  -9.187967 -18.321551
## [3,]   1.258681   0.7424546  -0.3724368   1.679853   3.816194   2.484909
## [4,]  18.534126  18.2498504  13.6873289  14.548386  11.037513  13.648590
## [5,]  45.628397  46.8416404  43.6873289  43.658054  33.073572  44.929954
##           [,89]      [,90]      [,91]      [,92]        [,93]      [,94]
## [1,] -34.958852 -19.805199 -39.947862 -42.603289 -37.10127577 -35.372035
## [2,] -10.954713  -9.805199 -12.171129 -16.058166 -15.43208702 -15.217934
## [3,]  -2.015279  -4.591398   5.334073  -1.780416   0.08497085  -2.587538
## [4,]  13.588283   5.873706  12.187466  13.621489  13.38284345   2.013507
## [5,]  38.444757  22.245323  34.067934  50.180897  35.80445393   3.817344
##            [,95]       [,96]      [,97]      [,98]      [,99]     [,100]
## [1,] -38.1152841 -39.8744299 -33.639039 -28.704611 -32.251232 -29.998296
## [2,] -10.8072120 -12.3415501 -14.667602 -15.834612  -8.114912 -16.922512
## [3,]   0.8157068  -0.7186313  -5.349453  -4.623659   0.803708  -5.951428
## [4,]   7.8064874  14.3925687   8.264512   9.867347   8.655619  13.932695
## [5,]  31.8847159  49.5389184  42.550164  32.134967  30.973762  60.001704
##          [,101]     [,102]    [,103]     [,104]     [,105]     [,106]
## [1,] -38.555616 -29.886680 -12.32963 -28.815224 -44.354317 -34.478498
## [2,] -12.421637  -9.613752 -12.32963 -11.453386 -11.438584 -12.825575
## [3,]  -1.068063  -3.544743 -12.32963  -3.003765   1.039542  -4.745537
## [4,]  10.114786   6.965418 -12.32963   9.474361  12.064803   9.985263
## [5,]  40.525766  26.976412 -12.32963  39.369814  46.083290  39.513042
##          [,107]     [,108]      [,109]     [,110]      [,111]     [,112]
## [1,] -47.716092 -48.390232 -54.9007516 -38.169837 -29.4777411 -37.606114
## [2,] -13.346644 -16.285813 -16.2993906 -10.818017 -13.2752800 -16.656729
## [3,]  -1.741014  -1.216233  -0.2887214  -2.270882  -0.7264127  -1.440088
## [4,]  11.912333  15.227127  17.3391579  15.894293   9.9118852  13.793851
## [5,]  43.739091  57.199523  55.1213385  35.502490  32.4973250  53.501764
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
## [1,] -3.676703 -6.686173 -5.361892 -6.690136 -7.429318 -12.41764 -6.846674
## [2,]  2.141562  3.956527  3.147367  5.809484  7.435737  10.94455  9.408523
##            [,8]       [,9]     [,10]     [,11]     [,12]     [,13]
## [1,] -14.441357 -0.9117024 -3.915791 -13.49997 -4.256495 0.1306944
## [2,]   7.231472  4.7376402  8.097743   7.61527  7.173875 3.7236697
##          [,14]     [,15]     [,16]       [,17]      [,18]      [,19]
## [1,] 0.2344389 -4.031864 -3.832017 -11.8556528 -12.646365 -15.313815
## [2,] 4.5963337 12.668175  6.279925   0.9361706   5.298934  -6.855385
##            [,20]      [,21]     [,22]     [,23]     [,24]     [,25]
## [1,] -16.0300705 -16.010572 -3.753226 -8.082900 -2.336792 -6.616545
## [2,]  -0.3846123  -2.123551  2.662758  7.714006  8.141129  9.104442
##           [,26]     [,27]     [,28]      [,29]     [,30]      [,31]
## [1,] -12.256514 -7.872821 -1.299031 -0.8476351 -1.913260 -10.329288
## [2,]   5.610648  7.231543  6.988052 13.7376012  6.255453   6.279545
##          [,32]     [,33]      [,34]     [,35]      [,36]     [,37]
## [1,] -13.98819 -6.576785 -10.423895 -10.76114 -16.557516 -15.38518
## [2,]   1.10009  8.463324   7.203282  20.69223  -5.160551  11.01800
##          [,38]     [,39]     [,40]     [,41]       [,42]     [,43]
## [1,] -2.244123 -4.695275 -2.593174 -4.117786 -11.4784282 -6.835504
## [2,] 12.576361  5.003742  4.228344  3.389025   0.1402698  3.961631
##          [,44]     [,45]     [,46]     [,47]      [,48]     [,49]
## [1,] -8.534259 -5.433319 -5.007053 -7.242378 -5.1593794 -7.224546
## [2,]  3.715234 10.733185  8.931712  2.953387  0.6924724  3.238030
##          [,50]     [,51]      [,52]     [,53]     [,54]     [,55]
## [1,] -6.900087 -4.749647 -10.018626 -4.936797 -5.764539 -9.667152
## [2,]  7.735844  5.916176   1.924909  5.815136  5.488717  5.706603
##          [,56]     [,57]      [,58]     [,59]     [,60]      [,61]
## [1,] -3.013160 -10.74521 -13.640842 -4.948515 -4.676545 -23.625910
## [2,]  4.201548  14.09572   7.048105  5.415731  8.066751  -4.239911
##          [,62]     [,63]       [,64]     [,65]     [,66]     [,67]
## [1,] -8.675024 -6.181491 -15.2377089 -5.193720 -5.306366 -5.131942
## [2,]  8.388310  7.951258   0.9630173  1.210696  6.830297  6.202695
##          [,68]     [,69]     [,70]     [,71]     [,72]     [,73]     [,74]
## [1,] -4.660071 -3.987907 -2.560074 -2.744217 -3.890209 -6.921248 -9.347489
## [2,]  2.327012  2.190018  3.441914  2.049405  7.597399  0.419731  8.058654
##          [,75]      [,76]     [,77]     [,78]     [,79]     [,80]
## [1,] -12.13733 -6.8653622 -9.982221 -8.569349 -5.487668 -7.544392
## [2,]  12.55693  0.8255162  4.890928  3.666771  2.243164  2.592349
##          [,81]     [,82]     [,83]     [,84]     [,85]     [,86]
## [1,] -4.908820 -7.856758 -6.576545 -8.333074 -6.405101 -1.502466
## [2,]  5.984032  1.765971  9.093908  9.817983  5.660228  4.862172
##           [,87]    [,88]     [,89]      [,90]     [,91]     [,92]
## [1,] -0.8451055 -5.50187 -8.224718 -8.7787431 -0.909275 -9.489718
## [2,]  8.4774931 10.47169  4.194160 -0.4040521 11.577420  5.928886
##          [,93]      [,94]     [,95]     [,96]       [,97]      [,98]
## [1,] -8.676818 -12.213269 -1.869014 -8.071655 -11.1513349 -16.346495
## [2,]  8.846760   7.038192  3.500428  6.634393   0.4524286   7.099178
##          [,99]     [,100]    [,101]     [,102]    [,103]     [,104]
## [1,] -1.585486 -13.139405 -7.362648 -7.8511896 -12.32963 -6.9844242
## [2,]  3.192902   1.236548  5.226522  0.7617029 -12.32963  0.9768949
##         [,105]    [,106]     [,107]    [,108]    [,109]    [,110]
## [1,] -3.715153 -8.246160 -3.9581905 -4.067225 -5.333392 -9.209415
## [2,]  5.794236 -1.244913  0.4761632  1.634760  4.755950  4.667651
##         [,111]    [,112]
## [1,] -8.537175 -6.069659
## [2,]  7.084349  3.189483
## 
## $out
##   [1]  57.23895  40.26649 -27.70135  26.93377 -42.76788 -83.68659 -64.86258
##   [8] -56.43903  47.13970  50.20637  37.22124  38.16657 -41.43272 -39.45563
##  [15] -49.90122 -47.23328 -58.66581 -45.59469 -56.12678 -55.37398 -60.40710
##  [22] -54.04010 -43.11149 -66.65220  61.93806 -45.31429 -44.99196 -63.04769
##  [29]  38.96422  46.97302  36.01015  44.33973  65.58207  17.54333  50.38347
##  [36]  21.00299  30.94646  58.79029  52.41431  37.10547 -46.88383  57.11280
##  [43]  35.53809  30.75189 -44.44308  55.59392  53.50759  50.69041 -41.79776
##  [50]  38.14004  34.42069 -38.99826  37.93660 -35.93195  40.58071 -42.55470
##  [57] -38.73449 -39.90165  52.76446  38.66851  37.94345 -45.28221 -53.11982
##  [64]  35.95011 -41.03240 -51.32376 -26.46909  38.17169 -41.96388 -38.38198
##  [71]  27.65111 -26.50127  31.96379  28.80139  34.07594  40.69981 -36.70854
##  [78]  31.92250  44.44892 -50.14892  49.69028  47.38221 -54.27942 -63.65110
##  [85] -58.41295 -44.02494  39.86310  35.83621  34.48925  28.99213  40.85945
##  [92]  49.52424 -43.92605  45.01620 -41.05499 -55.81893 -54.00421  64.24987
##  [99]  58.85180  59.82918 -52.16430 -60.63002  55.77437  55.77437 -67.22244
## 
## $group
##   [1]   1   3   3   3   9   9   9   9   9   9  13  13  13  13  13  13  13
##  [18]  13  13  14  14  14  14  14  14  14  14  14  16  16  16  17  19  19
##  [35]  20  21  22  23  24  28  29  32  39  39  40  41  42  43  43  48  48
##  [52]  48  48  48  48  56  56  56  68  68  68  68  68  69  70  70  71  71
##  [69]  71  71  71  71  71  71  73  76  76  76  78  82  82  82  84  86  87
##  [86]  87  90  90  90  94  95  97  99  99  99 101 105 107 107 107 107 107
## [103] 107 107 108
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##             [,1]       [,2]       [,3]       [,4]        [,5]        [,6]
## [1,] -49.0941480 -43.470169 -16.408417 -31.244218 -15.9615500 -22.6494681
## [2,] -14.3491163 -14.803633  -8.333018 -18.006080  -5.4685949  -9.6004198
## [3,]  -0.8373246   2.080794  -1.350929  -4.788879   0.1652414  -0.3979381
## [4,]  14.8777774  16.568308   5.641032  14.012205   8.0096190  13.2698635
## [5,]  56.2460355  60.379473  14.695517  37.187011   9.9584058  39.1593411
##            [,7]       [,8]       [,9]       [,10]      [,11]      [,12]
## [1,] -25.063626 -27.395313 -43.496820 -38.0543298 -35.725398 -37.421415
## [2,] -13.043108 -12.443142 -11.575813 -14.7731511 -20.536526 -12.150124
## [3,]  -7.785661  -5.521222   1.431297  -0.9221918   2.242768  -1.411297
## [4,]   7.103369  -1.186577  12.316862  13.7935231  12.404465  13.709046
## [5,]  19.718648  -1.186577  46.227799  36.8070127  30.849222  46.152739
##           [,13]      [,14]      [,15]       [,16]      [,17]      [,18]
## [1,] -34.615857 -40.738412 -32.476346 -25.7904142 -37.326309 -29.484853
## [2,]  -9.484127  -8.712480 -12.373247  -5.1079727 -17.123527 -14.156916
## [3,]   1.081841   2.986876   5.133803   0.6211951  -6.303771  -1.702258
## [4,]  10.347297  13.399933  15.817944  10.4156434   9.067655  10.093937
## [5,]  36.345971  38.734234  30.288417  32.6251713  45.188784  25.429771
##           [,19]      [,20]      [,21]      [,22]      [,23]      [,24]
## [1,] -22.991024 -25.643198 -20.802287 -11.463394 -42.298876 -35.015235
## [2,] -17.492870 -21.139704 -17.299467  -4.021960 -15.166738 -13.203324
## [3,] -11.504156  -6.420577  -9.277762  -1.360295  -4.210784   4.142962
## [4,]   5.181761   7.590668   1.263290   4.710996  10.018621  12.033605
## [5,]  25.792152  35.678968  20.437952  14.965279  40.557989  49.218975
##            [,25]      [,26]      [,27]      [,28]      [,29]      [,30]
## [1,] -27.6825493 -36.322436 -31.845140 -27.822732 -30.043947 -32.249852
## [2,]  -7.7935727 -11.152331 -16.159881  -8.344541  -8.132391  -7.508550
## [3,]   0.4193811  -1.851963  -1.266264   3.760297   6.057817   4.067523
## [4,]  13.6452101  11.420739  16.612462  11.750841  12.866627  11.948712
## [5,]  37.3485698  30.782108  43.777504  37.035126  21.427700  31.786537
##           [,31]      [,32]      [,33]      [,34]     [,35]      [,36]
## [1,] -43.773320 -38.926741 -35.184466 -51.345564 -32.93158 -15.042653
## [2,] -21.431673 -15.137152 -10.326609 -15.451737 -18.26309 -15.042653
## [3,]  -3.967591  -2.234199   1.340058  -2.118404  -3.33550 -11.851089
## [4,]  15.474774  14.722904  15.384104  15.219932  11.73038  -8.659525
## [5,]  56.058982  38.009874  42.050770  47.313026  22.81371  -8.659525
##           [,37]      [,38]      [,39]        [,40]       [,41]      [,42]
## [1,] -19.952883 -29.191758 -23.130210 -30.42797663 -34.0857044 -21.324525
## [2,] -19.952883 -15.379325 -12.694837 -10.27119820 -13.5354826 -10.996042
## [3,] -12.125090   2.380884  -6.851635  -0.08385861   0.5601128  -2.133182
## [4,]  -4.297298  17.034732  13.553064   9.86630392  13.1311841   6.492825
## [5,]  -4.297298  38.425883  34.927783  39.38581894  44.5713017  16.526594
##           [,43]      [,44]     [,45]      [,46]      [,47]       [,48]
## [1,] -33.727289 -40.907143 -32.12426 -29.903821 -20.914745 -33.8265219
## [2,] -12.972126 -14.107270 -12.86659 -11.906322  -8.963793  -8.7051400
## [3,]  -3.691913  -6.484649   1.20113   1.904988  -1.118425  -0.5139491
## [4,]   9.242035  15.142974  11.49266   8.623204   7.962518  10.3967323
## [5,]  38.188597  34.536371  30.24750  22.189273  23.810607  38.6331956
##           [,49]      [,50]       [,51]      [,52]      [,53]      [,54]
## [1,] -27.352356 -40.273389 -30.3036585 -23.289059 -33.590893 -36.197488
## [2,]  -9.178777 -10.170291 -10.2625409 -11.480250  -7.402883 -11.052951
## [3,]  -3.112176  -1.117206   0.3060292  -3.573344   1.381276   4.275179
## [4,]   9.792288  12.775659  11.4041258  10.044274  14.588289  12.664378
## [5,]  31.814653  34.799640  31.5462683  30.612844  31.894344  34.601001
##           [,55]       [,56]       [,57]      [,58]      [,59]       [,60]
## [1,] -34.271614 -30.4406790 -16.4707366 -27.942119 -34.169403 -34.7520699
## [2,] -18.390307  -6.8658269  -1.4936765 -13.937117 -12.779209  -9.1319504
## [3,]  -4.386656   0.9751686   0.7057585  -4.608786  -2.644879   0.7202734
## [4,]  15.614695  10.0851249  10.8689379  13.866690  10.961274   9.9508315
## [5,]  45.217343  35.2156250  18.1557599  53.866690  34.306073  36.0733867
##           [,61]      [,62]      [,63]       [,64]      [,65]      [,66]
## [1,] -12.200641 -28.511574 -20.082090 -34.5211509 -48.465498 -32.836751
## [2,]  -6.841394 -11.361915 -11.219230 -16.4238913 -12.023935  -9.771883
## [3,]  -1.482147  -4.660785   1.442434  -0.9077239   1.164629   1.054284
## [4,]  15.124721   8.748726   8.033080  13.0558925  14.209727  14.407685
## [5,]  31.731589  22.973518  31.726720  40.7385761  39.752676  35.299273
##           [,67]      [,68]       [,69]       [,70]      [,71]      [,72]
## [1,] -34.984506 -30.635292 -22.2258684 -24.8656711 -29.696302 -20.547243
## [2,] -16.585178  -9.090440 -10.2358729 -10.0329305  -9.269217 -14.579395
## [3,]  -4.636477  -1.664427  -0.9465856  -0.3731693  -2.602550   2.786090
## [4,]  16.564718   7.095334   9.0509615   7.2494513  13.420747   9.220021
## [5,]  37.757780  28.051288  33.3488456  31.9423371  39.480302  18.039475
##           [,73]      [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -32.674865 -36.444228 -27.040261 -27.133416 -36.465323 -27.743873
## [2,] -11.424846 -12.568506 -15.310145 -10.151447 -10.893286 -16.243115
## [3,]  -2.955419  -3.620989  -8.333823  -3.022515  -1.713067  -2.021293
## [4,]   7.044581  12.936509  11.844928   5.520027   6.467843  15.630568
## [5,]  29.838858  48.205157  44.676601  25.520027  21.757958  38.963901
##           [,79]      [,80]       [,81]      [,82]        [,83]      [,84]
## [1,] -34.061788 -46.369547 -34.3328838 -50.503558 -39.67556447 -33.944083
## [2,]  -8.721488 -12.168278  -7.9104760 -14.081945 -12.65761565 -14.637146
## [3,]  -1.878744  -2.487964  -0.1466952  -2.841705   0.04991837   1.585443
## [4,]  10.178882  13.717543  12.4759254  12.687821  20.78536221  17.967825
## [5,]  31.950181  34.029077  30.3532961  49.819959  49.04809476  45.521443
##            [,85]      [,86]       [,87]      [,88]      [,89]      [,90]
## [1,] -31.2735027 -44.771848 -38.5698188 -39.979786 -34.622749 -18.563965
## [2,] -13.8817374 -11.772591 -10.0427935 -17.869977 -13.329496  -8.421822
## [3,]  -0.5112406   1.882477  -0.6629131   1.216872  -1.820084  -4.653855
## [4,]  13.2485203  14.201213  10.2839504  14.777005  14.137012   2.198297
## [5,]  42.3441157  47.744237  32.6531367  45.708742  37.186061  15.441038
##           [,91]       [,92]      [,93]      [,94]       [,95]       [,96]
## [1,] -29.663293 -29.7621646 -37.507802 -21.349417 -36.5610861 -38.8950658
## [2,] -14.521151 -16.7131163 -19.174468 -18.712016 -10.3839635 -11.1436525
## [3,]   4.200558   0.6252193   6.109817  -6.169851   0.5835576   0.6651567
## [4,]  12.527553  11.8811573  14.150702   2.825249   9.2527254  13.8334859
## [5,]  33.385755  49.6692653  35.541247   9.208859  34.9181870  47.3318234
##           [,97]      [,98]       [,99]     [,100]     [,101]     [,102]
## [1,] -33.441129 -18.866021 -33.4834387 -44.650312 -38.104010 -31.007543
## [2,] -11.920691 -13.925488  -8.5320510 -16.611215 -14.553032 -12.049746
## [3,]  -4.578270  -4.107432   0.3876521  -5.616218  -2.755851  -4.368568
## [4,]   8.031845  10.247487   8.3036159  11.192592  11.117245   5.769907
## [5,]  26.938495  25.628797  30.2526242  50.908307  49.599155  27.623887
##         [,103]      [,104]     [,105]     [,106]     [,107]     [,108]
## [1,] -12.04209 -39.9610204 -40.240737 -33.607244 -43.909140 -56.943203
## [2,] -12.04209 -11.8529505 -10.697451 -13.188098 -12.584351 -17.569442
## [3,] -12.04209  -0.8969965   1.515681  -4.034373  -3.010489  -1.818375
## [4,] -12.04209   9.7892052  13.307551   8.880139  12.091006  16.276054
## [5,] -12.04209  35.4853851  46.537785  38.122998  48.162977  62.274195
##          [,109]     [,110]      [,111]      [,112]
## [1,] -54.712392 -37.381061 -28.8881713 -44.8675236
## [2,] -16.909052  -9.459299 -12.4112704 -17.5156388
## [3,]  -2.196893   1.780940  -0.5027626  -0.9914101
## [4,]  21.512454  13.849658  10.1514999  16.6018332
## [5,]  64.019761  29.016177  32.9954149  53.3697508
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]       [,4]      [,5]      [,6]
## [1,] -3.882244 -3.327482 -6.870679 -16.100899 -5.741089 -9.431700
## [2,]  2.207595  7.489070  4.168821   6.523141  6.071572  8.635824
##           [,7]       [,8]      [,9]     [,10]     [,11]     [,12]
## [1,] -16.61411 -13.475082 -1.661338 -7.242404 -10.76892 -6.365995
## [2,]   1.04279   2.432638  4.523933  5.398021  15.25446  3.543401
##          [,13]     [,14]     [,15]     [,16]      [,17]     [,18]
## [1,] -0.895922 0.7683354 -3.137454 -3.784040 -13.859068 -9.523550
## [2,]  3.059603 5.2054168 13.405060  5.026431   1.251527  6.119033
##           [,19]      [,20]      [,21]     [,22]      [,23]     [,24]
## [1,] -20.754368 -16.098608 -20.363141 -4.922942 -11.137833 -1.868322
## [2,]  -2.253944   3.257454   1.807616  2.202351   2.716265 10.154247
##          [,25]     [,26]     [,27]     [,28]     [,29]     [,30]
## [1,] -5.075590 -9.455862 -9.352989 -1.026301 -0.860367 -1.056222
## [2,]  5.914352  5.751937  6.820460  8.546896 12.976002  9.191269
##           [,31]      [,32]     [,33]      [,34]     [,35]      [,36]
## [1,] -14.440769 -11.313781 -6.076640 -10.085382 -17.62403 -18.982503
## [2,]   6.505587   6.845382  8.756755   5.848575  10.95303  -4.719675
##           [,37]     [,38]      [,39]     [,40]     [,41]     [,42]
## [1,] -29.615960 -8.073173 -15.901498 -3.475589 -3.254460 -8.472491
## [2,]   5.365779 12.834941   2.198229  3.307872  4.374686  4.206127
##           [,43]      [,44]     [,45]     [,46]     [,47]     [,48]
## [1,] -10.099970 -14.529715 -5.602584 -4.858521 -4.827092 -3.867389
## [2,]   2.716144   1.560417  8.004843  8.668497  2.590241  2.839491
##          [,49]     [,50]     [,51]     [,52]     [,53]    [,54]      [,55]
## [1,] -7.851527 -7.334819 -5.175691 -9.090292 -4.182545 -1.57718 -13.600915
## [2,]  1.627174  5.100406  5.787750  1.943603  6.945097 10.12754   4.827604
##          [,56]     [,57]      [,58]     [,59]     [,60]     [,61]
## [1,] -2.076981 -6.676995 -13.394789 -8.175421 -3.631619 -21.51993
## [2,]  4.027318  8.088512   4.177217  2.885664  5.072166  18.55564
##           [,62]     [,63]     [,64]     [,65]     [,66]      [,67]
## [1,] -13.833383 -5.727310 -9.015923 -2.231020 -5.699243 -12.718398
## [2,]   4.511814  8.612179  7.200475  4.560279  7.807811   3.445444
##          [,68]     [,69]     [,70]     [,71]     [,72]     [,73]
## [1,] -4.578800 -4.158743 -3.407187 -7.351018 -5.054694 -7.305569
## [2,]  1.249947  2.265572  2.660849  2.145917 10.626874  1.394731
##           [,74]      [,75]      [,76]     [,77]      [,78]     [,79]
## [1,] -10.978350 -18.739818 -6.6342676 -6.349682 -13.010837 -5.260016
## [2,]   3.736371   2.072171  0.5892376  2.923547   8.968251  1.502528
##          [,80]     [,81]     [,82]     [,83]     [,84]     [,85]     [,86]
## [1,] -7.447769 -4.845082 -9.146855 -8.103456 -8.717728 -6.698387 -1.265041
## [2,]  2.471841  4.551691  3.463444  8.203293 11.888614  5.675906  5.029995
##          [,87]     [,88]     [,89]     [,90]     [,91]     [,92]     [,93]
## [1,] -5.347550 -6.938996 -8.769179 -7.490157 -2.732291 -6.802155 -4.023406
## [2,]  4.021724  9.372739  5.129010 -1.817553 11.133407  8.052593 16.243039
##           [,94]     [,95]     [,96]      [,97]      [,98]     [,99]
## [1,] -18.200877 -2.248713 -6.204623 -9.6263109 -15.132888 -2.010821
## [2,]   5.861174  3.415828  7.534936  0.4697717   6.918025  2.786125
##           [,100]    [,101]     [,102]    [,103]    [,104]    [,105]
## [1,] -12.0933449 -9.925744 -8.9972301 -12.04209 -5.013543 -3.340489
## [2,]   0.8609097  4.414042  0.2600938 -12.04209  3.219550  6.371852
##          [,106]     [,107]    [,108]    [,109]    [,110]    [,111]
## [1,] -7.4210344 -5.1764365 -4.880395 -7.958850 -4.273571 -8.103193
## [2,] -0.6477112 -0.8445405  1.243645  3.565063  7.835451  7.097667
##         [,112]
## [1,] -6.178479
## [2,]  4.195659
## 
## $out
##  [1]  66.45749  45.25463 -33.97804  29.67402  34.34046  28.12643 -83.18055
##  [8]  48.13655 -70.12297 -45.74712 -39.86029 -39.92063 -56.68642 -52.55334
## [15] -48.39942 -42.83506 -43.89056  50.67965 -69.26237 -45.08315 -47.42149
## [22] -61.26173 -45.69937 -53.98596 -28.69960  50.92886  63.00515 -18.26956
## [29]  20.60986  61.21564 -41.25062  61.77705  41.00864  49.13273 -40.68703
## [36] -53.95369  47.04584 -38.74359 -44.56322 -39.31301 -37.43001 -36.69243
## [43] -36.95309  30.53774 -58.40759  57.03818  62.34650 -53.38479 -55.76200
## [50]  54.43367 -40.74256  43.93223  33.08066  33.05856 -42.73061 -63.67507
## [57] -41.62926  35.32057 -60.01881 -48.20781  52.59452  46.13695  30.83462
## [64]  33.42693  36.19384  38.14483  47.33098 -41.29510  40.99423 -64.60352
## [71] -62.82326  40.94397  42.17681 -40.99010 -40.55255  20.86698  23.91651
## [78]  36.86246  27.75825  56.85449  42.19867 -47.31178  39.54337 -45.64120
## [85] -48.96051  49.22054 -40.80325  36.20377  59.66807  37.96958  35.74191
## [92] -52.40316  59.26241  54.05499  56.20106  56.20106  68.10779
## 
## $group
##  [1]   1   3   3   5   5   8   9  13  13  13  13  13  14  14  14  14  14
## [18]  14  14  14  14  14  14  14  16  18  19  22  22  23  30  32  40  40
## [35]  40  41  43  47  47  48  48  56  56  57  58  58  65  65  65  68  68
## [52]  68  68  68  69  70  70  70  71  71  71  73  76  76  77  77  77  79
## [69]  79  83  86  87  87  87  87  90  90  90  90  91  95  95  95  95  95
## [86]  97  99  99 100 102 102 105 107 107 107 107 108
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##           [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
## [1,] -46.47264 -36.844330 -11.867113 -17.386901 -21.771452 -26.027156
## [2,] -14.76330 -12.043730  -6.077989  -9.798931 -12.100538 -15.936296
## [3,]  -1.53332  -2.249466   1.414849  -4.268681  -3.866353  -3.912514
## [4,]  15.14061  11.243404   9.072168  10.373129  18.684095  13.870421
## [5,]  40.76913  40.039928  25.216768  32.468238  39.557056  45.059465
##            [,7]       [,8]       [,9]       [,10]     [,11]      [,12]
## [1,] -1.7564861 -30.613100 -35.762869 -38.9478975 -39.99718 -34.925117
## [2,] -1.7564861  -8.575534  -9.057046  -9.6255911 -17.37752 -11.677575
## [3,] -0.8341226  -7.670252   1.353660  -0.7546702  -3.39217   1.809253
## [4,]  8.3615377  31.605456  12.694594  13.8195661  14.00471  13.685756
## [5,] 10.4974539  39.587771  39.200910  30.6991146  23.86637  43.209694
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -42.490228 -41.249948 -31.454886 -22.270292 -40.203833 -21.890880
## [2,] -10.487644  -8.853899 -12.129919  -9.029003 -13.635110 -10.518309
## [3,]   1.391905   2.632226   4.536747   3.761760  -4.181101  -3.233861
## [4,]  11.568399  13.573558  15.139005   8.916544   4.857642   9.684975
## [5,]  43.660767  39.463613  31.026920  23.264364  28.102789  27.341298
##           [,19]      [,20]      [,21]     [,22]      [,23]      [,24]
## [1,] -17.671745 -25.756402 -22.594575 -8.253097 -45.401909 -34.916448
## [2,] -13.044398 -14.767749 -16.887103 -5.301694 -14.726876 -11.611450
## [3,]  -8.699766  -6.996353 -11.140790 -1.482713  -2.454899   3.138557
## [4,]  -3.630179   7.105198   8.433446  2.613901  15.449618  12.175982
## [5,]   4.939261  34.877428  20.464167 14.338416  60.245360  47.814628
##            [,25]      [,26]      [,27]      [,28]      [,29]        [,30]
## [1,] -26.3967903 -36.574346 -42.152031 -26.798747 -17.765579 -29.80694206
## [2,] -15.6965747 -15.297055 -15.783290  -8.979425  -4.946256 -13.52248582
## [3,]  -0.8891244  -1.449711   2.154279   2.900124   3.864700   0.05278131
## [4,]  23.1711014   8.123272  16.188927  11.711080  10.744248   8.47436161
## [5,]  39.1161052  32.639907  45.048213  37.988371  21.902094  33.51006902
##           [,31]      [,32]      [,33]      [,34]      [,35]      [,36]
## [1,] -52.504310 -25.980759 -37.807398 -31.342595 -26.599353 -11.642049
## [2,] -19.773234 -15.300002  -9.453733 -15.454680 -24.881158 -11.642049
## [3,]  -1.975783  -3.188002   1.449654  -3.575131  -2.823011  -9.899615
## [4,]   9.023939  12.301118  12.000052  14.863682  15.686166  -8.157181
## [5,]  50.672179  38.233448  43.277343  45.147578  21.306494  -8.157181
##           [,37]      [,38]      [,39]       [,40]      [,41]      [,42]
## [1,] -24.106889 -30.962892 -18.977523 -30.6569122 -39.031665 -19.007975
## [2,] -24.106889  -6.985926 -10.857072 -10.3995713 -13.210944 -13.826517
## [3,] -14.034074   6.349277  -3.764641   0.7604848   1.008162  -3.676785
## [4,]  -3.961259  14.283227  10.689789   9.8518812  12.994462   3.964201
## [5,]  -3.961259  29.252593  31.811182  37.6342052  45.047221  27.009026
##           [,43]      [,44]      [,45]      [,46]       [,47]     [,48]
## [1,] -26.124906 -37.581386 -50.149333 -29.360324 -33.2711339 -31.40888
## [2,] -11.221415 -12.469732 -15.242439 -11.115238 -10.0544970 -10.06838
## [3,]  -1.025969  -6.228829  -1.580003   1.672121  -0.6381386  -2.94015
## [4,]   6.378111   4.799193  12.256502   8.033235   8.0446109   7.67457
## [5,]  30.982448  21.398510  36.273378  22.944988  31.4195814  31.36907
##           [,49]     [,50]       [,51]      [,52]       [,53]      [,54]
## [1,] -42.059988 -38.67250 -37.6472655 -31.876537 -35.7258688 -29.361067
## [2,] -12.218680 -11.15430  -9.0426874 -10.343965  -7.7945303 -11.524499
## [3,]  -0.335693  -0.21453   0.2304202  -5.769729  -0.9609287   2.671133
## [4,]  11.692882  12.35552  11.7386198  14.796141  14.6353246  11.110917
## [5,]  25.957069  33.03056  31.0819475  40.766777  31.8294726  33.905441
##           [,55]      [,56]      [,57]      [,58]      [,59]        [,60]
## [1,] -22.481917 -22.812294 -6.1855134 -57.747801 -33.124014 -31.75862866
## [2,] -10.176605  -3.679712 -0.5414741 -17.072767 -12.272487 -15.68923557
## [3,]  -1.275363   2.025521  6.4573553  -2.887907  -1.203288  -0.09196199
## [4,]   9.658064  11.051789 17.9520324  15.658308   9.835318  11.45333778
## [5,]  26.490061  20.997500 37.9432654  55.658308  35.422201  33.93945205
##            [,61]        [,62]      [,63]      [,64]       [,65]      [,66]
## [1,] -21.5021142 -27.15949524 -19.703171 -40.117636 -40.0254915 -41.769076
## [2,] -11.4374058  -9.79627244 -12.434247 -15.922573 -12.2272669 -10.311365
## [3,]  -1.3726973  -0.05499052   1.714226  -5.795758   0.6323538   5.023588
## [4,]   0.7665791  11.37945067  11.016383  13.977317  12.9376659  16.935948
## [5,]   2.9058555  30.81044626  32.176377  46.448331  46.3202272  31.521894
##            [,67]       [,68]      [,69]       [,70]      [,71]      [,72]
## [1,] -39.7362182 -27.0330341 -32.133545 -32.1731863 -21.000809 -28.212072
## [2,] -14.1935484  -8.8973922 -10.718673  -9.7904803  -7.953836 -12.973999
## [3,]  -0.6717389   0.5003502  -1.536944   0.2095197  -2.201753   4.419560
## [4,]  14.2902806   8.6208017   9.707090   8.3299712   5.888742   8.552845
## [5,]  36.5522432  29.2230593  27.563038  31.5981989  20.128493  15.198311
##           [,73]      [,74]      [,75]      [,76]      [,77]      [,78]
## [1,] -23.766367 -31.943806 -34.104313 -36.539692 -32.154344 -38.243384
## [2,] -12.312582 -10.808279 -10.167716 -12.067813 -10.228457 -13.866969
## [3,]  -3.038714  -3.290786   4.246338  -1.687675  -1.517177  -3.675647
## [4,]   5.382106  12.282472   9.173061   6.779531   6.658938  12.030045
## [5,]  28.113181  34.794862  15.700123  33.321563  29.870529  47.392150
##           [,79]     [,80]       [,81]      [,82]       [,83]       [,84]
## [1,] -25.395352 -40.86827 -34.8809976 -31.251888 -48.4307411 -33.8787328
## [2,]  -9.206408 -11.88323  -9.4256735 -12.121069 -11.7146999 -13.9122218
## [3,]  -1.722212   1.61585  -0.2703733  -2.972596  -0.8134479   0.6591157
## [4,]  10.461138  12.68480  12.8136903  11.565252  18.1679485  28.8295371
## [5,]  31.063396  29.61226  22.2114327  46.850910  48.4082919  37.1492884
##          [,85]      [,86]      [,87]      [,88]      [,89]      [,90]
## [1,] -45.34647 -52.129790 -46.434636 -37.361901 -37.164038 -17.677387
## [2,] -11.03894 -12.631304 -12.864718 -18.113128 -11.525392  -8.617162
## [3,]  -1.27415   2.373997   1.724432   2.412592   2.410199  -6.223603
## [4,]  12.78608  14.710054  10.115122  12.359117  12.368887   3.049505
## [5,]  42.78608  39.632237  34.209277  48.901217  44.071657  15.832440
##           [,91]      [,92]      [,93]      [,94]        [,95]      [,96]
## [1,] -42.365698 -35.364666 -39.950207 -21.054781 -39.01049428 -27.419976
## [2,] -11.467848 -14.643884 -19.348785 -15.998545 -11.02689257 -10.140429
## [3,]   2.627278  -4.557137   6.057161  -5.845047   0.07079385   1.126239
## [4,]  13.487367  13.050804  15.100757   0.824767   8.59594662   9.167188
## [5,]  34.773055  51.597019  36.732194  15.576881  30.98950572  36.941379
##           [,97]      [,98]      [,99]     [,100]      [,101]      [,102]
## [1,] -34.367417 -31.580939 -35.547359 -29.717601 -43.7037574 -21.0904317
## [2,] -12.102201  -9.365094  -9.034297 -16.633537 -13.7020466 -12.1067528
## [3,]  -2.062105  -2.908552   1.292195  -9.152111  -0.7995561   0.7891168
## [4,]   9.303433   8.856239   8.796488  11.417304  13.6610844   8.3073107
## [5,]  32.052022  19.415225  31.543342  51.310420  40.7715979  22.0136158
##         [,103]     [,104]     [,105]     [,106]     [,107]    [,108]
## [1,] -12.30868 -35.023129 -41.699999 -34.584552 -43.804302 -51.99996
## [2,] -12.30868 -11.741449 -10.987176 -13.533268 -13.946389 -15.63443
## [3,] -12.30868  -1.492179   1.319112  -4.038354  -2.651908  -2.01875
## [4,] -12.30868  10.740357  10.819596  11.680431  12.074160  17.67371
## [5,] -12.30868  35.350981  41.737262  37.371454  51.044171  65.58862
##          [,109]      [,110]     [,111]     [,112]
## [1,] -58.263242 -37.4306558 -29.882892 -46.693469
## [2,] -15.689666  -9.4866983 -13.676907 -17.586683
## [3,]  -0.116984  -0.5874952  -1.136261  -1.347396
## [4,]  18.650586  11.7905926   9.502037  13.738940
## [5,]  51.942767  29.3087866  35.093135  56.372345
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]       [,4]       [,5]       [,6]
## [1,] -4.648773 -6.263984 -4.569463 -11.395444 -17.356584 -15.686167
## [2,]  1.582134  1.765053  7.399161   2.858082   9.623878   7.861139
##           [,7]      [,8]      [,9]     [,10]      [,11]     [,12]
## [1,] -5.267974 -36.06204 -1.461843 -5.941777 -15.788150 -3.050441
## [2,]  3.599729  20.72153  4.169163  4.432436   9.003809  6.668947
##           [,13]     [,14]     [,15]     [,16]     [,17]     [,18]
## [1,] -0.8077161 0.3820765 -3.463918 -1.330762 -9.515654 -9.749746
## [2,]  3.5915255 4.8823746 12.537412  8.854283  1.153452  3.282025
##           [,19]       [,20]      [,21]     [,22]      [,23]     [,24]
## [1,] -12.540337 -14.3644112 -26.261814 -4.711913 -10.754723 -2.527466
## [2,]  -4.859195   0.3717052   3.980233  1.746488   5.844925  8.804580
##           [,25]     [,26]     [,27]     [,28]    [,29]     [,30]
## [1,] -10.851292 -9.339016 -5.735011 -2.028230 -1.30458 -5.739722
## [2,]   9.073044  6.439593 10.043569  7.828477  9.03398  5.845284
##           [,31]      [,32]     [,33]      [,34]     [,35]     [,36]
## [1,] -10.147740 -11.580706 -4.739060 -11.450338 -22.14879 -13.79301
## [2,]   6.196175   5.204702  7.638367   4.300075  16.50277  -6.00622
##           [,37]      [,38]      [,39]     [,40]     [,41]      [,42]
## [1,] -36.541351 -0.5103686 -11.193660 -2.650438 -2.740429 -10.125508
## [2,]   8.473203 13.2089221   3.664377  4.171408  4.756753   2.771939
##          [,43]      [,44]     [,45]     [,46]     [,47]     [,48]
## [1,] -6.102855 -10.978520 -9.260654 -4.636396 -4.603772 -6.055024
## [2,]  4.050918  -1.479137  6.100649  7.980639  3.327495  0.174724
##          [,49]     [,50]     [,51]       [,52]     [,53]     [,54]
## [1,] -6.309278 -6.584935 -5.027302 -12.2133855 -6.635738 -2.914259
## [2,]  5.637892  6.155875  5.488143   0.6739279  4.713880  8.256525
##          [,55]      [,56]     [,57]      [,58]     [,59]     [,60]
## [1,] -6.649919 -0.6269995 -4.586668 -13.230927 -6.353485 -6.281917
## [2,]  4.099193  4.6780409 17.501379   7.455113  3.946908  6.097993
##           [,61]     [,62]     [,63]      [,64]     [,65]     [,66]
## [1,] -12.505336 -9.713380 -7.019013 -14.019504 -2.624961 -2.586782
## [2,]   9.759941  9.603399 10.447465   2.427989  3.889669 12.633958
##          [,67]     [,68]     [,69]     [,70]      [,71]     [,72]
## [1,] -7.616076 -2.653936 -4.938786 -2.971626 -5.0986741 -2.672519
## [2,]  6.272599  3.654636  1.864898  3.390666  0.6951674 11.511638
##          [,73]     [,74]     [,75]     [,76]     [,77]      [,78]
## [1,] -7.206388 -9.951710 -3.165169 -6.031360 -6.027273 -12.604530
## [2,]  1.128960  3.370139 11.657846  2.656009  2.992918   5.253236
##          [,79]     [,80]     [,81]     [,82]     [,83]     [,84]     [,85]
## [1,] -5.240732 -3.091464 -5.395806 -8.551495 -8.098816 -12.84728 -6.707524
## [2,]  1.796308  6.323163  4.855059  2.606302  6.471921  14.16551  4.159224
##           [,86]     [,87]     [,88]     [,89]     [,90]     [,91]
## [1,] -0.9392411 -3.571656 -5.199983 -3.635113 -9.339405 -3.768989
## [2,]  5.6872360  7.020519 10.025166  8.455511 -3.107801  9.023545
##           [,92]    [,93]      [,94]     [,95]     [,96]     [,97]
## [1,] -11.750843 -4.41795 -15.242791 -2.759479 -4.184180 -7.477783
## [2,]   2.636569 16.53227   3.552696  2.901067  6.436658  3.353574
##           [,98]     [,99]     [,100]    [,101]    [,102]    [,103]
## [1,] -11.219424 -1.248046 -15.686787 -8.442275 -4.513444 -12.30868
## [2,]   5.402321  3.832436  -2.617435  6.843163  6.091678 -12.30868
##         [,104]    [,105]     [,106]     [,107]     [,108]    [,109]
## [1,] -5.768436 -3.092361 -7.9077276 -4.9359338 -5.0321564 -5.266887
## [2,]  2.784077  5.730584 -0.1689797 -0.3678817  0.9946555  5.032919
##        [,110]    [,111]    [,112]
## [1,] -6.11428 -8.944254 -6.110005
## [2,]  4.93929  6.671731  3.415213
## 
## $out
##   [1]  70.06659  68.24708  62.78004  40.68747 -41.69548 -50.83947  45.85723
##   [8] -19.87694  31.81938  25.91364 -19.68336 -26.10679 -52.49867 -43.50605
##  [15] -84.40891 -47.80639  46.39918  47.50753 -49.08684 -47.58007 -54.44491
##  [22]  50.53344 -52.33638 -54.57456 -46.19141 -43.48248 -64.77680 -52.30574
##  [29] -48.34055 -45.36203 -64.61575 -37.31261  37.59363  32.62987  42.34022
##  [36] -44.19008  47.25173  64.03131 -30.70767  23.07027  31.74507 -19.77862
##  [43]  55.91102  62.18123 -41.93420  45.90546  30.72262  44.32207  34.62270
##  [50]  34.76091  47.03408 -48.49040  71.15990  47.22503  38.27683  37.40071
##  [57] -54.34215  52.64151 -48.40211  45.03628 -33.56309 -41.40289 -36.35167
##  [64] -48.00942 -35.89362 -49.36634  55.46396 -35.38864 -43.96157 -35.24504
##  [71] -36.91333 -36.65985  38.14057  37.07003 -28.77040  27.00978  30.83768
##  [78]  27.47677  62.21341 -43.05037  41.56207  35.19281  35.22411 -36.06819
##  [85]  46.92716 -56.61916  40.46114 -52.86976 -43.10436 -47.85543  49.15622
##  [92] -64.55688  35.79690 -39.45644  24.37865  37.96126  42.72280  39.96148
##  [99] -41.98172  55.15747  47.79291 -41.79325  48.78942  60.28240  46.10525
## [106] -53.96752  58.74162  56.37903  56.37903  51.17877  54.90026
## 
## $group
##   [1]   1   1   1   3   3   4   4   7   7   7   7   7   9   9   9   9   9
##  [18]  13  13  13  13  13  14  14  14  14  14  14  14  14  14  16  16  17
##  [35]  17  17  17  19  19  19  22  22  31  32  40  40  42  43  44  44  44
##  [52]  44  44  48  48  48  51  54  55  55  56  56  56  56  56  59  68  68
##  [69]  68  68  68  68  70  71  71  71  71  71  73  76  76  76  77  77  77
##  [86]  78  79  79  79  82  82  86  90  90  90  90  90  95  95  95  96  96
## [103]  97 100 105 105 107 107 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"

## $stats
##             [,1]        [,2]       [,3]       [,4]       [,5]       [,6]
## [1,] -51.1109594 -41.7460876 -35.485920 -25.045155 -11.423689 -21.169536
## [2,] -14.0104848 -14.4258355 -14.051713 -11.078409  -2.420533 -12.180313
## [3,]   0.2988843   0.7993225   2.615615  -1.958957   4.168650  -3.506937
## [4,]  13.7517791  12.0957970  18.199002   9.005416  10.419639  14.157794
## [5,]  47.1565631  37.3328940  43.526758  33.242714  21.978982  41.067319
##            [,7]       [,8]       [,9]       [,10]      [,11]      [,12]
## [1,] -21.169516 -10.985766 -43.514533 -38.0776373 -37.551630 -35.204645
## [2,]  -6.901423 -10.985766 -12.981610 -11.4783120 -20.861072 -11.732722
## [3,]   1.149138  -6.037687   1.820361  -0.2532469  -7.124223   3.963988
## [4,]  15.547656  -1.202044  12.787964  13.1687839  16.528777  15.322189
## [5,]  22.115238  -1.202044  46.271358  30.9149980  35.813147  33.410792
##           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
## [1,] -36.139041 -39.890566 -38.609470 -23.937358 -36.104189 -49.371233
## [2,]  -8.745702  -8.625719 -11.076781  -8.027722 -19.976701 -13.110501
## [3,]   2.492255   2.818516   5.257605   2.064746  -7.422838   5.481303
## [4,]  10.473134  12.558499  16.425850   7.230643  19.680893  13.224143
## [5,]  37.466765  39.743432  31.752870  29.378813  43.626387  30.228093
##           [,19]      [,20]      [,21]     [,22]      [,23]      [,24]
## [1,] -20.165830 -28.537036 -21.650905 -8.711320 -46.512451 -53.531308
## [2,] -13.112214 -11.625165 -17.401760 -3.113350 -16.351570 -11.719282
## [3,]  -6.697388  -3.126405 -10.150379  1.265725  -1.014029   3.584278
## [4,]  -3.403863  11.792522   4.065151  4.150363  15.955640  16.585330
## [5,]   8.523758  38.812869  19.517340 11.506418  57.314044  52.083752
##           [,25]      [,26]      [,27]      [,28]      [,29]       [,30]
## [1,] -40.404078 -35.831268 -51.560627 -40.231252 -16.315489 -26.5490424
## [2,]  -8.368218  -8.958679 -18.271787  -8.634871  -6.149349  -7.4755635
## [3,]   1.783716  -1.168810   6.412849   1.868099   3.271497   0.5077485
## [4,]  13.284452   9.334873  16.740058  12.985025  10.935300   9.2389840
## [5,]  36.797923  30.571327  46.315198  43.994304  17.189245  33.0633598
##          [,31]     [,32]      [,33]      [,34]      [,35]      [,36]
## [1,] -38.01346 -40.18988 -34.252673 -30.976423 -43.838597 -12.280469
## [2,] -15.22773  -9.69626 -16.921444 -16.082585 -14.897036 -12.280469
## [3,]  -7.72452  -6.19738   1.835131  -2.084689   6.534165 -10.327194
## [4,]  10.27841  17.71125  16.411889  14.753379  12.633519  -8.373919
## [5,]  44.56425  57.21897  43.078556  46.918467  30.422910  -8.373919
##           [,37]      [,38]      [,39]       [,40]       [,41]       [,42]
## [1,] -28.632610 -24.604942 -24.336052 -39.4160854 -41.0907786 -16.7072250
## [2,] -28.632610  -4.921412 -11.080110  -9.6679266 -13.2688343  -9.7421298
## [3,] -15.495050  -1.188940  -4.035525  -0.1663484  -0.7918471  -0.8513953
## [4,]  -2.357491  11.247755   8.486445  10.3320734  13.5692335   3.5550683
## [5,]  -2.357491  30.995174  36.135876  34.5818104  43.8257507  11.0938338
##            [,43]     [,44]      [,45]      [,46]      [,47]      [,48]
## [1,] -38.0759437 -48.11625 -62.261756 -29.730710 -25.694215 -38.612622
## [2,] -12.5516703 -14.34154 -17.407614 -11.649218  -8.403614 -10.671831
## [3,]  -0.9616218  -2.51515   1.277244   2.102098  -1.355366  -1.890415
## [4,]   7.7715775  10.05130  15.526981   8.688324   6.501143   8.109585
## [5,]  20.2531327  34.18075  47.505741  22.434379  20.542220  35.944497
##           [,49]      [,50]       [,51]      [,52]      [,53]      [,54]
## [1,] -27.879911 -39.082833 -30.8822768 -24.424464 -37.083494 -29.380356
## [2,] -11.711667 -11.751604  -9.9862972  -9.258324  -8.127000 -11.722796
## [3,]  -1.127544   1.115341  -0.2177143  -3.973447   1.323783   2.877668
## [4,]   9.895668  12.085411  11.2023717   9.241150  13.665552  15.081685
## [5,]  31.562467  42.459185  34.0799313  29.573431  37.485335  38.446791
##          [,55]       [,56]     [,57]      [,58]      [,59]      [,60]
## [1,] -41.20270 -22.8254041 -6.876031 -59.686451 -31.661198 -29.654053
## [2,] -13.18202  -4.8185541 -2.651281 -15.023992 -12.101433 -11.436161
## [3,]  -1.26351  -0.8206583  1.723106  -2.194341  -2.577601   2.861277
## [4,]  12.99149  10.5084666  9.020887  20.142148   8.510203  12.071405
## [5,]  41.98412  29.1411967 20.229044  52.971799  36.838276  33.738072
##           [,61]       [,62]      [,63]      [,64]     [,65]      [,66]
## [1,] -18.679898 -14.1027107 -19.395613 -28.357644 -45.20554 -28.371272
## [2,] -11.040002  -4.5637828 -11.892982 -16.890575 -10.67321  -6.929831
## [3,]  -3.400107   0.2676552  -2.215937  -6.020576   1.42856   2.715308
## [4,]  19.261907   7.7246144   8.107018  11.971914  13.65696  11.590379
## [5,]  41.923920  13.7606535  32.437195  46.856129  37.47897  31.072146
##           [,67]      [,68]      [,69]       [,70]      [,71]      [,72]
## [1,] -28.971579 -26.459681 -28.577006 -24.5202040 -31.475069 -20.368895
## [2,] -13.348738  -9.618794  -9.952611 -10.3130569  -8.269206 -11.853507
## [3,]  -3.812465  -1.955284  -1.238378  -0.8843668  -2.070444  -1.204858
## [4,]  15.780074   5.718747   8.194059   7.6713522   8.499856   8.886103
## [5,]  31.771657  28.630799  34.854598  29.9064457  29.924534  21.914819
##           [,73]      [,74]      [,75]      [,76]       [,77]      [,78]
## [1,] -22.573850 -36.826560 -24.582064 -24.216515 -27.2249231 -32.885979
## [2,] -10.610206 -11.149956 -14.224484  -9.717193 -10.5667831 -14.009381
## [3,]  -2.775294  -3.148016  -3.352877  -3.320873  -0.5419137  -4.429504
## [4,]   5.207642   9.274960   5.496611   5.721157   6.4379653   8.562864
## [5,]  25.038271  26.915855  26.160057  24.181758  30.6767570  37.823427
##           [,79]      [,80]       [,81]      [,82]       [,83]      [,84]
## [1,] -45.020743 -40.650554 -33.6976779 -46.447766 -53.4704488 -31.069163
## [2,] -12.163409 -10.736255  -8.8537080 -13.260346 -13.0298657 -13.303233
## [3,]  -1.241155  -1.593312   0.3002179  -3.924908  -0.7239007  -2.970951
## [4,]  14.342968  13.765323  12.2991658  10.270208  17.1415355  23.740805
## [5,]  43.655395  39.668678  23.3012700  29.740706  46.1341708  43.363434
##            [,85]      [,86]       [,87]       [,88]     [,89]      [,90]
## [1,] -36.8374855 -43.011652 -33.7194440 -45.1963072 -35.88561 -19.470095
## [2,] -12.7846081 -12.521993  -9.3325477 -17.9204921 -10.45592  -8.157958
## [3,]  -0.9386215   2.472947  -0.2064572   0.5789818   2.29646  -5.335089
## [4,]  14.1470790  13.374684   8.9261244  14.5016981  12.61386   3.164385
## [5,]  41.9767300  43.057786  32.3907979  45.3324012  35.94720  17.001401
##           [,91]      [,92]      [,93]      [,94]      [,95]        [,96]
## [1,] -28.274730 -28.634126 -40.451981 -20.155554 -36.867713 -31.14099400
## [2,] -13.108589 -14.351640 -18.665444 -18.801668 -10.120606 -10.35022865
## [3,]   1.377274  -3.183396   5.841394  -2.508664   1.219040   0.07709967
## [4,]  12.893515  12.818709  14.725275   4.674236   9.220092   8.09588979
## [5,]  34.726322  51.982745  36.002274   8.790275  33.132287  29.24113608
##           [,97]       [,98]       [,99]     [,100]      [,101]     [,102]
## [1,] -34.532618 -22.1983636 -39.0146482 -29.836830 -21.2317547 -22.801378
## [2,] -12.123765 -11.5312797  -9.7030472 -16.050837 -10.1506558 -12.125875
## [3,]  -5.861743  -3.4976762   0.3826534  -5.860196   0.4762775  -1.129031
## [4,]   7.163387   0.4992378   9.8634537   8.557305   9.8248027   6.980680
## [5,]  31.921041  16.9997639  37.9552854  44.332467  36.5101190  29.798984
##         [,103]    [,104]     [,105]     [,106]     [,107]    [,108]
## [1,] -11.17206 -43.79119 -46.687512 -33.347703 -48.875712 -56.72168
## [2,] -11.17206 -11.75158 -12.247769 -12.816795 -13.186246 -16.69629
## [3,] -11.17206  -3.58439   2.918372  -3.721218  -3.079567  -2.13734
## [4,] -11.17206  11.08859  12.415835   9.308178  12.993771  17.46239
## [5,] -11.17206  35.41877  44.751179  38.459129  47.342913  60.50333
##           [,109]     [,110]    [,111]     [,112]
## [1,] -56.6208218 -34.162349 -28.49722 -54.211647
## [2,] -14.7629537  -7.750430 -12.11223 -16.880332
## [3,]   0.5467251   0.167828  -1.50712  -1.123373
## [4,]  16.1893881  14.498005  10.64903  16.005399
## [5,]  51.6105549  31.829234  33.41700  55.742142
## 
## $n
##   [1] 230  84  16  20  13  16  13   5 149  51  16  68 251 248  29  31  30
##  [18]  24  15  22   7  15  33  44  38  22  41  44  23  36  31  27  30  37
##  [35]  11   2   2  24  21  88 122  19  30  33  32  23  52  81  40  34  39
##  [52]  38  39  41  34  77   7  25  46  48   3  12  18  33 149  32  42  77
##  [69]  90  81  57  23  45  30  17  47  35  21  78  68  47  45  42  25  48
##  [86] 170  47  40  39  35  38  37  27   8 120  33  39  12 123  46  32  37
## [103]   1  69  61 106 324 305 111  37  22 108
## 
## $conf
##           [,1]      [,2]      [,3]      [,4]      [,5]       [,6]
## [1,] -2.593447 -3.772798 -10.12342 -9.054547 -1.458082 -13.910490
## [2,]  3.191216  5.371443  15.35465  5.136633  9.795383   6.896615
##           [,7]        [,8]      [,9]     [,10]      [,11]     [,12]
## [1,] -8.688344 -12.9508400 -1.515218 -5.706275 -21.893213 -1.219818
## [2,] 10.986620   0.8754665  5.155940  5.199781   7.644767  9.147794
##          [,13]     [,14]     [,15]     [,16]      [,17]     [,18]
## [1,] 0.5755856 0.6931017 -2.811629 -2.265218 -18.862754 -3.012045
## [2,] 4.4089251 4.9439310 13.326839  6.394709   4.017078 13.974651
##           [,19]     [,20]      [,21]     [,22]     [,23]    [,24]
## [1,] -10.657951 -11.01482 -22.970071 -1.697537 -9.899891 -3.15771
## [2,]  -2.736825   4.76201   2.669313  4.228988  7.871834 10.32627
##          [,25]     [,26]     [,27]     [,28]     [,29]     [,30]
## [1,] -3.766076 -7.331123 -2.226483 -3.281630 -2.357088 -3.893749
## [2,]  7.333509  4.993503 15.052181  7.017828  8.900082  4.909246
##            [,31]      [,32]     [,33]      [,34]     [,35]      [,36]
## [1,] -14.9625598 -14.531211 -7.780442 -10.094343 -6.581058 -14.691704
## [2,]  -0.4864795   2.136452 11.450705   5.924965 19.649389  -5.962685
##          [,37]     [,38]      [,39]     [,40]     [,41]     [,42]
## [1,] -44.85037 -6.403757 -10.781765 -3.534920 -4.630938 -5.671322
## [2,]  13.86027  4.025878   2.710716  3.202223  3.047244  3.968531
##          [,43]     [,44]     [,45]     [,46]     [,47]     [,48]     [,49]
## [1,] -6.824212 -9.224223 -7.921625 -4.598162 -4.621096 -5.187597 -6.525487
## [2,]  4.900969  4.193923 10.476113  8.802357  1.910365  1.406767  4.270399
##          [,50]     [,51]      [,52]     [,53]     [,54]     [,55]
## [1,] -5.343722 -5.578500 -8.7150444 -4.189786 -3.736459 -8.355687
## [2,]  7.574404  5.143072  0.7681502  6.837353  9.491796  5.828666
##          [,56]     [,57]      [,58]     [,59]     [,60]     [,61]
## [1,] -3.580407 -5.247324 -13.306842 -7.379253 -2.499703 -31.04192
## [2,]  1.939090  8.693537   8.918159  2.224050  8.222256  24.24170
##          [,62]     [,63]      [,64]     [,65]     [,66]      [,67]
## [1,] -5.337165 -9.664128 -13.958993 -1.720705 -2.457520 -10.914049
## [2,]  5.872476  5.232254   1.917841  4.577825  7.888135   3.289119
##           [,68]     [,69]     [,70]     [,71]     [,72]      [,73]
## [1,] -4.7169267 -4.260645 -4.041630 -5.579808 -8.037580 -6.5009115
## [2,]  0.8063587  1.783888  2.272896  1.438920  5.627864  0.9503227
##          [,74]     [,75]      [,76]     [,77]      [,78]     [,79]
## [1,] -9.039935 -10.91012 -6.8788981 -5.083351 -12.212059 -5.983139
## [2,]  2.743902   4.20437  0.2371526  3.999523   3.353051  3.500830
##          [,80]     [,81]     [,82]     [,83]      [,84]     [,85]
## [1,] -6.287891 -4.574815 -9.467118 -8.079667 -14.676867 -7.080483
## [2,]  3.101268  5.175251  1.617302  6.631866   8.734964  5.203240
##           [,86]     [,87]     [,88]     [,89]     [,90]     [,91]
## [1,] -0.6652241 -4.414473 -7.520728 -3.540252 -8.358933 -5.287321
## [2,]  5.6111189  4.001558  8.678691  8.133173 -2.311245  8.041869
##           [,92]    [,93]     [,94]     [,95]    [,96]       [,97]
## [1,] -10.240905 -4.31176 -15.62264 -1.570539 -4.99637 -10.7414404
## [2,]   3.874113 15.99455  10.60531  4.008619  5.15057  -0.9820457
##          [,98]     [,99]      [,100]    [,101]    [,102]    [,103]
## [1,] -8.984876 -2.404865 -11.5928670 -5.103011 -6.091966 -11.17206
## [2,]  1.989524  3.170171  -0.1275255  6.055566  3.833904 -11.17206
##          [,104]    [,105]     [,106]     [,107]     [,108]    [,109]
## [1,] -7.9288110 -2.071033 -7.1165868 -5.3775906 -5.2276936 -4.095104
## [2,]  0.7600312  7.907776 -0.3258499 -0.7815431  0.9530135  5.188554
##         [,110]    [,111]    [,112]
## [1,] -5.611212 -9.174412 -6.123174
## [2,]  5.946868  6.160172  3.876428
## 
## $out
##   [1]  80.65101 -26.45781  53.41210 -22.72859 -32.28878  31.26456 -83.61185
##   [8] -43.09646 -38.94845  40.61370  43.70175 -37.60026 -39.35926 -40.46378
##  [15] -50.12358 -53.22215 -46.51090 -49.93636 -47.65463 -59.09113 -47.95714
##  [22] -52.77754 -51.26268 -42.90659 -41.85807 -48.18416 -49.65215 -45.70300
##  [29] -54.36161 -43.69088 -56.56224  32.47163  34.22983 -32.38101 -54.57167
##  [36] -38.72030  66.63594 -55.21491 -26.10492  51.15932  48.60748  51.24253
##  [43] -56.32018 -36.25483  44.53680  49.59187  34.79066 -39.48201  47.51555
##  [50] -30.15850 -40.98358 -43.07971 -26.31376 -58.36025 -58.53305  49.71664
##  [57]  43.94459  40.93448  29.72208  53.55505  50.64174  32.01667  36.22828
##  [64]  32.08279  42.00588  47.75793  47.73650  45.57141 -52.33216 -67.15458
##  [71]  52.62142 -55.27087  61.52515 -65.77291  54.78151 -58.91474 -38.04137
##  [78] -25.79033  25.50087  38.66281  26.07782  41.96404 -60.87585  39.24114
##  [85]  49.90570 -44.85011  48.29703  22.65027  44.29096  51.33141  49.16633
##  [92]  60.16317 -55.23303 -52.18851  43.64064 -52.92071  50.05297 -54.57374
##  [99]  59.66996  59.20040  55.96605  55.96605
## 
## $group
##   [1]   1   5   5   5   8   8   9  13  13  13  13  13  13  13  13  13  13
##  [18]  13  13  14  14  14  14  14  14  14  14  14  14  14  14  16  16  16
##  [35]  18  19  19  20  22  31  31  31  41  42  43  44  47  48  48  56  56
##  [52]  56  57  60  60  68  68  68  73  74  76  76  77  77  77  78  82  82
##  [69]  86  86  86  86  86  86  86  86  87  90  90  90  90  95  95  96  96
##  [86]  96  97  98  99 100 100 100 101 101 101 101 101 105 107 107 107 107
## 
## $names
##   [1] "17"  "18"  "57"  "58"  "59"  "60"  "61"  "70"  "77"  "129" "130"
##  [12] "132" "133" "135" "136" "138" "139" "141" "143" "145" "147" "150"
##  [23] "152" "153" "154" "155" "157" "159" "161" "162" "163" "164" "190"
##  [34] "191" "192" "197" "198" "202" "203" "233" "240" "244" "245" "246"
##  [45] "248" "251" "257" "260" "265" "266" "267" "268" "277" "278" "282"
##  [56] "305" "312" "313" "316" "333" "362" "363" "366" "375" "388" "389"
##  [67] "402" "406" "407" "408" "409" "414" "415" "419" "420" "424" "429"
##  [78] "430" "431" "432" "433" "434" "435" "436" "442" "445" "457" "458"
##  [89] "463" "467" "471" "476" "477" "505" "510" "511" "519" "527" "528"
## [100] "532" "536" "537" "546" "554" "573" "576" "603" "604" "605" "631"
## [111] "635" "693"
for (i in 1:10) {
  D1 <- lmer(mod3,data=MIdata[[i]])
#Assumption of Normality or residuals: want points to be near the line
print(qqmath(D1))
}

This produces a table of the pvalues for the correlations between the student level variables for pretest and the residuals.

res_corr <- data_frame(it =NA,
                       with_pre_score = NA,
                       with_cent_pre = NA)

for (i in 1:10) {
res_corr[i,1] <- i
D1 <- lmer(mod3,data=MIdata[[i]])
temp <- cor.test(resid(D1), MIdata[[i]]$pre_scor)
res_corr[i,2] <-  temp$p.value 
temp <-  cor.test(resid(D1), MIdata[[i]]$stud_pre_cent)
res_corr[i,3] <- temp$p.value
}
lev_p <- data_frame(it =NA,
                       p_values = NA)

for (i in 1:10) {
lev_p[i,1] <-i
D1 <- lmer(mod3,data=MIdata[[i]])
MIdata[[i]]$Model.F.Res<- residuals(D1) #extracts the residuals and places them in a new column in our original data table
MIdata[[i]]$Abs.Model.F.Res <-abs(MIdata[[i]]$Model.F.Res) #creates a new column with the absolute value of the residuals
Levene.Model.F <- lm(Model.F.Res ~ crse_id, data=MIdata[[i]]) #ANOVA of the residuals
temp <-anova(Levene.Model.F) #displays the results: want a p>0.05
lev_p[i,2] <- temp$`Pr(>F)`[[1]]
}

Making the figures

#linearity: Shouldn't see a pattern
D1 <- lmer(mod3,data=MIdata[[1]])
plot(D1, xlab="Fitted Value", ylab="Residual Variance")

#variables are not correlated to the residuals: want a p-value>0.05
cor.test(resid(D1), MIdata[[1]]$pre_scor) 
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$pre_scor
## t = 0.036379, df = 5957, p-value = 0.971
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02491982  0.02586189
## sample estimates:
##          cor 
## 0.0004713377
cor.test(resid(D1), MIdata[[1]]$stud_pre_cent)
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$stud_pre_cent
## t = 3.7072e-14, df = 5957, p-value = 1
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02539086  0.02539086
## sample estimates:
##          cor 
## 4.803239e-16
#quantitative homogeneity of variance
MIdata[[1]]$Model.F.Res <- residuals(D1) #extracts the residuals and places them in a new column in our original data table
MIdata[[1]]$Abs.Model.F.Res <-abs(MIdata[[1]]$Model.F.Res) #creates a new column with the absolute value of the residuals
#MIdata[[1]]$Model.F.Res2 <- MIdata[[1]]$Abs.Model.F.Res^2 #squares the absolute values of the residuals to provide the more robust estimate

Levene.Model.F <- lm(Model.F.Res ~ crse_id, data=MIdata[[1]]) #ANOVA of the residuals
anova(Levene.Model.F) #displays the results: want a p>0.05
## Analysis of Variance Table
## 
## Response: Model.F.Res
##             Df  Sum Sq Mean Sq F value Pr(>F)
## crse_id      1      47   46.98  0.1397 0.7085
## Residuals 5957 2002574  336.17
#Levene.Model.F2 <- lm(Model.F.Res2 ~ crse_id, data=MIdata[[1]]) #ANOVA of the squared residuals
#anova(Levene.Model.F2) #displays the results

#visual homogeneity of variance
boxplot(MIdata[[1]]$Model.F.Res ~ MIdata[[1]]$crse_id, xlab = "Course", ylab = "Residuals" )

#boxplot(MIdata[[1]]$Model.F.Res2 ~ MIdata[[1]]$crse_id)

#Assumption of Normality or residuals: want points to be near the line
qqmath(D1)

Linearity… check residuals following this guys solution? https://lib.dr.iastate.edu/cgi/viewcontent.cgi?article=4284&context=etd https://stackoverflow.com/questions/33859440/residual-plots-for-multiple-imputation-using-mice-package-in-r https://biologyforfun.wordpress.com/2014/04/16/checking-glm-model-assumptions-in-r/

Linearity (we pass)

plot(D1)

X variables are not correlated to the residuals (we pass)

cor.test(resid(D1), MIdata[[1]]$pre_scor) 
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$pre_scor
## t = 0.036379, df = 5957, p-value = 0.971
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02491982  0.02586189
## sample estimates:
##          cor 
## 0.0004713377
cor.test(resid(D1), MIdata[[1]]$stud_pre_cent)
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$stud_pre_cent
## t = 3.7072e-14, df = 5957, p-value = 1
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02539086  0.02539086
## sample estimates:
##          cor 
## 4.803239e-16
cor.test(resid(D1), MIdata[[1]]$used_las) 
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$used_las
## t = 4.5456e-13, df = 5957, p-value = 1
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02539086  0.02539086
## sample estimates:
##         cor 
## 5.88948e-15
cor.test(resid(D1), MIdata[[1]]$collabnla) 
## 
##  Pearson's product-moment correlation
## 
## data:  resid(D1) and MIdata[[1]]$collabnla
## t = 1.8433e-14, df = 5957, p-value = 1
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.02539086  0.02539086
## sample estimates:
##          cor 
## 2.388256e-16

https://ademos.people.uic.edu/Chapter18.html homogenety of variance we pass with variance, but not variance^2 (don’t pass?) Maybe not required for HLM: https://stats.stackexchange.com/questions/77891/checking-assumptions-lmer-lme-mixed-models-in-r https://stats.stackexchange.com/questions/255546/test-homogeneity-in-lmer-models

MIdata[[1]]$Model.F.Res<- residuals(D1) #extracts the residuals and places them in a new column in our original data table
MIdata[[1]]$Abs.Model.F.Res <-abs(MIdata[[1]]$Model.F.Res) #creates a new column with the absolute value of the residuals
MIdata[[1]]$Model.F.Res2 <- MIdata[[1]]$Abs.Model.F.Res^2 #squares the absolute values of the residuals to provide the more robust estimate

Levene.Model.F <- lm(Model.F.Res ~ crse_id, data=MIdata[[1]]) #ANOVA of the residuals
anova(Levene.Model.F) #displays the results
## Analysis of Variance Table
## 
## Response: Model.F.Res
##             Df  Sum Sq Mean Sq F value Pr(>F)
## crse_id      1      47   46.98  0.1397 0.7085
## Residuals 5957 2002574  336.17
Levene.Model.F2 <- lm(Model.F.Res2 ~ crse_id, data=MIdata[[1]]) #ANOVA of the squared residuals
anova(Levene.Model.F2) #displays the results
## Analysis of Variance Table
## 
## Response: Model.F.Res2
##             Df     Sum Sq Mean Sq F value   Pr(>F)   
## crse_id      1    2682790 2682790  10.834 0.001003 **
## Residuals 5957 1475167082  247636                    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
boxplot(MIdata[[1]]$Model.F.Res ~ MIdata[[1]]$crse_id)

boxplot(MIdata[[1]]$Model.F.Res2 ~ MIdata[[1]]$crse_id)

Assumption of Normality or residuals (we pass)

qqmath(D1) #id: identifies values that may be exerting undue influence on the model (i.e. outliers) 

Try using Dharma: https://cran.r-project.org/web/packages/DHARMa/vignettes/DHARMa.html

simulationOutput <- with(MIdata[[1]],{simulateResiduals(fittedModel = D1, n = 250)})
## Model family was recoginzed or set as continous, but duplicate values were detected in the simulation - changing to integer residuals (see help for details)
plot(simulationOutput)

testUniformity(simulationOutput = simulationOutput)

## 
##  One-sample Kolmogorov-Smirnov test
## 
## data:  simulationOutput$scaledResiduals
## D = 0.027775, p-value = 0.0002032
## alternative hypothesis: two-sided

HLMdiag version of diagnostics

http://r-statistics.co/Assumptions-of-Linear-Regression.html

Creating groups for final model

WM=c(1,0,0,0,0,0,0,0,0,0,0,0,0)
WMC=c(1,0,1,0,0,0,0,0,0,0,0,0,0)
WMLA=c(1,0,0,1,0,0,0,0,0,0,0,0,0)
BM=c(1,0,0,0,0,1,0,0,0,0,0,0,0)
BMC=c(1,0,1,0,0,1,0,0,0,0,1,0,0)
BMLA=c(1,0,0,1,0,1,0,0,0,1,0,0,0)
WF=c(1,0,0,0,0,0,1,0,0,0,0,0,0)
WFC=c(1,0,1,0,0,0,1,0,0,0,0,0,1)
WFLA=c(1,0,0,1,0,0,1,0,0,0,0,1,0)
BF=c(1,0,0,0,0,1,1,0,0,0,0,0,0)
BFC=c(1,0,1,0,0,1,1,0,0,0,1,0,1)
BFLA=c(1,0,0,1,0,1,1,0,0,1,0,1,0)

contrast.formulas <-rbind('White Male'=WM, 'White Male Collab'=WMC, 'White Male LA'=WMLA,
                          'Black Male'=BM, 'Black Male Collab'=BMC, 'Black Male LA'=BMLA,
                          'White Female'=WF,'White Female Collab'=WFC, 'White Female LA'=WFLA,
                          'Black Female'=BF,'Black Female Collab'=BFC, 'Black Female LA'=BFLA)

Centering analysis

The code below generates models to compare the three different methods for centering the prescores.

HLM1.1.1<-with(MIdata,{lmer(gain~1 +  stud_pre_grand + (1|crse_id))})
HLM1.2.1<-with(MIdata,{lmer(gain~1 +  stud_pre_grand + class_pre_cent + (1|crse_id))})
HLM3.1<-with(MIdata,{lmer(gain~1 + stud_pre_grand+  used_las + collabnla + (1|crse_id))})
HLM4.1<-with(MIdata,{lmer(gain~1 + stud_pre_grand+  used_las + collabnla + (1+ stud_pre_grand|crse_id))})
HLM5.1<-with(MIdata,{lmer(gain~1 + stud_pre_grand+  used_las + collabnla + class_pre_cent + (1|crse_id))})
HLM6.1<-with(MIdata,{lmer(gain~1 + stud_pre_grand+  used_las + collabnla + class_pre_cent + (1+ stud_pre_cent|crse_id))})
HLM1.1.2<-with(MIdata,{lmer(gain~1 +  pre_scor + (1|crse_id))})
HLM1.2.2<-with(MIdata,{lmer(gain~1 +  pre_scor + class_pre_cent + (1|crse_id))})
HLM3.2<-with(MIdata,{lmer(gain~1 + pre_scor+  used_las + collabnla + (1|crse_id))})
HLM4.2<-with(MIdata,{lmer(gain~1 + pre_scor+  used_las + collabnla + (1+ stud_pre_grand|crse_id))})
HLM5.2<-with(MIdata,{lmer(gain~1 + pre_scor+  used_las + collabnla + class_pre_cent + (1|crse_id))})
HLM6.2<-with(MIdata,{lmer(gain~1 + pre_scor+  used_las + collabnla + class_pre_cent + (1+ stud_pre_cent|crse_id))})
testEstimates(HLM1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    18.284     0.872    20.977   669.237     0.000     0.131     0.119 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   60.620 
## Residual~~Residual            409.799 
## ICC|crse_id                     0.129 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      18.246     0.874    20.871   634.548     0.000     0.135     0.122 
## stud_pre_cent     0.440     0.023    18.844    19.164     0.000     2.178     0.714 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   62.777 
## Residual~~Residual            343.477 
## ICC|crse_id                     0.154 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       18.259     0.873    20.925   743.465     0.000     0.124     0.112 
## stud_pre_cent      0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## class_pre_cent    -0.030     0.096    -0.313   266.451     0.754     0.225     0.190 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   63.060 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.155 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.1.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.1.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       18.363     0.971    18.919  1042.418     0.000     0.102     0.095 
## stud_pre_grand    -0.432     0.023   -18.598    19.053     0.000     2.198     0.716 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   82.562 
## Residual~~Residual            343.443 
## ICC|crse_id                     0.194 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.1.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.1.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    33.911     1.357    24.985    52.280     0.000     0.709     0.436 
## pre_scor       -0.432     0.023   -18.598    19.053     0.000     2.198     0.716 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   82.562 
## Residual~~Residual            343.443 
## ICC|crse_id                     0.194 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.2.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.2.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       18.620     0.869    21.423   851.208     0.000     0.115     0.105 
## stud_pre_grand    -0.440     0.023   -18.844    19.165     0.000     2.177     0.714 
## class_pre_cent    -0.470     0.098    -4.796   238.554     0.000     0.241     0.201 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   63.060 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.155 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM1.2.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM1.2.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       34.468     1.286    26.799    44.526     0.000     0.817     0.473 
## pre_scor          -0.440     0.023   -18.844    19.165     0.000     2.177     0.714 
## class_pre_cent    -0.470     0.098    -4.796   238.554     0.000     0.241     0.201 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   63.060 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.155 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    13.696     2.131     6.428   891.880     0.000     0.112     0.102 
## used_las        6.014     2.352     2.557  1215.587     0.011     0.094     0.088 
## collabnla       3.512     3.008     1.168   269.731     0.244     0.223     0.189 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   56.570 
## Residual~~Residual            409.830 
## ICC|crse_id                     0.121 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM3, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM3, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      13.760     2.144     6.419   812.980     0.000     0.118     0.107 
## stud_pre_cent     0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## used_las          5.936     2.363     2.512  1160.173     0.012     0.097     0.090 
## collabnla         3.256     3.023     1.077   247.880     0.283     0.235     0.197 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   58.827 
## Residual~~Residual            343.497 
## ICC|crse_id                     0.146 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM3.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM3.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       11.272     2.322     4.854  1624.496     0.000     0.080     0.076 
## stud_pre_grand    -0.432     0.023   -18.660    19.111     0.000     2.187     0.715 
## used_las           8.805     2.573     3.422  2010.566     0.001     0.072     0.068 
## collabnla          7.076     3.196     2.214   568.343     0.027     0.144     0.129 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   73.324 
## Residual~~Residual            343.509 
## ICC|crse_id                     0.176 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM3.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM3.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    26.846     2.461    10.911   393.338     0.000     0.178     0.156 
## pre_scor       -0.432     0.023   -18.660    19.111     0.000     2.187     0.715 
## used_las        8.805     2.573     3.422  2010.566     0.001     0.072     0.068 
## collabnla       7.076     3.196     2.214   568.343     0.027     0.144     0.129 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   73.324 
## Residual~~Residual            343.509 
## ICC|crse_id                     0.176 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM4, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM4, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)      13.347     2.052     6.505  1402.786     0.000     0.087     0.081 
## stud_pre_cent     0.426     0.025    16.833    33.212     0.000     1.086     0.547 
## used_las          6.662     2.289     2.911   951.492     0.004     0.108     0.099 
## collabnla         2.952     2.951     1.001   201.219     0.318     0.268     0.219 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           58.825 
## Intercept~~stud_pre_cent|crse_id        0.437 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.331 
## ICC|crse_id                             0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM4.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM4.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       11.298     2.276     4.963  2878.877     0.000     0.059     0.057 
## stud_pre_grand    -0.413     0.024   -17.007    32.604     0.000     1.107     0.552 
## used_las           9.155     2.543     3.601  1232.006     0.000     0.093     0.087 
## collabnla          6.780     3.104     2.184   442.351     0.029     0.166     0.146 
## 
##                                        Estimate 
## Intercept~~Intercept|crse_id             73.946 
## Intercept~~stud_pre_grand|crse_id        -0.584 
## stud_pre_grand~~stud_pre_grand|crse_id    0.009 
## Residual~~Residual                      340.496 
## ICC|crse_id                               0.178 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM4.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM4.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##              Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)    26.190     2.446    10.706  1066.612     0.000     0.101     0.094 
## pre_scor       -0.413     0.024   -17.007    32.604     0.000     1.107     0.552 
## used_las        9.155     2.543     3.601  1232.042     0.000     0.093     0.087 
## collabnla       6.780     3.104     2.184   442.355     0.029     0.166     0.146 
## 
##                                        Estimate 
## Intercept~~Intercept|crse_id             73.946 
## Intercept~~stud_pre_grand|crse_id        -0.584 
## stud_pre_grand~~stud_pre_grand|crse_id    0.009 
## Residual~~Residual                      340.496 
## ICC|crse_id                               0.178 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM5, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM5, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       13.659     2.209     6.184  1545.906     0.000     0.083     0.077 
## stud_pre_cent      0.440     0.023    18.844    19.165     0.000     2.177     0.714 
## used_las           6.036     2.445     2.469  1426.050     0.014     0.086     0.081 
## collabnla          3.388     3.135     1.081   306.195     0.281     0.207     0.177 
## class_pre_cent     0.015     0.098     0.158   337.055     0.874     0.195     0.168 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   59.394 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM5.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM5.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       14.021     2.213     6.335  1410.968     0.000     0.087     0.081 
## stud_pre_grand    -0.440     0.023   -18.844    19.165     0.000     2.177     0.714 
## used_las           6.036     2.445     2.469  1426.050     0.014     0.086     0.081 
## collabnla          3.388     3.135     1.081   306.195     0.281     0.207     0.177 
## class_pre_cent    -0.424     0.099    -4.270   285.200     0.000     0.216     0.183 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   59.394 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM5.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM5.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       29.869     2.393    12.481   302.391     0.000     0.208     0.178 
## pre_scor          -0.440     0.023   -18.844    19.165     0.000     2.177     0.714 
## used_las           6.036     2.445     2.469  1426.051     0.014     0.086     0.081 
## collabnla          3.388     3.135     1.081   306.195     0.281     0.207     0.177 
## class_pre_cent    -0.424     0.099    -4.270   285.200     0.000     0.216     0.183 
## 
##                              Estimate 
## Intercept~~Intercept|crse_id   59.394 
## Residual~~Residual            343.494 
## ICC|crse_id                     0.147 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM6, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM6, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       12.932     2.122     6.093  2448.080     0.000     0.065     0.061 
## stud_pre_cent      0.426     0.025    16.876    33.600     0.000     1.073     0.544 
## used_las           7.091     2.405     2.948   636.817     0.003     0.135     0.122 
## collabnla          3.493     3.032     1.152   279.923     0.250     0.218     0.185 
## class_pre_cent     0.060     0.096     0.623   218.282     0.534     0.255     0.210 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           59.644 
## Intercept~~stud_pre_cent|crse_id        0.455 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.323 
## ICC|crse_id                             0.149 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM6.1, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM6.1, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       13.281     2.132     6.229  1958.181     0.000     0.073     0.069 
## stud_pre_grand    -0.426     0.025   -16.876    33.600     0.000     1.073     0.544 
## used_las           7.091     2.405     2.948   636.863     0.003     0.135     0.122 
## collabnla          3.493     3.032     1.152   279.921     0.250     0.218     0.185 
## class_pre_cent    -0.366     0.097    -3.770   277.193     0.000     0.220     0.186 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           59.645 
## Intercept~~stud_pre_cent|crse_id        0.455 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.323 
## ICC|crse_id                             0.149 
## 
## Unadjusted hypothesis test as appropriate in larger samples.
testEstimates(HLM6.2, var.comp=TRUE)
## 
## Call:
## 
## testEstimates(model = HLM6.2, var.comp = TRUE)
## 
## Final parameter estimates and inferences obtained from 10 imputed data sets.
## 
##                 Estimate Std.Error   t.value        df   P(>|t|)       RIV       FMI 
## (Intercept)       28.620     2.319    12.343  1187.049     0.000     0.095     0.089 
## pre_scor          -0.426     0.025   -16.876    33.600     0.000     1.073     0.544 
## used_las           7.091     2.405     2.948   636.843     0.003     0.135     0.122 
## collabnla          3.493     3.032     1.152   279.912     0.250     0.218     0.185 
## class_pre_cent    -0.366     0.097    -3.770   277.173     0.000     0.220     0.186 
## 
##                                      Estimate 
## Intercept~~Intercept|crse_id           59.645 
## Intercept~~stud_pre_cent|crse_id        0.455 
## stud_pre_cent~~stud_pre_cent|crse_id    0.009 
## Residual~~Residual                    340.323 
## ICC|crse_id                             0.149 
## 
## Unadjusted hypothesis test as appropriate in larger samples.

This gets the citations for the important packages

citation(package = "hmi")
## 
## To cite package 'hmi' in publications use:
## 
##   Matthias Speidel, Joerg Drechsler and Shahab Jolani (2018). hmi:
##   Hierarchical Multiple Imputation. R package version 0.9.13.
##   https://CRAN.R-project.org/package=hmi
## 
## A BibTeX entry for LaTeX users is
## 
##   @Manual{,
##     title = {hmi: Hierarchical Multiple Imputation},
##     author = {Matthias Speidel and Joerg Drechsler and Shahab Jolani},
##     year = {2018},
##     note = {R package version 0.9.13},
##     url = {https://CRAN.R-project.org/package=hmi},
##   }
citation(package = "mitml")
## 
## To cite package 'mitml' in publications use:
## 
##   Simon Grund, Alexander Robitzsch and Oliver Luedtke (2018).
##   mitml: Tools for Multiple Imputation in Multilevel Modeling. R
##   package version 0.3-6. https://CRAN.R-project.org/package=mitml
## 
## A BibTeX entry for LaTeX users is
## 
##   @Manual{,
##     title = {mitml: Tools for Multiple Imputation in Multilevel Modeling},
##     author = {Simon Grund and Alexander Robitzsch and Oliver Luedtke},
##     year = {2018},
##     note = {R package version 0.3-6},
##     url = {https://CRAN.R-project.org/package=mitml},
##   }
## 
## ATTENTION: This citation information has been auto-generated from
## the package DESCRIPTION file and may need manual editing, see
## 'help("citation")'.
citation(package = "lme4")
## 
## To cite lme4 in publications use:
## 
##   Douglas Bates, Martin Maechler, Ben Bolker, Steve Walker (2015).
##   Fitting Linear Mixed-Effects Models Using lme4. Journal of
##   Statistical Software, 67(1), 1-48. doi:10.18637/jss.v067.i01.
## 
## A BibTeX entry for LaTeX users is
## 
##   @Article{,
##     title = {Fitting Linear Mixed-Effects Models Using {lme4}},
##     author = {Douglas Bates and Martin M{\"a}chler and Ben Bolker and Steve Walker},
##     journal = {Journal of Statistical Software},
##     year = {2015},
##     volume = {67},
##     number = {1},
##     pages = {1--48},
##     doi = {10.18637/jss.v067.i01},
##   }